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EXECUTIVE  SUMMARY 


This  study  was  conducted  to  Investigate  the  persistence  of  2,4,6- 
trinitrotoluene  (TNT)  and  hcxahydro-1 , 3, 5- trinl tro-1 , 3, 7-tr iazlne 
(RDX)  In  waste  disposal  lagoons  at  Louisiana  Army  Ammunition  Plant. 

The  waste  disposal  lagoons  are  complex  environments  due  to  a  hetero¬ 
geneous  distribution  of  solid  TNT  and  RDX  munitions  on  the  lagoon 
sediment  surface.  These  solids  serve  as  a  constant  input  source  of 
munitions  to  the  lagoon  waters  through  dissolution.  Flux  rates  for 
the  dissolution  of  munitions  were  measured  at  1.08  x  10-14  cm  sec-1  for 
TNT  and  7.75  x  10  5  cm  sec  1  for  RDX  at  25°C  in  the  lagoon. 

The  roajor  lose  processes  for  the  transformation  of  TNT  were 
phototransformation  and  biotransformation.  A  first-order,  depth- 
independent  rate  constant  for  the  indirect  photolysis  of  TNT  in  mid- 
March  was  measured  at  4.3  cm  d  1  which  yields  an  estimated  half-life 
of  eight  days  in  a  body  of  water  50  cm  deep.  A  major  stable  TNT 
phototransformation  product  was  2-araino-4 , 6-dinltrobenzoic  acid.  Bio- 
transformation  occurred  in  the  lagoon  with  a  pseudo-first-order  rate 
constant  of  O.OlOd  *  and  a  second-order  rate  constant  of  8.95  *  10  9 
ml  cell  1  d  The  half-life  of  TNT  resulting  from  blotransformation 
was  calculated  to  be  69  days.  The  major  blotransformation  products 
were  2-amino-4 , 6-dinitrotoluene  and  4-amino-2,6-dinitrotoluene. 

The  major  transformation  process  for  RDX  was  found  to  be  pho¬ 
tolysis.  RDX  transformed  with  first-order  rate  constants  ranging  from 
0.016  cm  d  1  in  winter  to  0.076  cm  d-*  in  summer.  The  half-life  of 
RDX  was  estimated  to  range  from  over  2000  days  in  winter  to  456  days 
in  summer  in  a  lagoon  50  cm  deep. 

Although  photolysis  was  found  to  be  the  major  transformation  pro¬ 
cess  for  TNT  and  RDX,  the  rate  of  transformation  was  slowed  due  to  the 
high  absorptivity  of  the  lagoon  waters,  whose  absorbances  ranged  from 
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3.0  in  winter  to  9.0  in  summer.  These  high  absorbances  allow  little 
light  to  penetrate  for  a  significant  distance  below  the  lagoon  sur- 
V  face.  The  light  absorption,  however,  heats  the  lagoon  water,  thus 

accelerating  the  dissolution  process  of  solid  TNT  and  RDX  on  the 
lagoon  sediment  surface.  Thus,  TNT  and  RDX  concentrations  in  the 
lagoon  water  were  found  to  increase  from  the  winter  to  the  summer 
months. 

A  computer  model  was  constructed  to  simulate  the  lagoon  environ¬ 
ment.  Using  the  laboratory-developed  transformation  and  dissolution 
rate  constants,  the  persistence  of  TNT  and  RDX  can  be  estimated  as  a 
function  of  time  from  initial  inputs  of  lagoon  depth  and  aqueous  and 
sediment-surface  munition  concentrations. 


CONTENTS 


EXECUTIVE  SUMMARY 


LIST  OF  ILLUSTRATIONS 


LIST  OF  TABLES, 


ACKNOWLEDGMENTS . 


I.  INTRODUCTION. 


II.  SITE  DESCRIPTION. 


III.  CHEMICALS  IDENTIFIED  IN  LAGOON  NO.  9. 


IV.  SITE  EVALUATION. 


V.  TRANSPORT  AND  TRANSFORMATION  STUDIES. 


Physical  Transport. 


1.  Water  Balance . 

2.  TNT  and  RDX  Solubility . 

3.  Sediment  Sorption . 

4.  Distribution  of  TNT  and  RDX 

in  Sediment  Cores . 

5.  Mass  Transfer  Rate  from  Sediment . 

6.  Determination  of  TNT  and  RDX  Flux  Values  for 

Sediment  from  LAAP  Lagoon  No.  9 . . 


Photochemistry. 


1.  TNT  Photochemical  Rate  Constant- 

General  Remarks . . . . 

2.  Effect  of  Depth  on  Photolysis  Rate  Constant 

of  TNT  In  LAAP  Water . 

3.  Experimental  Verification  of  Equation  27 . 

4.  Correction  of  Measured  Race  Constants 

for  Weather  Conditions..... . 

5.  Effect  of  Time  of  Year  on  Photochemical  Rate 

Constant  of  TNT  in  LAAP  Water . 

6.  UV-Visible  Spectra  of  LAAP  Water . 

7.  Experimental  Studies  of  the  Photolysis 

of  TNT  in  LAAF  Water . . . 

8.  RDX  Photolysis  Rate  Constant . 


CONTENTS  (Concluded) 


C.  Blotransforraatlon . 59 

1.  Bacterial  Count .  59 

2.  Biotransformation  Screening  Tests .  61 

a.  Aerobic  TNT  Blotransforraatlon .  62 

b.  Anaerobic  TNT  Biotransformation . 64 

c.  Aerobic  TNT  Blotransforraatlon  at 

High  and  Low  Concentrations  of  TNT......  66 

d.  TNT  Metabolites .  67 

e.  Biological  Oxygen  Demand 

of  the  Lagoon  Water . 70 

>  .  ^nx  Aerobic  Blotransf ormation .  70 

g.  RDX  Anaerobic  Blotransforraatlon .  71 

3.  Detailed  Biotransformation  Rate  Study . .  73 

VI.  DEVELOPMENT  OF  THE  SIMULATION  MODEL .  78 

A.  General  Equations . . . *  78 

B.  Program  Structure .  79 

C.  Simulation  Results.... .  81 

1.  Simulation  of  the  Water  Surface .  81 

2.  TNT  Simulation . 84 

a.  Case  Is  Complete  Mixing . 84 

b.  Case  IT:  Zero  Mixing . 90 

3.  RDX  Simulation .  93 

VII.  DISCUSSION .  96 

VIII.  LITERATURE  CITED .  144 

APPENDIX  A:  Analysis  of  LAAP  Lagoon  No.  9  Sediment.... .  100 

APPENDIX  B:  Mass  Transfer  Problems — Assumptions 

and  Solutions .  105 

APPENDIX  C:  SRI  Computer  Program  for  Modeling  the 

Fate  of  TNT  and  RDX  in  Lagoons  containing 

Munition  Wastes . 130 

DISTRIBUTION  LIST .  146 


iv 


ILLUSTRATIONS 


Lagoon  Area  at  Louisiana  Army 

Ammunition  Plant . 

Diagram  of  Lagoon  No.  9  and  TNT-RDX 
Distributions  (%)  in  the  Sediment . 

HPLC  Profile  of  Lagoon  Water 

(Sampled  7/22/82) . 

HPLC  Profile  of  the  Methylated  Ether  II 
Fraction  of  Lagoon  Water.. . . . 

HPLC  Profile  of  Lagoon  Water  Resolving  2-Amino 
4,6-Dinitro  Benzoic  Acid  and  Other  Components. 

Estimation  of  Lagoon  Depth  Profile  Based  on 
Variable  Seepage  Rates . 

Lagoon  Water  Rests  on  Sediment . . 

TNT  (•)  and  RDX  (o)  Flux  for  LAAP  Lagoon 
No.  9  Sediment . 

Absorption  Spectrum  for  Lagoon  9,  Site  1  Water 
Filtered,  0.1-cm  Cell.. . . 

Aerobic  Biotransformation  of  TNT  and  RDX 
in  Lagoon  Water.... . . . . 

Anaerobic  Blotransforraation  of  TNT  and  RDX  in 
Lagoon  Water..... . . . 

Radioautograph  of  Thin-layer  Chromatogram  of 
TNT  Metabolites  in  Lagoon  Water  after  19  Days 
of  Incubation . 

Lineweaver-Burke  Plot  of  TNT  Biotransformation 
Rate  in  Lagoon  Water.... . . 

Simulation  Flow  Chart  of  the  Lagoon  Model . 

Program  Interactions. . . . 


ILLUSTRATIONS  (Concluded) 


6.  Approximate  Water  Contour  Lines  During 
the  Simulation . 


7.  Simulated  TNT  Concentration  with  Respect  to 
Time  of  Year . . . 


8.  Simulated  TNT  Concentration  as  a  Function  of 
Time  of  Year . . . . . . 


9.  Simulated  TNT  Concentration  from  the  Case 
of  Sole  Input  from  Dissolution . . . 


Concentration-Time-Dependent  Profile  of 
TNT  Measured  in  Lagoon  No.  9 . . 


Simulation  Results  for  TNT  in  Compartments  1 
and  13  for  the  Case  of  Zero-Mixing . 


Comparison  of  the  Simulation  Results  from  the 
Case  with  Dissolution  As  the  Only  Input  Source 
of  TNT  to  Measured  Data  from  the  Lagoon . 


Simulation  Results  for  RDX  Concentration 
in  LAAP  Lagoon  Water . 


HPLC  Profile  of  Lagoon  No.  9 
Sediment  Extract . 


Lagoon  Water  Rests  on  Sediment. 


Boundary  Layer  Approximation  in  Sediment/ 
Water  Phase . . . . 


TNT  (•)  and  RDX  (o)  Flux  Data  Compared  with  Curve 
Fit  (Equation  62) . . . 


TABLES 


Lagoon  Water  Temperatures  (°C)  for  the  Bottom, 
Middle,  and  Surface  on  December  9,  1981... . 

Concentrations  of  TNT  and  RDX  in  Water  at 
Lagoon  Center  and  Southeast  Edge  As  a  Function 
of  Time  of  Day . 

Concentrations  (ppm)  of  Munitions  Found  In 
Water  Samples  on  Various  Dates . 

Levels  of  2-Amino-4 , 6-dlnitrobenzolc  Acid 

Observed  In  Lagoon  No.  9  on  Various 

Sampling  Dates . . . 

Absorbance  Values  of  Lagoon  No.  9  Water  at 
290  and  400  ran  at  Various  Times  of  the  Year.... 

Average  Monthly  Precipitation  and  Evaporation 
Data  from  1979  through  1981  at  Red  River 
Weather  Station,  Restin,  LA . . . 

Monthly  Precipitation  and  Evaporation  Data 

from  January  1982  through  August  1982  at 

Red  River  Weather  Station,  Restin,  LA . 

Solubility  of  TNT  and  ROX  in  Water  at  10,  20, 
and  30°C . 

Measured  and  Calculated  Solubilities 

(mg  liter-1)  of  TNT  and  RDX  at  283,  293, 

and  303°K . 

Sorption  Partition  Coefficients  for  TNT  and 
RDX  in  LAAP  Lagoon  No.  9  Sediment . 

Sediment  Core  Sample  Dry  Weight  and  Percent  TNT 
and  RDX  Weights..... . . . 

TNT  and  RDX  in  Sectioned  Cores  from 

LAAP  Lagoon  No.  9 . 

TNT  and  RDX  Flux  Measurements  from  LAAP  Lagoon 
No.  9  Sediment . . . . . . 


vii 


TABLES  (Continued) 


14.  Depth  for  99%  Absorption  of  Light  of  Specific 

Wavelengths  for  LAAP  Water .  3° 

15.  Photolysis  of  4.97  ppm  TNT  In  LAAP  Water  In  an 

MGR  Apparatus  and  In  Sunlight .  43 

16.  Photolysis  of  4.97  ppm  TNT  in  LAAP  Water  As  a 

Function  of  Depth  Using  a  Xenon  Lamp . . .  45 

17.  Sunlight  Photolysis  of  4.97  ppm  TNT  in  LAAP  Water 

and  Actlnometers . . . . . .  46 

18.  Variation  in  TNT  Photolysis  Rates  in  LAAP  Water 

As  a  Function  of  Time  of  Year.... .  50 

19.  Absorption  Coefficient  of  December  and  July 

LAAP  Waters  at  Selected  Wavelengths .  51 

20.  Sunlight  Photolysis  of  TNT  in  LAAP  Water . .  55 

21.  Photochemical  Rate  Constants  for  RDX  in  July 

LAAP  Water  as  a  Function  of  Season . . .  58 

22.  Total  Aerobic  and  Anerobic  Bacteria  in  Lagoon  9 

Water  Collected  in  December .  60 

23.  Total  Aerobic  and  Anaerobic  Bacteria  in  Lagoon  9 

Sediment  Collected  in  December . . .  60 

24.  Bacterial  Plate  Count  of  LAAP  Lagoon  Water 

Collected  in  July....... . . . . .  61 

25.  RDX  Concentration  on  Different  Days  after 

Anaerobic  Incubation  in  LAAP  Lagoon  Water .  71 

26.  Anaerobic  Transformation  of  RDX  with  Selected 

Organisms  and  Yeast  Extract  Concentrations. .  72 

27.  Transformation  Rates  Determined  for  Various 

TNT  Concentrations*.. . .  75 

28.  Initial  TNT  Concentrations  and  Biotransformation 

Rates  over  5  Days. . . .  77 


viii 


TABLES  (Concluded) 


A-l.  Percent  Composition  of  Munition  Components  In 

LAAP  Lagoon  No.  9  Sediment . . . .  103 

A-2.  Percent  Water  in  Lagoon  No.  9  Sediment  Samples .  104 

B-l.  Calculated  Flux  Values  Considering  All  Parameters 
(Nj^),  No  Chemical  Reaction  (N2),  No  Convection 
(Nj),  No  Reaction  or  Convection  (N^),  and  Their 
Ratios  to  (N*) . . . . . .  113 

B-2.  Flux  Calculations  for  Sample  Lagoon  Problem .  120 

B-3.  Total  Mass  per  Unit  Area  Transferred  During  a 

Time  Interval  for  Sample  Lagoon  Problem... .  120 

B-4.  Parameters  for  Analyzing  Flux  Data . 128 


lx 


ACKNOWLEDGMENTS 


The  studies  reported  herein  could  not  have  been  performed  without 
the  assistance  of  an  excellent  support  staff.  We  therefore  gratefully 
acknowledge  the  help  of  Doris  Tse  (photochemistry),  Daniel  Haynes 
(physical  transport),  Philip  Alferness  (analytical  chemistry),  and 
Justine  Whaley  (microbiology).  We  also  want  to  thank  Dr.  John 
Winterle  for  his  contribution  to  the  development  of  equations 
describing  sensitized  photochemical  processes  in  LAAP  water  and  the 
observed  effect  of  depth  on  the  photochemical  rate  constant. 


x 


I.  INTRODUCTION 


The  U.S.  Army  Is  currently  involved  in  assessing  the  overall 
hazard  of  munitions-related  warer  pollutants  to  aquatic  and  mammalian 
systems.  Part  of  this  assessment  involves  the  loss  and  movement 
(fate)  of  militarily  unique  chemicals  when  they  are  discharged  into 
various  aquatic  environments.  Munitions  have  been  produced  and 
handled  for  many  years  with  minimal  pollution-abatement  controls,  but 
it  has  now  become  important  to  estimate  the  environmental  persistence 
of  selected  chemicals  used  in  present  and  past  military  operations. 

Environmental  persistence  can  be  estimated  through  the  use  of 
computer  modeling,  in  which  the  loss  and  movement  processes  as  well  as 
critical  environmental  parameters  can  be  described  by  mathematical 
equations  and  approximations.  In  most  environmental  situations  or  in 
cases  of  low  chemical  concentrations,  simple  first-order  or  pseudo- 
first-order  rate  expressions  and  equilibrium  constants  adequately 
describe  the  loss  and  movement  processes.  Once  the  rate  and  equi¬ 
librium  constants  are  known  for  the  dominant  processes,  environments 
can  be  simulated  under  any  number  of  specified  conditions.  Thus,  it 
becomes  possible  to  estimate  the  concentration  of  a  chemical  as  a 
function  of  tivae  when  the  chemical  is  subjected  to  various  environ¬ 
mental  stresses.  In  a  previous  study  (Spanggord  et  al. ,  1981),  we 
applied  the  modeling  concept  to  raunltions-related  chemicals  in  the 
Holston  River  and  generated  concentration-  time-dependent  profiles 
that  allowed  the  prediction  of  chemical  concentrations  20  km  from  the 
discharge  point.  Actual  data  on  samples  taken  from  the  river  verified 
the  applicability  of  the  modeling  concept. 

On  the  basis  of  that  experience,  this  study  was  undertaken  to 
estimate  the  persistence  of  2,4,6-trinitrotoluene  (TNT)  and  hexahydro- 
1 , 3, 5-trinitro-l , 3, 5-triazine  (RDX)  in  disposal  lagoon  environments. 
The  lagoon  disposal  of  munitions  wastes  has  been  practiced  for  many 
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years  and  the  Army  would  like  to  dispose  of  these  lagoons  in  a  safe 
manner.  The  studies  reported  herein  describe  the  major  transport  and 
transformation  processes  that  affect  the  persistence  of  TNT  and  RDX  In 
lagoons.  We  also  developed  a  computerized  simulation  model  to  des¬ 
cribe  the  loss  and  movement  of  these  chemicals  as  a  function  of  time, 
concentration,  and  location  In  the  lagoon. 
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II.  SITE  DESCRIPTION 


This  Investigation  covered  evaluation  of  the  persistence  of  TNT 
and  RDX  in  the  waste  disposal  lagoons  at  Louisiana  Army  Ammunition 
Plant  (LAAP)  in  Minden,  Louisiana.  LAAP  is  a  load,  assemble,  and  pack 
facility  that  generates  munitions  wastewaters  during  shell-washout 
operations  in  the  loading/unloading  process.  The  wastewaters  are 
transported  by  truck  to  disposal  lagoons,  where  natural  evaporation 
concentrates  the  residues. 

We  conducted  a  field  study  at  LAAP  on  December  7-11,  1981.  At 
that  time,  no  wastes  had  been  transported  to  the  lagoons  since  May  and 
a  new  carbon-column  waste-treatment  facility  was  to  be  Installed  prior 
to  further  waste  disposal. 

The  lagoon  disposal  area  at  LAAP  is  outlined  In  Figure  1.  Lagoon 
No.  9  was  selected  for  our  investigations  mainly  because  it  is  one  of 
the  larger  lagoons  and  because  previous  reports  had  indicated  that  TNT 

and  RDX  were  still  being  found  in  the  lagoon  waters  and  sediments 

(Wentzel,  et  al.,  1981). 

Lagoon  No.  9  is  a  nearly  square,  flat-bottomed  reservoir  mea¬ 
suring  56.7  m  x  52.1  m  (186  ft  *  171  ft).  During  the  week  of  December 

7-11,  1981,  it  was  filled  to  a  depth  of  60  cm  (2.0  ft)  with  a  dark 

tea-colored  water.  Water  input  into  the  lagoon  is  from  seasonal 
rains,  and  water  output  is  due  to  seepage  and  evaporation.  Trees  and 
shrubbery  line  the  lagoon  enbankments,  but  no  plant  or  aquatic  life 
was  evident  in  the  lagoon  water. 

To  bc*in  our  field  study,  a  tow  line  was  secured  across  the 
middle  of  the  lagoon  and  the  center  position  was  tagged.  We  used  an 
Inflatable  raft  to  transport  monitoring  and  sampling  equipment  to  the 
center  of  the  lagoon  for  sample  collection  and  monitoring  measure* 


ments 


Water  samples  were  collected  with  a  variable-water-depth  sample 
collector  in  50-ral  amber  glass  jars  with  Teflon-lined  caps.  We 
analyzed  these  samples  on  site,  using  a  high-pressure  liquid  chro¬ 
matograph  (HPLC).  We  collected  other  water  samples  in  4-liter  amber 
glass  bottles  and  transported  them  by  air  overnight  to  SRI  for  total 
microbial  counts.  Sediment  samples,  removed  from  the  lagoon  with  a 
hand  corer,  were  collected  at  various  locations  and  transported  to  SRI 
for  analysis.  These  samples  indicated  a  heavy  clay  bottom  that  was 
difficult  to  penetrate.  Except  for  sorption  properties,  we  did  not 
evaluate  the  morphological  characteristics  of  the  sediment. 

To  determine  the  ability  of  light  to  penetrate  into  the  lagoon 
water,  we  placed  photochemical  "trees"  at  various  depths  in  the 
lagoon.  These  trees  consisted  of  aluminum  rods,  cross  bars,  and  flo¬ 
tation  supports  and  were  equipped  with  sealed  PNA/PYR*  actlnometer 
solutions  attached  to  the  ends  of  the  crossbars.  After  selected  time 
intervals,  the  PNA/PYR  tubes  were  removed  from  the  apparatus,  wrapped 
in  aluminum  foil,  and  transported  to  SRI  for  analysis. 

Meteorological  conditions  that  were  measured  periodically  at  the 
site  were  temperature,  pressure,  wet  and  dry  bulb  temperatures  (rela¬ 
tive  humidity),  and  wind  velocity.  In  the  lagoon  water,  we  measured 
temperature,  pH,  and  dissolved  oxygen.  LAAP  personnel  obtained  a 
water  sample  in  March  1982.  During  a  second  visit  to  the  site  on 
May  25,  1982,  we  collected  water  and  sediment  samples.  In  addition, 
sample  bottles  were  left  with  plant  personnel  to  collect  samples  at 
two-week  intervals  and  to  mail  them  to  SRI.  Depth  at  the  center  of 
the  lagoon  was  measured  during  the  December  visit  and  again  on  July 
24,  1982,  to  estimate  water  loss  over  an  8-month  period.  The  lagoon 
was  completely  dry  in  late  August. 


*PNA/PYR,  jr-nitroacetophenone/pyridine,  sunlight  actlnometer  developed 
at  SRI. 


The  yearly  meteorological  conditions  for  the  lagoon  site  were 
estimated  from  dally  data  collected  by  the  Red  River  experimental 
weather  station  at  Restln,  LA,  approximately  30  miles  east  of  LAAP. 

The  data  collected  at  this  site  include  dally  maximum  and  minimum  air 
and  soil  temperatures,  24-hour  pan-evaporation,  and  24-hour  precip¬ 
itation.  We  averaged  the  data  collected  from  this  station  for  three 
years  (1979  through  1981)  to  estimate  water  balances  within  the  lagoon 
at  LAAP. 

A  graphic  sketch  of  Lagoon  No.  9  is  shown  in  Figure  2.  A 
concrete  ramp  in  the  southeast  corner  of  the  lagoon  serves  as  an 
access  point  for  vehicles  to  dump  liquid  and  solid  wastes.  Analysis 
of  sediment  samples  (Appendix  A)  collected  on  May  25,  1982,  from  a 
5x5  sectioned  matrix  of  the  lagoon  showed  very  high  concentrations 
of  munitions  that  decreased  as  a  function  of  distance  from  the  dump 
ramp  (Figure  2).  When  the  lagoon  dried  up  in  late  August  1982,  pure 
munitions  compounds  were  visible  on  the  sediment  surface  near  the 
ramp.  These  solids  serve  as  the  primary  input  source  of  munitions  to 
the  lagoon  waters.  The  concentration  of  munitions  also  decreased  as  a 
function  of  depth  in  the  sediment.  The  distribution  of  munitions  in 
small  sediment  cores  taken  from  various  locations  in  the  lagoon  is 
described  in  Section  V.A.4. 

Because  the  lagoon  is  a  highly  colored  body  of  water,  it  absorbs 
much  of  the  sunlight  energy  and  converts  it  to  heat,  thus  establishing 
a  thermocline  in  the  lagoon  water.  Temperature  fluctuations  during 
the  course  of  a  day  are  shown  in  Table  1  for  the  bottom,  middle,  and 
surface  of  the  lagoon.  Over  the  nine-month  course  of  this  study,  the 
lagoon  temperatures  ranged  from  10°C  in  December  to  354C  in  May  and 
41.5°C  in  August.  This  temperature  change  affects  rates  of 
dissolution,  diffusion,  and  microbial  activity  as  well  as  solubility. 

The  lagoor.  waters  appeared  to  be  well  mixed  in  a  vertical 
direction,  as  evidenced  by  similar  concentrations  of  TNT,  RDX,  and 
dissolved  oxygen  measured  at  the  surface  and  at  depths  of  25  and  50 


i 


cm.  Horizontal  mixing,  however,  was  not  complete;  concentrations  of 
TNT  and  RDX  were  higher  in  the  south  bank  waters  than  those  in  the 
center  of  the  lagoon  (Table  2).  The  higher  concentrations  near  the 
edge  of  the  lagoon  may  be  the  result  of  the  water  there  being  warmer 
which  increases  the  rates  of  dissolution  and  diffusion  of  munitions 
from  the  near-shore  sediments  relative  to  bulk  diffusion  rates  in  th 


main  lagoon. 


Table  1 
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LAGOON  WATER  TEMPERATURES  (°C)  FOR  THE  BOTTOM,  MIDDLE, 
AND  SURFACE  ON  DECEMBER  9,  1981 


Time 

Bottom  (50  cm) 

Middle  (25  cm) 

Surface  (0  cm) 

0820 

11.6 

11.8 

11.9 

0925 

11.6 

11.9 

12.1 

1020 

11.8 

12.0 

14.0 

1120 

11.7 

12.0 

13.5 

1240 

11.7 

11.7 

14.0 

1350 

12.0 

12.0 

16.5 

1440 

12.0 

12.3 

17.0 

1550 

12.0 

13.8 

16.0 

Table  2 


CONCENTRATIONS 

OF  TNT  AND 

RDX  IN  WATER 

AT  LAGOON 

CENTER  AND  SOUTHEAST  EDGE  AS 

A  FUNCTION  OF 

TIME  OF 

DAY 

A 

TNT 

(ppm) 

RDX 

(ppm) 

Time  of  Day 

Edge 

Center 

Edge 

Center 

» 

«  , 

0800 

8.34 

4.81 

12.7 

12.9 

0900 

7.11 

4.35 

13.2 

12.5 

e) 

1000 

7.00 

4.95 

13.6 

12.9 

1100 

7.62 

5.02 

13.0 

13.4 

1200 

6.80 

4.88 

12.7 

12.6 

1300 

6.54 

4.78 

12.9 

13.2 

(4 

1400 

7.27 

4.99 

12.8 

12.9 

;; 

1500 

7.28 

4.71 

12,8 

13.2 

", 
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III.  CHEMICALS  IDENTIFIED  IN  LAGOON  NO.  9 


Water  samples  were  collected  from  the  lagoon  during  site  visits 
In  December  1981,  May  1982,  and  July  1982.  In  addition,  samples  were 
sent  to  SRI  by  LAAP  personnel  at  two-week  Intervals  from  June  1982 
until  the  lagoon  became  dry  In  late  August  1982. 

The  samples  were  analyzed  by  HPLC,  using  direct  aqueous  Injection 
of  the  sample  after  the  addition  of  an  Internal  standard  (3 , 5-dinitro- 
toluene) . 

The  following  conditions  were  employed: 

Instrument:  Spectra-Physics  Model  3500B  Liquid  Chromatograph 

Column:  Waters  Assoc  C,0  Radial-Pak  A 

Solvent:  30%  B  in  A  to  80%  B  in  A.  A  *  H20; 

B  ■  methanol/acetonitrile  (50/50) 

Flow  Rate:  2.0  ml  min”1 

Detection:  UV  @  254  nm 

A  typical  HPLC  profile  appears  in  Figure  3.  The  chemicals 
identified  in  this  profile  and  confirmed  by  probe-mass  spectrometry 
were  octahydro-l,3,5,7-tetranitro-l,3,5,7-tetrazocine  (HMX),  RDX, 
1,3,5-trlnitrobenzene  (TNB),  3,5-dinitroanlllne  (DNA),  TNT,  2-araino- 
4,6-dinitrotoluene  (2-A-4,6-DNT),  and  4-amino-2 , 6-dlnitrotoluene 
(4-A~2,6-DNT) .  Concentrations  of  each  chemical  in  water  samples  at 
various  dates  are  shown  in  Table  3. 

Investigations  were  also  conducted  to  identify  stable  munitions- 
related  chemicals  in  the  lagoon  water.  A  400-ml  sample  of  lagoon 
water  was  extracted  with  benzene  (2  x  400  ml)  (benzene  fraction)  and 
ether  (2  x  200  ml)  (ether  I  fraction);  then  the  water  was  acidified  to 
pH  1  with  HC1  and  re-extracted  with  ether  (2  x  200  ml)  (ether  II 
fraction).  Each  fraction  was  dried  over  anhydrous  sodium  sulfate  and 
rotary-evaporated  to  dryness.  The  residues  were  dissolved  in  a  small 
amount  of  ether  for  chromatographic  analysis. 
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Table  3 


CONCENTRATIONS  (ppm)  OF  MUNITIONS  FOUND 
IN  WATER  SAMPLES  ON  VARIOUS  DATES 

Date 


Sampled 

Day 

HMX 

RDX 

TNB 

4-A-2 , 6-DNT 

2-A-4 , 6-DNT 

TNT 

DNA 

12/08/81 

0 

— 

13.2 

0.33 

1.10 

1.50 

6.98 

1.0 

3/15/82 

74 

2.87 

10.9 

0.58 

— 

1.01 

5.49 

<0.10 

5/25/82 

145 

3.15 

18.3 

0.65 

0.91 

1.33 

19.8 

0.54 

6/09/82 

160 

4.31 

22.4 

0.91 

2.35 

1.33 

11.5 

0.54 

6/23/82 

174 

4.40 

20.9 

0.75 

2.31 

1.39 

5.65 

0.68 

7/07/82 

188 

5.71 

28.9 

0.84 

0.83 

1.48 

6.31 

0.99 

7/22/82 

203 

4.03 

17.1 

0.20 

0.42 

0.77 

2.81 

0.73 

8/04/82 

216 

6.36 

42.1 

4.17 

0.58 

1.37 

12.3 

0.87 

8/20/82 

232 

6.03 

12.6 

2.06 

0.65 

0.83 

0.  50 

0.34 

9/03/82 

246 

5.07 

5.6 

<0.10 

— 

1.16 

0.20 

<0.10 

We  identified  TNT,  2-amlno-4 ,6-dlnitrotoluene,  4-amino-2 ,6- 
dinitrotoluene,  RDX,  and  trace  amounts  of  1 , 3, 5-trlnitrobenzene  and 
2 ,4,6-trlnltrobenzyl  alcohol  by  thin-layer  chromatographic  (TLC) 
analysis  on  silica  gel,  using  benzene  as  a  solvent  and  an 
ethylenediamlne/DMSO  spray  reagent.  In  the  ether  I  fraction,  only  one 
spot  migrated,  which  coincided  with  RDX.  In  the  ether  II  fraction, 
several  spots  were  observed  close  to  the  origin  but  could  not  be 
ldentlfed  by  cochromatography  with  standards  or  color  development  with 
the  spray  reagent. 

The  ether  II  fraction  was  treated  with  dlazomethane  and  subjected 
to  HPLC  analysis  under  the  following  conditions: 

Column:  Water’s  Clg  Radial  Pak  A 

Solvent:  0-992  methanol  in  water  with  a  10  rain  delay  and  a 

SO  min  linear  gradient 

Flow  rate:  2.0  ml  raln”^- 

Detection:  UV  @  254  run 


A  typical  HPLC  chromatogram  appears  in  Figure  4.  Components  from 
the  HPLC  profile  in  Figure  4  were  collected  and  evaluated  by  probe 
mass  spectrometry.  The  fraction  57.34  component  was  identified  as  the 
methyl  ether  of  3, 5-dinitrophenol.  Fraction  54.31  wa6  tentatively 
Identified  as  methyl  2-hyd roxy-4 , 6-dinitroben2oic  acid,  methyl  ether 
and  methyl  ester,  but  could  not  be  confirmed  due  to  the  lack  of  a 
reference  standard.  Fraction  54.31  also  contained  a  component  that 
was  identified  as  methyl  2-amino-4 ,6-dinitrobenzoic  acid,  methyl 
ester,  and  was  confirmed  by  the  synthesis  of  a  reference  standard  in 
our  laboratory.  The  other  component  could  not  be  identified  due  to 
multicomponent  spectra  and  elevated  background  signals.  Fraction 
43.50  was  identified  as  methyl  2,4-dinitrobenzoic  acid,  methyl  ester. 

Thus,  the  majority  of  chemicals  'n  the  LAAP  lagoon  water  are 
related  to  the  photochemical  transformation  of  TNT  with  the  exception 
of  2-amino-4,6-dinltrotoluene  and  4-amino-2 , 6-dinitrotoluene,  which 
arise  from  microbial  sources.  Based  on  the  work  of  Kaplan  et  al. 
(1975),  who  showed  that  many  of  the  identified  components  are  the 
result  of  photochemical  processes,  Scheme  I  indicates  the  probable 
route  of  formation  of  the  identified  products. 
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Scheme  I 
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57.34  3.5-Oiniwophenol.  Methyl  Etl.ei 


A  key  step  In  this  scheme  is  the  series  of  oxidative-reductive 
disproportionations  leading  to  2-araino-A ,6-dinitroben2oic  acid,  as 
shown  in  Scheme  II. 


CHO  CHC 


NO,  NO, 


|  diiP'OPO'. 


NO,  NO, 


Scheme  II 


All  the  chemicals  identified  in  the  lagoon  water  were  also  found 
in  the  sediment.  In  addition,  the  explosive,  tetryl,  was  found  in  the 
sediment. 

Because  2-amino-A, 6-dinitrobenzoic  acid  appeared  to  be  a  major 
transformation  product  in  the  lagoon  water,  we  developed  an  HPLC 
procedure  to  monitor  this  compound.  The  following  conditions  were 
employed: 

Column:  250  mm  *  A. 6  mm  Alltech  Clg  10  pm 

Solvent:  AOZ  B  in  A  +  50-’’  B  in  A  in  10  rain  after  a 

5-min  hold.  A  ■  0.1M  phosphate  buffer  (pH  7) 
with  0.05  M  tributylamraonium  phosphate  (TBAP); 

B  ■  0.005  M  (TBAP)  in  methanol. 

Flow  rate:  1.8  ral  mln”^ 

Detector:  UV  @  254  nm 

Retention  time:  3.62  rain 
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A  typical  chromatographic  profile  of  lagoon  water  appears  in 
Figure  5.  Lagoon  waters  that  were  collected  throughout  this  study 
were  analyzed  for  2-amino-4 , 6-dinitrobenzolc  acid.  The  results  (Table 
4)  suggest  that  2-amino-4,6-dinitrobenzoic  acid  is  a  major  trans¬ 
formation  product  in  lagoon  waters. 

The  lagoon  waters  were  found  to  be  nearly  saturated  with 
oxygen.  Measurements  taken  with  a  dissolved  oxygen  probe  during  the 
December  field  trip  Indicated  9.0  ppm  0j  at  the  surface  and  8.8  ppm 
near  the  bottom  sediment  (50  cm  deep).  In  such  a  highly  aerated 
environment,  nitrite  is  readily  oxidized.  In  fact,  when  the  lagoon 
water  collected  in  December  1981  was  analyzed  for  nitrite  and  nitrate 
according  to  the  method  of  Thayer  (1979),  only  nitrate  was  found  (16.7 
ppm);  the  pH  of  the  lagoon  waters  ranged  from  7.3  to  7.8,  and  the 
total  dissolved  solids  concentration  was  540  ppm. 

Absorbance  measurements  were  obtained  at  290  and  400  nm  as  indi¬ 
cators  of  the  light-absorbing  properties  of  the  lagoon  water  as  a 
function  of  time  of  year.  These  values,  reported  in  Table  5  along 
with  the  measured  water  volume  at  the  sampling  date,  indicate  that  the 
UV-absorbing  materials  are  concentrated  as  the  lagoon  waters 
evaporate,  as  evidenced  by  a  3-fold  Increase  In  the  absorbance  between 
December  and  August  with  the  concurrent  decrease  in  volume  by  a  factor 
of  three  over  the  same  period.  This  has  the  effect  of  reducing  the 
penetrability  of  light  into  the  lagoon  water  in  the  summer  months 
compared  to  winter. 
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LA-7934-32 

FIGURE  5  HPLC  PROFILE  OF  LAGOON  WATER  RESOLVING 

2-AMINO-4,6-DINITRO  BENZOIC  ACID  AND  OTHER 
COMPONENTS 


Table  4 

LEVELS  OF  2-AMINO-4.6-DINITROBENZOIC  ACID  OBSERVED  IN 
LAGOON  NO.  9  ON  VARIOUS  SAMPLING  DATES 


Date  of 
Sample 


2-Amlno-4,6-Dlnitrobenzolc  Acid 
(mg  A-1) 


12/18/81 

28.9 

03/15/82 

24.3 

03/15/82 

24.3 

05/25/82 

30.7 

06/09/82 

35.2 

06/23/82 

43.9 

07/07/82 

45.5 

08/04/82 

49.9 

08/20/82 

Table  5 

45.7 

ABSORBANCE 

VALUES  OP  LAGOON  NO. 9 

WATER  AT 

290  AND  400 

ran  AT  VARIOUS  TIMES  OF 

Absorbance  (1-cm  cell) 

THE  YEAR 

Date 

Lagoon  Volume 

of  Sample 

290  nm  400  nm 

m3 

12/08/81 

03/15/82 

06/09/82 

06/23/82 

07/07/82 

08/04/82 

08/20/82 

09/03/82 


IV.  SITE  EVALUATION 
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On  the  basis  of  on-site  observations  and  analysis  of  the  lagoon 
waters,  transformation  products,  and  sediments,  we  conclude  that  the 
lagoon  is  a  complex  environment  governed  by  chemical  inputs  that  are 
dependent  on  lagoon  water  temperature  and  transformation  rates,  which 
in  turn  are  dependent  on  time  of  year,  natural  sensitizers,  and  depth 
at  which  measurements  are  made. 

Near  the  dump  ramp  we  found  levels  of  TNT  and  RDX  in  the  sediment 
(15  and  8%,  respectively)  that  were  significantly  greater  than  amounts 
that  would  result  from  sorption  alone.  This  suggested  that 
particulate  munitions  were  present;  in  fact,  when  the  lagoon  surface 
became  visible,  chunks  of  pure  TNT  appeared  on  the  surface  sediment. 
The  percent  of  munitions  in  the  sediment  decreased  as  a  function  of 
distance  from  the  dump  ramp.  Therefore,  chemical  input  to  the  lagoon 
is  controlled  by  diffusion  of  chemical  through  the  sediment  and 
dissolution  of  pure  chemicals  from  the  sediment  surface  to  the  bulk- 
water  phase.  Because  the  bottom  sediment  is  a  heterogeneous  system, 
we  could  only  estimate  chemical  inputs  based  on  limited  laboratory 
investigations. 

Transformation  products  found  in  the  lagoon  suggested  that  both 
biological  and  photochemical  transformations  were  occurring  in  the 
lagoon  (see  above).  Most  of  the  light  is  absorbed  in  the  first  few 
millimeters  of  the  lagoon;  consequently,  good  vertical  mixing  is 
required  to  renew  depleted  chemicals  as  rapidly  as  they  are 
transformed.  If  good  mixing  is  achieved,  then  a  photochemical  rate 
constant  can  be  estimated  as  a  function  of  depth.  Moreover,  because 
absorption  of  the  sunlight  by  the  lagoon  produces  radiant  heat,  lagoon 
water  temperature  varies  widely  throughout  the  year;  for  example,  it 
was  10°C  in  December,  41.5°C  in  July.  Thi9  increase  in  temperature 
could  have  a  large  impact  on  microbial  populations  and  hence  on 
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biotransforraation  rates;  in  addition,  it  could  accelerate  diffusion 
and  dissolution  rates  of  chemicals  from  the  bottom  sediment. 

Estimation  of  the  concentration  of  a  chemical  at  any  particular 
time  is  further  complicated  by  the  water  balance  within  the  lagoon. 
Water  enters  the  lagoon  through  precipitation  and  is  lost  through 
evaporation  and  seepage.  Therefore,  an  understanding  of  water  input 
and  output  rates  is  an  integral  part  of  this  fate  assessment. 

On  the  basis  of  the  above  evaluations,  we  perceived  a  number  of 
processes  that  were  affecting  the  concentration  of  chemicals  in  the 
lagoon  environment.  Processes  that  would  increase  concentration  of  a 
chemical  in  the  lagoon  water  include  1)  water  evaporation  and  seepage 
2)  dissolution  of  a  pure  chemical  from  the  sediment  surface,  and  3) 
diffusion  of  munitions  from  the  sediment.  Processes  that  would 
decrease  the  concentration  in  the  water  include  1)  precipitation,  2) 
photochemical  transformation,  and  3)  biological  transformation.  We 
therefore  performed  laboratory  investigations  of  these  processes  to 
provide  rate  and  equilibrium  constants  that  are  applicable  to  the 
lagoon  environment.  From  these  data,  we  constructed  a  computer  model 
to  simulate  the  processes  in  the  lagoon  and  compared  the  simulation 
data  with  actual  data  obtained  for  the  lagoon. 


V. 


TRANSPORT  AND  TRANSFORMATION  STUTII'S 


A .  Physical  Transport 

1 .  Water  Balance 

The  volume  of  water  in  the  lagoon  is  controlled  primarily  by  pre¬ 
cipitation,  evaporation,  and  vertical  seepage.  Surface  runoff  or 
horizontal  seepage  was  considered  to  have  little  impact  on  the  net 
volume  changes  due  to  the  flat  bottom  and  banking  of  the  lagoon. 

We  created  a  water-volume  balance  expression  using: 


dV 

dt 


(D 


where  R„  ■  r  A 
P  P 

Rev  ‘  revA 

R*  “  rAA 

and  rp  represents  precipitation,  rgv  represents  evaporation,  r^  repre¬ 
sents  leaching  (all  in  era  day-1),  and  A  represents  the  surface  area  of 
the  lagoon.  Average  values  of  rgy  and  rp  were  obtained  from  three 
years  of  weather  data  (1979-1981)  collected  by  the  Red  River  weather 
station  near  LAAP;  these  are  shown  in  Table  6.  Similar  data  for  1982 
through  August  are  presented  in  Table  7. 

Leaching  is  usually  estimated  by  a  difference  determination  of 
volumes  when  the  volume  loss  from  evaporation  is  known.  During  the 
December  field  study  we  calculated  the  leaching  rate  to  be  0.17  cm 
day-1.  However,  we  used  this  parameter  as  a  variable  in  the  model  to 
adjust  the  overall  water  balance  to  approximate  that  observed  in  the 


Table  6 

AVERAGE  MONTHLY  PRECIPITATION  AND  EVAPORATION  DATA  FROM  1979 
THROUGH  1981  AT  RED  RIVER  WEATHER  STATION,  RESTIN,  LA 


Month 


Precipitation 

(cm/day) 


Evaporation 

(cm/day) 


1 

0.28 

0.13 

2 

3 

4 

0.30 

0.46 

0.34 

0.20 

0.35 

0.45 

5 

6 

7 

0.61 

0.28 

0.39 

0.52 

0.67 

0.69 

8 

9 

0.12 

0.18 

0.66 

0.50 

10 

0.44 

0.38 

11 

0.25 

0.20 

12 

0.16 

0.17 

Table  7 

MONTHLY  PRECIPITATION 
THROUGH  AUGUST  1982  AT 

AND  EVAPORATION  DATA  FROM  JANUARY  1982 
RED  RIVER  WEATHER  STATION,  RESTIN,  LA 

Month 

Precipitation 
( cm/day) 

Evaporation 

(cm/day) 

Jan 

0.31 

0.13 

Feb 

0.26 

0.15 

Mar 

0.19 

0.29 

Apr 

0.27 

0.34 

May 

0.17 

0.59 

June 

0.22 

0.61 

July 

0.25 

0.62 

Aug 

0.17 

0.62 
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Using  the  data  in  Tables  6  and  7,  we  plotted  the  depth  of  the 
lagoon  as  a  function  of  time  of  year  for  various  seepage  rates  (Figure 
6).  Fi om  depth  measurements  taken  at  various  time  points  in  the 
lagoon  we  found  that  a  seepage  rate  of  0.12  cm  day-1  gave  the  best 
approximation  to  the  actual  depths.  This  seepage  rate  projected  the 
lagoon  to  go  dry  on  Day  234  (August  22);  the  lagoon  actually  went  dry 
during  the  following  week.  We  therefore  used  the  depth  profile 
defined  by  the  0.12  cm  day  *  seepage  rate  curve  to  determine  the  water 
volume  on  any  particular  day  and  to  define  the  depth  that  was  critical 
in  estimating  the  magnitude  of  the  photochemical  rate  constant  at  any 
particular  time. 

2.  TNT  and  RDX  Solubility 

The  aqueous  solubilities  of  TNT  and  RDX  were  determined 
separately  at  10,  20,  and  30°C  and  simultaneously  at  20°C  by  the 
method  of  May  et  al.  (1978).  Briefly,  glass  beads  were  coated  with 
TNT  or  RDX  and  packed  into  a  column.  Then  water  was  percolated 
through  the  column.  After  an  equilibration  period,  the  concentration 
of  TNT  or  RDX  was  measured  in  the  column  effluent.  As  a  check  on  this 
method,  the  solubilities  of  TNT  and  RDX  were  measured  at  20°C  by  the 
method  of  Campbell  (1930),  in  which  the  solids  are  equilibrated  with 
water  in  a  therraostated  bath  and  filtered.  The  solubilities  using  the 
two  methods  are  shown  in  Table  8. 

Table  8 


SOLUBILITY  OF  TNT  AND  RDX  IN  WATER  AT  10,  20,  AND  30 °C 

TNT  (mg/liter) _ _ RDX  (mg/liter) _ 

Method  iO’C  20°C  30 0C  1Q°C _  20cC  30°C 

May  et  al.  67  i  1  104  ±  2  165  ±  2  21.9  ±  0.3  38.4  ±  0.4  67.0  ±  0.8 
Campbell  113  i  8  39  ±  2 
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Solubilities  were  measured  by  the  May  method  first  at  10  ,  then 
20°,  then  30°C,  and  finally  again  at  10°C.  No  significant  differences 
were  noted  between  the  two  measurements  at  10°C.  The  similarity  of 
the  measurements  at  10°C  demonstrates  the  repeatability  of  the  method 
after  equilibrating  the  column  and  showed  that  the  column  was  not 
exhausted.  Concentration  versus  flow-rate  measurements  were  made  to 
check  that  the  saturation  capacity  of  the  column  had  not  been 
exceeded.  Finally,  a  column  containing  both  RDX  and  TNT  was  used  to 
measure  the  solubility  of  one  compound  in  the  presence  of  a  saturated 
solution  of  the  other  compound.  Solubilities  obtained  for  RDX  and  TNT 
independently  were  in  agreement  within  the  standard  error  of  the 
values  for  RDX  and  TNT  in  the  presence  of  the  other  compound. 

The  results  of  the  two  methods  are  in  agreement  within  the  error 
of  each  measurement.  The  Campbell  solubility  values  show  a  greater 
error  than  those  using  the  May  method,  perhaps  because  the  Campbell 
method  requires  more  handling  of  the  samples. 

Because  the  temperature  of  the  lagoon  water  ranged  from  10°  to 
4l°C,  we  developed  mathematical  expressions  to  use  in  the  computer 
modeling  to  correlate  solubility  with  temperature.  These  expressions 
are  shown  in  Equations  2  and  3: 

7  O 

STNT  -  5.  567  x  LO  exp  (-"■—)  pg  ml"1  ;  (2) 


7  o 

SRDX  *  U-9U  X  10  exP  C479t-  P8  ml"1 


The  neasured  and  calculated  solubility  values  for  TNT  and  RDX  at 
10°,  20°,  and  30°C  (283,  293,  and  303°K)  are  shown  in  Table  9. 
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Table  9 


MEASURED  AND  CALCULATED  SOLUBILITIES  (rag  liter*1) 
OF  TNT  AND  RDX  AT  283,  293,  AND  303 *K 


Water  Temperature 
(°K) 


_ TNT _ 

Measured  Calculated 


_ RDX _ 

Measured  Calculated 


3.  Sediment  Sorption 

It  Is  generally  accepted  that  adsorption  of  organics  on  soils  and 
sediments  usually  fits  the  “mpirlcal  Freundlich  isotherm  (Equation  4): 


C  -  KC 

8  W 


where  Cg  is  the  equilibrium  concentration  on  the  sediment  (pg  g  ) ,  Cw 
is  the  equilibrium  concentration  in  the  aqueous  phase  (pg  nl*1  -  pg 
g  1 ) ,  and  K  and  n  are  constants.  Thus,  a  plot  of  the  logarithms  of  Cg 
versus  Cw  has  a  slope  of  n  and  an  intercept  of  K.  If  n  is  equal  to  1, 
the  adsorption  isotherm  is  linear  and  K  is  equal  to  Ks.  If  pure  TNT 
and  RDX  are  present,  as  in  the  LAAP  lagoon,  Cy  i3  the  saturated 
concentration  regardless  of  Cg  and  the  Freundlich  isotherm  may  not 
hold. 

We  performed  several  isotherm  measurements  and  fitted  the  data  to 
both  the  linear  (Cfl  *  KpCw)  and  nonlinear  (log  Cg  “  log  Kp  +  nlog  Cw) 
cases.  Very  high  values  of  Kp  were  obtained  from  the  nonlinear 
equation  and  were  rejected  in  favor  of  more  reasonable  values  obtained 
from  the  linear  equation.  These  Kp  values,  along  with  the  calculated 
Koc  values  (based  on  a  measured  organic  carbon  content  of  1.5%), 
appear  in  Table  10. 
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Table  10 


SORPTION  PARTITION  COEFFICIENTS  FOR  TNT  AND 
RDX  IN  LAAP  LAGOON  NO.  9  SEDIMENT 


Che.nical  Kp  KQC 


TNT 


7  ±  2  470  ±  130 


RDX 


4  i  7  270  ±  470 


4 .  Distribution  of  TNT  and  RDX  in  Sediment  Cores 

Several  core  samples  were  taken  throughout  the  lagoon  during  the 
December  1982  visit  to  determine  the  horizontal  and  vertical  distri¬ 
bution  of  RDX  and  TNT  in  the  lagoon  sediment.  A  51-cm  core  sediment 
sample  from  the  southeast  corner  of  the  lagoon  was  cut  into  5-cm  seg¬ 
ments.  Each  segment  was  analyzed  for  percent  TNT,  RDX,  and  dry 
sediment.  The  results  are  presented  in  Table  11. 


Table  11 

SEDIMENT  CORE  SAMPLE  DRY  WEIGHT  AND  PERCENT 
TNT  AND  RDX  WEIGHTS 


Core  Segment 
(cm  from  top) 

%  Dry  Weight 

%  TNT3 

%  RDXa 

5 

71.5 

9.34 

12.6 

5-10 

76.8 

2.74 

1.28 

10  -  15 

79.4 

2.38 

0.83 

15  -  20 

80.3 

1.25 

0.39 

20  -  25 

83.2 

0.14 

0.004 

25  -  31 

82.8 

0.010 

0.005 

31  -  36 

82.3 

0.0052 

0.0065 

36  -  41 

82.3 

0.0052 

0.0065 

41  -  46 

83.2 

0.023 

0.029 

46  -  51 

74.2 

0.04 

0.04 

aBased  on  dry 

weight  of  sediment 

. 

Saturated  sediment  Is  represented  by  a  concentration  of  >0,60% 
for  TNT  and  >0.033%  for  RDX.  Although  the  TNT  and  RDX  concentrations 
in  the  sediment  decrease  with  depth,  the  data  suggest  that  the  amounts 
exceed  the  absorption  capacity  for  TNT  and  RDX  at  depths  less  than 
25  cm  and  greater  than  50  cm,  respectively.  The  high  amounts  of  TNT 
and  RDX  in  the  surface  sediments  suggest  that  the  sediment  is  a  source 
of  these  munitions  in  the  lagoon  water. 

To  determine  the  horizontal  distribution  of  TNT  and  RDX,  we  took 
another  set  of  core  samples — three  from  the  southeast  (SE)  corner  and 
one  each  from  the  northwest  (NVJ),  southwest  (SW),  and  northeast  (NE) 
corners  of  the  lagoon.  The  cores  were  about  12  cm  in  depth  and  were 
cut  into  1-cra  segments,  which  were  extracted  and  analyzed  by  HPLC  for 
their  RDX  and  TNT  content.  The  results  are  presented  in  Table  12. 

Only  the  top  and  bottom  three  sections  of  each  core  were 
analyzed.  Concentrations  of  TNT  and  RDX  were  highest  in  the  SE  cores, 
which  are  near  the  dump  ramp,  and  they  rapidly  decrease  with  depth,  as 
shown  by  comparing  the  top  three  sections  with  the  bottom  three.  The 
TNT  and  RDX  concentrations  in  the  NW,  NE,  and  SW  cores  are  low  and 
uniform  throughout  the  core  as  compared  with  the  SE  cores. 


5 .  MasB  Transfer  Rate  From  Sediment 

The  presence  of  sorbed  and  solid  munitions  chemicals  in  the  LAAP 
sediment  presents  a  complex  problem  in  establishing  the  mass  transfer 
rate  of  chemicals  to  the  bulk-water  phase.  The  processes  that  may  be 
occurring  in  a  sediment  pore,  described  pictorially  in  Figure  7, 
Include  diffusion,  seepage,  sorption-desorption  equilibria, 
transformation,  and  dissolution  of  solid  munitions.  These  proceses 
can  be  described  mathematically  and  their  sum  can  be  used  to  compute 
the  overall  mass  transfer  rate. 

Diffusion  is  generally  accepted  as  a  second-order  differential 
expression  relating  changes  in  the  chemical  concentration  in  the 
sediment  (Cgw)  with  respect  to  depth  (z)  times  a  diffusion  constant 
D,.  Seepage  is  a  first-order  differential  expression  relating 
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FIGURE  7  LAGOON  WATER  RESTS  ON  SEDIMENT 

Pond  and  sediment  are  considered  to  be  semi-infinite. 
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concentration  to  depth  tinies  a  seepage  velocity  tactor,  v.  The 
sorption-desorption  equation  was  defined  in  Equation  4.  The  only 
transformation  process  believed  to  be  occurring  in  the  sediment  is 

Table  12 

TNT  AND  RDX  IN  SECTIONED  CORES  FROM  LAAP  LAGOON  NO.  9 


Section  Identification 


Sample 

No.a 

SE-A 

SE-B 

SE-C 

NW 

sw 

NE 

TNT  (pg  g“ 

ii 

12 

0.05 

11 

501 

6.4 

367 

0.04 

0.55 

0.03 

10 

525 

71 

38 

0.05 

0.23 

0.1" 

9 

— 

109 

280 

0.02 

0.08 

— 

3 

91 

0.4 

5.5 

0.14 

0.16 

0.009 

2 

27 

1.2 

2.9 

0.09 

0.28 

0.03 

1 

26 

0.7 

1.4 

0.06 

0.13 

0.05 

RDX  (pg  g‘ 

ill 

12 

0.10 

11 

89 

3.0 

88 

0.01 

0.01 

0.02 

10 

48 

22 

226 

0.02 

0.02 

0.01 

9 

— 

51 

88 

0.01 

0.01 

— 

3 

4.7 

0.2 

2.1 

0.05 

0.05 

0.001 

2 

1.3 

0.3 

0.8 

0.05 

0.05 

0.02 

1 

1.3 

0.3 

0.9 

0.01 

0.01 

0.04 

aSample  12  is  at  surface,  sample  1  is  lowest  below  surface. 
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biotransformation.  Because  of  the  high  concentrations  of  munitions  in 
the  sediment,  biotransfortuatlon  is  believed  to  be  a  zero-order  process 
(kQ)  for  TNT  (see  Biotransformation,  Section  V.C). 

To  account  for  the  dissolution,  we  envision  solid  TNT/RDX  dis¬ 
persed  homogeneously  throughout  the  sediment.  The  dissolution  rate  of 

dw 

any  solid  TNT  or  RDX,  — ,  can  be  written  as  Equation  5: 


dw 

dt 


k  <C  . 
r  sat 


-  V 


(5) 


where  kf  is  an  experimentally  determined  dissolution  coefficient,  C8at 
is  the  saturation  concentration  of  the  dissolving  species  in  the 
liquid  phase,  and  is  the  actual  concentration  of  the  dissolving 
species  at  the  solid/liquid  interface.*  If  diffusion  from  the 
Interface  into  the  surrounding  liquid  is  fast  compared  to  dissolution, 
then  is  approximately  equal  to  the  concentration  of  the  dissolved 
species  in  the  surrounding  liquid.  In  this  limit, 


dw 

dt 


k  (C 
r  sat 


C  ) 

sw' 


(6) 


Therefore,  in  the  sediment,  the  change  in  chemical  concentration  with 

ac 

time  ( .  8W)  Is  equal  to  the  sum  of  the  dissolution  rate  (Equation  6), 

at  2 

a  c  ac 

8W  8  V 

diffusion  rate  (D ^  ^ — )>  and  seepage  rate  (-^y  )  minus  the  bio- 

dz* 

transformation  rate  (kQ,  zero-order  process  in  the  sediment)  and  the 

.  ac 

n-1  gy 

amount  of  chemical  that  is  sorbed  to  the  sediment  (K  nC  — - )  or 

p  sw  8t 


*Thi8  is  a  standard  "concentration  driving  force”  model  that  has  been 
adopted  by  many  workers  (Wurster  and  Kildsig,  1965). 
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m. 

*  V 

•*,V 
-■  S 


dt 


^dis^sat 


C  ]  +  D 
sw1 


1  dz2 


+  V 


.11  A  J  w 


dz 


k  -  K  nC 

O  P  sw  dt 


(7) 


Rearranging  Equation  7  yields  Equation  8: 


1  dC 

(1  +  K  nCn  )  -  k  ..  [C  -  C  ]  +  D 

p  sw  dt  dis  sat  sw 


d2C  5C 

sw  ,  sw  , 

1.2  dz  o 

dz 


(8) 


where  C  *  chemical  concentration  In  the  sedlraent/water  phase 
sw 

K  “  sediment  partition  coefficient  (C  ■  K  Cn  ) 
p  r  s  p  sw 

t  ■  time 

■  diffusion  coefficient  in  the  sediment /water  phase 

v  *  seepage  velocity 

z  *  distance  from  the  sediment 

k  ■  zero-order  rate  constant  for  blotransformatlon  In 
o 

sediment. 

Similarly,  for  the  free-water  phase,  the  change  In  concentration 
of  a  chemical  with  time  can  be  described  by 


where  Cw  =  concentration  of  chemical  in  the  free-water  phase 

t  ■  time 

D2  •  diffusion  coefficient  in  the  free-water  phase 

v  "  seepage  velocity 

Ek.C,  ■  first-order  rate  expression  for  the  transformation 
process. 


Solutions  to  these  equations  and  the  assumptions  that  are  used  to 
arrive  at  solutions  appear  in  Appendix  B.  However,  laboratory  studies 
were  necessary  to  determine  so  that  a  flux  from  the  lagoon  sedi¬ 

ment  could  be  computed. 


Determination  of  TNT  and  RDX  Flux  Values 
for  Sediment  from  LAAP  Lagoon  No.  ^ 


From  the  sediment  analyses  and  the  visual  observations  of  pure 
TNT  and  RDX  on  the  sediment  surface,  we  recognized  that  the  bottom 
sediment  was  the  predominant  Input  source  of  TNT  and  RDX  to  the  lagoon 
waters  during  the  period  of  our  study.  Therefore,  we  designed  experi¬ 
ments  to  determine  the  flux  of  these  chemicals  to  the  bulk-water 


phase. 


To  determine  the  fluxes,  460  g  of  LAAP  Lagoon  No.  9  sediment 
containing  12%  TNT  and  9%  RDX  was  placed  in  a  1000-ml  beaker  (surface 
area,  780  cm)  equipped  with  a  water-inlet  port  directly  above  the 
sediment  and  a  water-outlet  port  near  the  top  of  the  beaker.  The 
beaker  was  wrapped  in  aluminum  foil.  Water  was  pumped  into  the  beaker 
from  a  Milton-Roy  precision  pump  at  various  flow  rates  and  stirred 
slowly  with  an  overhead  stirrer.  The  effluent  was  collected  to  verify 
total  water  volumes  passing  through  the  beaker.  TNT  and  RDX  in  the 
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effluent  were  monitored  until  equilibrium  values  were  obtained 
(usually  In  3  to  4  days).  We  calculated  the  flux,  N,  from  Equation 
10: 


where  C  ■  concentration  of  chemical  at  equilibrium 
Q  -  flow  rate  through  the  beaker 
A  ■  surface  area  of  the  sediment. 


(10) 


m 


l-y. 


Values  for  the  fluxes  of  TNT  and  RDX  at  various  flow  rates  are 
shown  In  Table  13.  We  would  expect  that  as  the  flow  rate  Increased, 
the  concentration  of  chemical  would  decrease  and  the  flux  would 
Increase.  The  changes  In  concentration  between  runs  A  and  B  for  TNT 
and  B  and  C  for  RDX  do  not  follow  the  proper  concentration  trend. 
Also,  the  changes  In  flux  between  runs  D  and  E  for  TNT  and  RDX  do  not 
follow  the  proper  flux  trend.  These  observations  cannot  be  explained 
by  any  mass  transfer  model;  however,  we  cannot  be  certain  whether  any 
of  the  data  points  are  more  or  less  valid  than  others  and  must 
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allow  for  the  Imprecision  of  the  data  in  our  interpretations.  A  plot 
of  TNT  and  RDX  flux  versus  concentration  appears  in  Figure  8. 

From  this  plot  and  linear  least-squares  analysis  (Appendix  B)  we 
calculated  the  TNT  flux  to  be  1.08  x  10  **  era  sec  *  and  the  RDX  flux  to 
be  7.75  x  10  “5  cm  sec-*  in  the  lagoon.  The  first-order  dissolution 
rate  constants  were  calculated  to  be  1.18  x  10  ®  cm  sec  *  for  TNT  and 
1.16  x  io"5  ^  8ec"l  for  RDX  (see  Appendix  B) . 

B.  Photochemistry 

Photochemical  transformation  of  TNT  in  the  lagoon  environment  was 

evident  from  the  identification  of  photolytic  transformation  products 

and  the  deep  orange-red  color  of  the  water  that  is  reminiscent  of  the 

"pink  water"  at  TNT  production  facilities.  Because  lagoon  water 

absorbs  all  UV  light  near  the  surface,  the  photochemical  rate  constant 

has  a  depth  dependence  of  the  form  k'  ■  k  /D.  Moreover,  the  value  of 

P  P 

kp/D  also  is  subject  to  variations  owing  to  changes  in  light  intensity 
with  changes  in  cloud  cover  and  season. 

1.  TNT  Photochemical  Rate  Constant— -General  Remarks 

Usually,  photochemical  kinetic  studies  are  performed  using 
solutions  that  absorb  less  than  5%  of  the  incident  light.  The  amount 
of  light  absorbed  at  wavelength  X  by  a  solution  is  given  by  the  Beer- 
Lambert  Law  (Equations  11  and  12): 


-  Io(l  -  io’A) 

> 

(11) 

exc*  *  ax* 

» 

(12) 

where  IQ  and  1^  are  the  incident  and  absorbed  light  intensities, 
respectively,  A  is  the  absorbance,  i  is  the  path  length  of  light  in 
solution,  C  is  the  concentration  of  chemical,  is  the  molar  absor¬ 
ption  coefficient  of  a  chemical  at  wavelength  X,  and  is  the 
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CONCENTRATION  IN  BULK  WATER  (jig/ml) 

LA-7934-26 

FIGURE  8  TNT  (•)  ANO  RDX  (O)  FLUX  FOR  LAAP  LAGOON  NO.  9  SEDIMENT 
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absorbance  of  solution  without  added  chemical  at  wavelength  X.  in  a 
1-cm  cell.  The  terra  is  useful  for  characterizing  a  natural  water 
because  the  molar  concent ation  of  humic  acids  and  their  values  in 
water  are  difficult  to  determine.  Because  1^  is  a  function  of  path- 
length  A,  the  value  for  1^  can  be  calculated  for  any  pathlength  in 
solution. 

If  the  solution  absorbs  less  than  2%  of  the  incident  light  (A  < 
0.02),  Equation  11  takes  the  form  of  Equation  13: 

lX  "  2*3lnA  *  (13) 


and  the  rate  constant  for  photolysis,  kp,  is  given  by  Equation  14: 


k  =  2. 3$e  1  A 
p  \  o 


(14) 


where  is  the  photolysis  reaction  quantum  yield  (Mabey  «t  al., 


1982). 


However,  if  the  solution  absorbs  a  significant  amount  of  light 
due  to  some  combination  of  long  pathlength  and/or  large  eC  or 
values,  then  the  photolysis  rate  constant  will  show  a  dependence  on 
pathlength  due  to  attenuation  of  I0  with  increased  absorbance,  and  the 
resulting  equation  for  kp  is  more  complicated. 


The  solar  photolysis  rate  constant  kpE  for  chemicals  in  aquatic 
systems  in  optically  thin  solutions  (l.e.,  near  the  surface)  can  be 
calculated  from  Equation  15: 


(15) 


where  b  is  a  unit  conversion  factor  and  e^I^  values  for  each 
wavelength  Interval  are  summed  over  the  entire  solar  spectral 
region.  If  the  water  has  a  large  a ^A  value,  Equation  15  becomes  more 
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complex  because  the  values  decrease  as  wavelengths  become  longer, 
and  therefore  longer  wavelength  light  penetrates  natural  waters  to 
greater  depths  than  does  shorter  (more  energetic)  wavelength  light. 

Table  14  shows  the  (a^  +  values  at  selected  wavelengths  for 

the  LAAP  water  collected  in  December  1981  and  the  depths  at  which  99% 
of  incident  light  at  each  wavelength  Is  absorbed.  These  data  clearly 
show  that  light  of  wavelengths  below  500  ran  will  not  penetrate  below 
the  top  few  centimeters  of  LAAP  water.  The  absorbance  of  the  water  is 
almost  entirely  due  to  substances  other  than  TNT  (a^  >>  e^C); 
absorbance  of  6  ppm  TNT  In  a  1.0-cra  cell  is  0.02  AU  at  320  nm  (e^  a 
8'i0  M  1  cm  ^3,  which  is  negligible  compared  to  the  2.60  AU  measured 
for  the  water  at  320  nm. 


Table  14 

DEPTH  FOR  99%  ABSORPTION  OF  LIGHT 
OF  SPECIFIC  WAVELENGTHS  FOR  LAAP  WATER 


Wavelength, 

\,  (nm) 

(“A  +  £XC>a 
(cm-1) 

99%A 

300 

3.32 

0.60 

320 

2.60 

0.77 

340 

2.32 

0.86 

360 

2.10 

0.95 

380 

1.93 

1.04 

400 

1.74 

1.15 

420 

1.47 

1.36 

440 

1.23 

1.62 

460 

0.95 

2.11 

480 

0.74 

2.71 

505 

0.58 

3.45 

585 

0.28 

7.04 

655 

0.07 

27.00 

aAbsorbance  measured  in  1.0-cm  cell;  light 
absorp  ion  is  entirely  due  to  substances  other 
than  TNT  in  the  water  (see  text). 
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The  TNT  photolysis  rate  constant  for  use  in  modeling  the  fate  of 
TNT  in  LAAP  water  must  account  for  the  loss  of  TNT  in  the  entire  vol¬ 
ume  of  the  lagoon.  Although  the  mass  of  TNT  undergoing  photolytic 
transformation  per  unit  time  is  dependent  only  on  light  intensity  and 
spectral  distribution  entering  the  water,  the  rate  of  change  of  TNT 
concentration  is  dependent  on  volume.  For  photolysis  of  TNT  in  solu¬ 
tion  out  of  the  photic  zone,  the  rate  constant  must  then  be  an  inverse 
function  of  depth.  The  theory  and  validation  of  the  procedure  for 
estimating  the  photolysis  rate  constant  in  LAAP  water  where  the  photic 
zone  is  shallow  compared  to  the  total  depth  of  the  water  are  developed 
and  discussed  in  the  next  section. 


2.  Effect  of  Depth  on  Photolysis  Rate  Constant  of  TNT 


in  LAAP  Water 


To  estimate  a  rate  constant  for  TNT  in  LAAP  water,  we  have 
developed  an  expression  that  relates  the  TNT  photolysis  rate  constant 
to  lagoon  depth  (D).  This  expression  is  valid  only  if  vertical  mixing 
in  the  lagoon  is  substantially  faster  than  the  photolysis  rate  of  TNT 
near  the  surface  so  that  the  TNT  concentration  is  the  same  throughout 
the  lagoon.  The  expression  also  assumes  that  TNT  decomposes  by 
indirect  photolysis;  the  data  in  Table  14  clearly  show  that  substances 
other  than  TNT  account  for  the  major  absorption  of  light,  and  our  pre¬ 
vious  studies  have  demonstrated  that  indirect  photolysis  of  TNT  in 
natural  waters  is  more  rapid  than  direct  photolysis  (Spanggord  et  al., 
1981).  Light  energy  is  absorbed  by  substances  present  in  LAAP  water 
(humic  and  fulvic  acids,  TNT  reaction  products,  and  other  chemical 
contaminants  in  the  lagoon  water)  and  is  then  transferred  to  or 
induces  reaction  of  TNT.  TNT  reacts  even  though  it  does  not  absorb 
light  directly. 


The  proposed  mechanism  for  TNT  photolysis  is  described  by 
Equations  16  to  19  below: 


S  +  TNT  -»  S  +  TNT 


(18) 


*  H 

TNT  -*•  (f)  (products)  +  (1  -  f)  (TNT)  ,  (19) 

where  S  Is  the  sensitizer,  S*  and  TNT*  are  excited  states,  and  f  is 
the  fraction  of  TNT  converted  to  products.  The  rate  of  photo¬ 
absorption,  R  (Equation  16),  per  unit  volume  is  shown  in  Equation  20: 

R  -  1qA/V  -  Io/D  ,  (20) 

where  IQ  is  the  incident  light  intensity  (Einstein’s  per  unit  area  per 
second),  A  is  the  exposed  surface  area,  V  is  the  volume  of  the 
container,  and  D  is  the  average  depth. 

The  quantum  yield,  4>,  for  TNT  loss  is  defined  in  Equation  21: 


4>  -  $s*k3(TNT]f/(k3[TNT]  +  k£)  ,  (21) 

where  is  the  quantum  yield  for  production  of  S*.  In  dilute 
solutions  of  TNT,  where  k2  »  k3(TNT], 


♦  *  ♦«.*!<, [TNTJf/k, 


(22) 


The  loss  rate  of  TNT  can  be  defined  as  the  TNT  quantum  yield 
times  the  rate  of  light  absorption  of  the  water  body.  Then, 


d(TNT]/dt  -  <>R  -  Jo^kjimif/ItjD  ,  (23) 


d(TNT]/{TNT)  -  fl  *  kj/k^ 


On  rearrangement  and  Integration, 


ln([TNT]0/[TNT]t)  -  fl^/k^t/D  -  k  t  (25) 


kp  ■  'VWW0  ‘ 


IQ  and  D  depend  on  sky  conditions  and  lagoon  geometry,  respectively. 

D  can  be  calculated  for  water  in  a  container  with  vertical  sines  and  a 
horizontal  bottom. 

By  making  the  substitution  k*  »  I  f (k- /k- >6- .  ir.to  Equation  26, 

p  O  J  Z  5" 

we  obtain  Equation  27,  which  states  that  kp  will  vary  Inversely  with 
depth : 


k  -  k'/D 
pD  p 


For  containers  of  depth  and  D^,  respectively,  relative  photolysis 
rate  constants  can  be  calculated  from  Equation  28: 
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Using  Equation  28,  we  can  calculate  the  TNT  photolysis  rate  constant 
In  the  lagoon  from  the  sunlight  photolysis  rate  constant  measured  for 
TNT  in  LAAP  water  in  a  shallow  container.  This  is  the  only  practical 
way  to  measure  kp  because  TNT  in  LAAP  lagoon  water  photolyses  too 
slowly  in  Che  lagoon  or  in  containers  of  the  same  depth  for  convenient 
measurement . 

3.  Experimental  Verification  of  Equation  27 


To  demonstrate  the  validity  of  Equation  27,  we  first  photolyzed 
LAAP  water  with  a  constant  intensity  source  in  a  merry-go-round 
apparatus  (MOR)  to  determine  the  TNT  concentretion-time  dependence  of 
Equation  25.  If  Equation  25  holds,  then  a  first-ordei  plot  (ln[TNT] 
versus  time)  will  be  linear.  For  these  experiments,  LAAP  water  con¬ 
taining  5  ppm  TNT  was  Irradiated  in  13  *  100-mm  borosllicatc  tubes 
with  monochromatic  light  at  313  or  366  nm  at  room  temperature. 

The  data  for  these  experiments  are  summarized  in  Table  15  along 
with  data  for  the  same  solution  in  tubes  exposed  to  sunlight  on 
January  22,  1982. 

Table  15 

PHOTOLYSIS  OF  4.97  ppm  TNT  IN  LAAP 
WATER  IN  AN  MGR  APPARATUS  AND  IN  SUNLIGHT 


Duration  of 


Final  TNT 


Light  Source 

Experiment 

(min) 

Concentration 

(ppm) 

103  k, 

j  (min" 

313  nm 

110 

3.29 

4.30 

O 

• 

o 

-H 

366  nm 

240 

2.34 

3.17 

±  0.11 

Sunlight 

80 

2.49 

8.51 

i  0.23 

(Jan.  22) 

First-order  treatment  of  these  data  gave  linear  plots,  as 
expected.  Similar  first-order  behavior  was  observed  in  all  other 
experiments . 
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We  then  verified  Equations  27  and  28  by  irradiating  LAAP  water  in 
cylinders  of  different  depths  using  a  solar  simulator  light  source — an 
intensity-regulated,  300-watt  xenon  lamp.  Radiation  from  the  lamp  was 
filtered  through  water  to  remove  infrared  wavelengths  and  then  through 
borosilicate  glass  to  remove  high-energy  UV  light.  Cylinders  were 
filled  to  depths  of  2.0,  4.2,  and  7.0  cm  and  Irradiated  from  above. 

The  sides  and  bottoms  of  the  cylinders  were  blackened  to  reduce 
reflection  and  to  prevent  light  from  entering  through  the  sides.  Data 
from  these  experiments  are  shown  in  Table  16. 

The  observed  photolysis  rate  constants  are  in  the  ratio 
4.05:1.75:1.0;  the  expected  values  are  3. 5:2. 1:1.0.  Stirring  had  no 
'ffect  on  the  photolysis  rate  constant:  when  the  stirrer  was  turned 
off  after  240  min  of  irradiation  of  the  dish  containing  4.2  cm  LAAP 
water,  the  rate  constant  did  not  change.  We  conclude  that  natural 
mixing  of  the  LAAP  water  in  our  experiments  was  sufficiently  fast  to 
keep  the  mixture  homogeneous  and  not  allow  depletion  of  TNT  in  the 
upper  irradiated  layers,  which  would  slow  down  the  reaction. 


,\ 
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Table  16 


PHOTOLYIS  OF  4.97  ppm  TNT  IN  LAAP  WATER 
AS  A  FUNCTION  OF  DEPTH  USING  A  XENON  LAMP® 


Duration  of 
Experiment 
(min) 

Water 

Depth 

(cm) 

Final  TNT 
Concentration 
(ppm) 

103  k„ 

(min~*)b 

236 

2.0 

1.81 

4.21 

i  0.C9 

542 

4.2 

r4 

00 

1.82 

*  0.02 

1075 

7.0 

1.33 

1.04 

±  0.03 

®300-W  controlled-current  lamp. 
br2  -  0.99. 
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Finally,  we  exposed  LAAP  water  to  sunlight  to  measure  TNT  rate 
constants  that  could  be  adjusted  for  depth.  Act*;-ometer  solutions  of 
PNA/PYR  were  exposed  tc  sunlight  in  parallel  studlep.  so  >,hat  cor¬ 
rections  for  sky  conditions  could  be  made  (discussed  in  the  next 
section).  These  data  permit  calculation  of  the  TNT  photolysis  rate 
constant  in  the  LAAP  lagoon  for  any  water  depth. 

To  measure  the  depth-independenc  photolysis  rate  constant  k^  for 
TNT  exposed  to  sunlight,  LAAP  water  filtered  through  a  0.45-micron 
filter  was  placed  in  opaque  glass  cylinders  whose  tops  were  exposed  to 
sunlight  on  March  12  and  15,  1982.  On  the  first  day  the  container  was 
filled  to  a  4.2  cm  depth  and  on  the  other  day  the  depth  was  3.5  cm. 

The  experimental  conditions  and  results  are  given  in  Table  17, 
uncorrected  for  sky  conditions. 

If  we  assume  that  an  Insignificant  amount  of  TNT  photolysis 
occurred  before  the  experiment  began  and  after  the  experiment  was 
ended  on  each  day,  the  measured  rate  constant  kp  (in  m_1)  can  be 
readily  converted  into  a  24-hour  day-averaged  rate  constant  kp  (in 
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d”1)  by  multiplying  the  observed  rate  constant  by  the  fraction  of  a 
24-hour  day  that  the  irradiation  period  represented.  Multiplication 
of  the  latter  by  the  depth  value  then  gives  a  value  of  3.0  cm  d”1, 
which  is  the  depth-independent  rate  constant,  k/,  used  in  Equation  26 
for  photolysis  of  TNT  in  LAAP  water  under  weather  conditions  of  mid- 
March.  To  obtain  a  rate  constant  for  a  specific  depth,  kV  is  divided 
by  the  lagoon  depth.  He  must  also  correct  for  sky  conditions 
during  the  experiments. 


Table  17 


IN  LAAP  WATER® »b 

AND  ACTINOMETERS 

Date 

Time 

3-12-82 

8:45-16:00 

3-15-82 

9:40-15:30 

Sky  Conditions 

bright,  sunny  w/ 
low  clouds 

cloudy  +  hazy, 
windy 

Depth,  cm 

4.2 

3.5 

Time,  min 

435 

350 

kp  x  103,  m"1 

1.67  *  0.045 

2.51  ±  0.041 

tl/2h 

6.9 

4.6 

kp  d-^  (day  averaged)0 

0.72 

0.89 

k'  cm  d~l 

P 

3.02 

3.12 

k '  cm  d” 1 , 

A  W 

4.14 

4.27 

clear  sky  conditions 


aPhotolyses  monitored  using  PNA/PYR  actlnometer  [PYR]  ■ 

1.25  *  10"3  M;  kp  (PNA)  -  4.8  x  10"3  m"1. 

bAll  rate  constants  are  uncorrected  for  sky  conditions. 

ckp  (day  averaged)  was  calculated  by  multiplying  the  observed 

photolysis  rate  constant  expressed  in  units  of  d~*  by  the 
fraction  of  a  24-hour  day  that  the  irradiation  represented. 
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Correction  of  Measured  Rate  Constants  for  Weather  Conditions 


Sunlight  photolysis  rate  constants  are  dependent  on  weather  con¬ 
ditions  as  well  as  the  diurnal  and  seasonal  light  Intensity  and  spec¬ 
tral  variation.  Diurnal  and  seasonal  variations  in  light  intensity 
(in  terras  of  photon  flux)  and  distribution  are  Included  In  procedures 
to  calculate  direct  photolysis  rate  constants  of  chemicals  in  opti¬ 
cally  thin  solutions  (Mill  et  al.,  1982).  These  rate  constants  are 
estimated  for  clear-sky  conditions  and  can  be  calculated  as  a  24-hour 
day-averaged  rate  constant.  Corrections  for  weather  conditions  (over¬ 
cast  or  cloudy  cover)  must  be  made  using  field  data;  they  were  per¬ 
formed  as  follows. 

During  the  LAAP  field  study  and  in  our  outdoor  studies  at  SRI,  we 
exposed  PNA/PYR  actlnoraeter  solutions  to  sunlight  simultaneously  with 
our  LAAP  water/TNT  experiments.  The  photolysis  rate  constants  for 
this  actlnoraeter  are  available  for  clear-sky  conditions  as  a  function 
of  season  (Dulin  and  Mill,  1982).  Rate  constants  calculated  for 
PNA/PYR,  with  (PYRJ  *  1.25  •'  10“3  M  ,  were  4.3  x  10“3  m_1  for  January 
and  9.6  x  10"3  nf 1  for  late  April  [winter  and  spring  season  dates, 
respectively,  using  the  solar  program  of  Zepp  and  Cline  (1977)].  By 
interpolation,  we  estimate  that  an  average  rate  constant  of  about  6.6 
x  10  3  m  1  would  then  be  expected  for  photolysis  of  PNA/PYR  in 
mid-March.  Because  the  actual  rate  constant  measured  for  photolysis 
of  PNA/PYR  on  both  days  was  4.8  *  10-3  m_1,  or  0.73  of  that  expected 
for  clear-sky  condition,  we  conclude  that  73Z  of  the  clear-sky  irrad- 
lance  was  available  during  our  experiments  with  LAAP  water. 

The  measured  rate  constant,  k' ,  should  therefore  be  divided  by 

P 

0.73  to  obtain  the  clear-sky  rate  constant.  The  depth-dependent  rate 
constant  for  photolysis  of  TNT  in  LAAP  water  in  mid-March  under  clear- 
sky  conditions  Is  3.1/0.73  “4.3  cm  d  For  LAAP  water  of  50-cm 
depth,  the  rate  constant  of  8.6  *  10“2  d_1  corresponds  to  a  half-life 
of  8  days. 
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The  above  calculations  assume  that  the  0.73  weather  correction 
factor,  fw,  obtained  from  the  actlnoraeter  experiment,  applies  to  light 
that  Induces  photolysis  of  TNT  In  LAAP  water.  Although  we  know  that 
this  assumption  Is  not  strictly  correct,  It  Is  sufficiently  accurate 
for  estimating  environmental  photolysis  rate  constants  of  chemicals 
that  absorb  In  the  same  region  of  the  spectrum.  The  assumption  may 
not  hold  in  the  case  of  LAAP  water  where  light  absorption  by  products 
occurs  Into  the  visible  region. 


5 .  Effect  of  Time  of  Year  on  Photochemical  Rate  Constant 
of  TNT  In  LAAP  Water 

Because  concentrations  of  TNT  were  found  to  persist  in  the  LAAP 
water  at  all  times  of  year.  In  part  due  to  dissolution  of  TNT  out  of 
the  sediments  (see  Section  V.A.5),  modeling  of  the  fate  of  TNT  In  the 
pond  water  required  information  on  how  the  photolysis  rate  constant  of 
TNT  in  the  water  will  vary  with  the  time  of  year  (l.e.,  seasonal  rate 
constants  are  needed).  Equation  25  shows  that  the  photolysis  rate 
constant  of  TNT  Is  seen  to  be  a  function  of  the  Incident  light  inten¬ 
sity  and  the  reaction  quantum  yield,  Although  we  have  no  Infor¬ 
mation  on  the  rate  of  photolysis  of  TNT  In  LAAP  water  as  a  function  of 
season,  the  relative  rate  constants  for  different  times  of  the  year 
can  be  calculated  if  we  assume  that  L^  values  of  Mabey  et  al.  (1982) 
are  useful  approximations  of  1^.  Then 


V^V25 
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(29) 


where  1  and  2  denote  different  times  of  the  year.  This  expression 
assumes  that  the  quantum  yield  for  photolysis  of  TNT  In  LAAP  water  is 
the  same  at  all  seasons.  Use  of  L^  values  to  estimate  the  1^  values 
for  this  calculation  is  not  strictly  correct  because  L^  values  are 
day-averaged  photon  flux  values  for  light  passing  through  optically 
thin  water,  and  therefore  they  are  corrected  for  path  length  and  the 
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refractive  index  of  water.  The  actual  flux  values  that  should  be  used 
in  the  calculation  are  the  flux  values  incident  on  the  surface  of  the 
water  with  no  path  length  correction  because  all  light  is  absorbed  by 
the  water.  This  flux  is  referred  to  as  by  Zepp  and  Cline  (1977), 
However,  since  the  day-averaged  analogs  of  are  not  known,  we  have 
used  the  values  as  approximations  of  the  1^  values  needed  for  the 
calculation;  comparison  of  and  values  listed  by  Zepp  and  Cline 
for  the  season  dates  at  40°  latitude  indicates  that  the  flux  values 
corrected  for  path  length  are  larger  than  the  uncorrected  values  by 
10-20Z  over  a  wide  wavelength  range. 

The  calculation  of  the  relative  photolysis  rates  of  TNT  in  LAAP 
water  also  requires  the  knowledge  or  an  assumption  that  the  absorption 
of  light  above  a  certain  wavelength  will  not  photolyze  TNT  due  to 
energy  limitations.  Therefore,  we  used  several  wavelength-cutoff  reg¬ 
ions  to  calculate  the  ratio  in  Equation  28  for  the  summer  and  winter 
solstices  and  the  ratio  for  the  September  and  March  dates  during  which 
outdoor  experiments  were  performed  (Table  18) .  The  data  in  the  table 
show  that  photolysis  of  TNT  in  LAAP  water  will  vary  by  a  factor  of 
about  3  between  winter  and  summer  seasons  if  only  sunlight  variations 
are  considered;  as  discussed  in  other  sections,  the  depth  of  the 
water,  weather  conditions,  and  possible  changes  in  the  nature  of  the 
lagoon  water  Itself  (absorbance  and  sensitizing  abilities)  will  also 
affect  the  photolysis  rates  of  TNT  in  the  LAAP  water  throughout  the 


Table  18 


VARIATION  IN  TNT  PHOTOLYSIS  RATES 
IN  LAAP  WATER  AS  A  FUNCTION  OF  TIME  OF  YEAR8 


Wavelength* 


S“\L\1/E(I\LX2 


Summer  season 


Cutoff 

Winter  season 

March 

400  nm 

3.3 

1.8 

500  nm 

2.9 

1.8 

600  nm 

2.8 

1,8 

Calculated  using  values  for  LAAP  water 
collected  in  July. 

and  values  above  this  wavelength 
not  used  in  calculation  of  rate  constant  ratio. 

CL^  value  used  for  summer  and  winter  solstices 
at  40°  latitude. 

values  used  were  interpolated  between 
season  dates. 


ft.  UV-Vlsihle  Spectra  of  LAAP  Water 

As  discussed  above,  the  photic  zone  of  the  December  LAAP  water 
was  thin  (1-2  cm)  compared  to  the  total  depth  of  the  water  (~60  cm) 
measured  during  the  December  site  visit.  In  subsequent  trips  in  May 
and  July  we  found  that  the  water  depth  was  decreasing;  therefore  the 
absorption  spectrum  of  the  water  was  measured  before  use  in  our 
experiments.  The  UV  spectrum  of  LAAP  water  in  December  is  shown  in 
Figure  9.  The  spectrum  of  the  May  LAAP  water  was  nearly  the  same  as 
that  of  the  December  water,  whereas  the  absorption  coefficients  of 
July  water  were  1.7  to  1.5  times  greater  than  those  of  the  December 
water  (Table  19).  From  the  Information  and  analyses  obtained  in  our 
field  and  laboratory  studies,  it  is  impossible  to  determine  how  the 
constituents  of  the  water  affecting  the  UV  spectra  changed  as  a  fun¬ 
ction  of  the  time  of  year,  although  it  is  clear  that  the  factor  of 
1.5-1. 7  shown  in  Table  19  is  less  than  a  3-fold  Increase  expected  from 
a  simple  concentration  effect,  considering  that  the  water  depth  of  ~60 
cm  in  December  decreased  to  approximately  18  cm  in  July. 


Table  19 

ABSORPTION  COEFFICIENT®  OF  DECEMBER  AND  JULY 
LAAP  WATERS  AT  SELECTED  WAVELENGTHS 


a^,  cm 


ax  (July) 


X  ran 


Dec  LAAP 


July  LAAP 


(Dec) 
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7.  Experimental  Studies  of  the  Photolysis  of  TNT  in  LAAP  Water 

The  photolysis  of  TNT  was  studied  In  different  LAAP  water  samples 
collected  during  the  field  trips  in  December  1981,  May  1982,  and  July 
1982;  no  additional  TNT  was  added  to  these  samples  for  our  experi¬ 
ments.  The  finding  that  the  character  of  the  lagoon  Itself  and  the 
LAAP  water  were  different  on  each  visit  (see  Section  IV)  required  that 
photolysis  of  TNT  be  studied  in  each  water  so  as  to  obtain  the  moat 
relevant  rate  constant  data  that  could  be  used  In  the  lagoon-modeling 
efforts.  The  results  from  the  several  studies  are  summarized  in  Table 
20.  The  studies  are  listed  In  the  chronological  order  in  which  they 
were  performed. 

The  rationale  for  the  sequence  is  as  follows.  The  studies  In 
March  (Expts.  2-4)  using  the  December  LAAP  water  were  conducted  after 
the  effect  of  depth  on  the  photolysis  rate  constant  was  found  to  be 
Important  (see  Section  V.B.3).  In  the  late  summer,  Experiments  6  and 
7  were  conducted  simultaneously  using  the  December  and  July  waters 
after  It  was  found  that  the  TNT  rate  constant  measured  in  the  July 
LAAP  water  In  Septerber  (Expt.  5)  was  nearly  the  same  as  that  measured 
In  March  (Expt.  3),  which  was  unexpected  because  photolysis  rates  are 
more  rapid  In  summer  than  In  spring  (see  Section  V.B.5).  When  these 
experiments  showed  that  the  photolysis  rates  of  TNT  were  different  in 
the  two  waters  at  the  same  time  of  year.  Experiments  8-10  were 
conducted  in  which  the  absorbances  of  December,  May,  and  July  waters 
were  made  the  same  by  dilution;  the  intent  of  these  experiments  was  to 
determine  whether  the  absorbances  of  waters  were  somehow  responsible 
for  the  different  rate  constants  observed,  or  whether  the  ability  of 
the  individual  waters  to  promote  TNT  photolysis  were  indeed  different 
among  the  water  samples  after  allowing  for  the  different 
absorbances.  Unfortunately,  the  error  in  the  measured  rate  constants 
1 8  too  large  to  clearly  distinguish  between  the  two  effects. 

From  the  data  presented  in  Table  20,  several  conclusions  can  be 
offered  that  are  important  in  estimating  the  photolysis  rate  constants 
of  TNT  in  LAAP  water.  On  the  basis  of  the  results  of  Experiments 
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8-10,  the  ability  of  the  LAAP  water  to  promote  photolysis  of  TNT 
appears  to  be  the  same  per  unit  absorbance  of  the  water,  although 
because  of  the  error  In  the  rate  constant  data,  we  cannot  entirely 
rule  out  the  possibility  that  the  difference  in  the  rate  constants  in 
Experiments  6  and  7  is  not  due  to  a  slight  effect  of  the  concentration 
of  natural  substances,  as  would  be  observed  by  their  absorbances. 


i 
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SUNLIGHT  PHOTOLYSIS  OP  TNT  IN  LAAP  WATER 
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The  data  in  Table  20  also  show  that  the  photolysis  rate  constant  for 
TNT  in  the  LAAP  waters  is  independent  of  TNT  concentration,  as  would 
be  expected  if  the  photolysis  of  TNT  were  a  first-order  process 
(Equation  26).  However,  ve  have  insufficient  daca  with  which  to 
evaluate  the  effect  of  changes  in  the  properties  of  the  lagoon  water 


on  the  photolysis  rate  constant. 

8.  RDX  Photolysis  Rate  Constant 

Concurrent  with  the  TNT  studies,  we  perfumed  calculations  to 
estimate  the  RDX  photolysis  rate  constant.  The  rate  of  RDX  photolysis 
can  be  described  by  Equation  30  for  wavelength  X: 


dC 

dt 


nji  -  io-(\  *  Cxc)t) 

D(ax  +  exC> 


(30) 


where  D  *  depth  (cm) 

IQX  •  incident  light  intensity 

aX  »  absorpt'fVlty  of  the  water 

eX  -  molar  absorptivity  coefficient  for  RDX  at  X 
i  ■  solution  path  length 


In  LAAP  water,  »  e^C.  Therefore,  Equation  30  above  reduces 
to  Equation  31: 


dc  «<,x(1  -  10'VKC 


dt 


Da, 


(31) 


*"(X  X. 

Since  the  exponential  function  10  X  is  also  very  small  compared  to  1 
in  waters  where  is  large.  Equation  31  reduces  to  Equation  32: 


Integration  of  Equation  32  gives  a  first-order  rate  expression  where 
kp  is  equal  to  the  expression  shown  below  at  wavelength  X.: 


♦1oXc\ 

Da, 


(33) 


The  environmental  rate  constant  can  now  be  summed  over  all  wavelengths 
in  the  solar  spectral  region  and  the  photochemical  rate  constant  can 
be  described  as  shown  in  Equation  34: 


♦^VoX 

Da, 


(34) 


Using  the  molar  absorptivity  values  measured  for  RDX  from  299  nm 
(e  ■  164)  to  400  nm  (e  »  0.02),  the  absorptivity  of  the  LAAP  water 
(a^)  over  the  same  spectral  range,  and  the  quantum  yield  ($  ■  0.16), 
we  calculated  the  photochemical  rate  constant  (kpE)  at  various  depths 
of  the  LAAP  water  using  a  modified  version  (GC  SOLAR)  of  the  computer 
program  developed  by  Zepp  and  Cline  (1977)  for  calculating  the  inten¬ 
sity  of  ultra-violet  radiation  (297.5-380  nm)  reaching  the  earth's 
surface  as  a  function  of  time  of  day,  longitude,  latitude,  and  sea¬ 
son.  These  data  are  shown  in  Table  21  for  the  four  seasons.  Also 
shown  are  the  depth  independent  rate  constants  for  RDX  computed  for 
LAAP  waters  collected  in  December  and  July.  From  the  rate  constants 
shown  in  parentheses,  a  comparison  of  the  seasonal  effects  on  photo¬ 
lysis  can  be  made;  however,  these  values  are  not  environmentally 
revelant  due  to  changes  in  the  water's  composition  as  a  function  of 
time  of  year. 
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Table  21 


PHOTOCHEMICAL  RATE  CONSTANTS  FOR  RDX 
IN  LAAP  WATERS  AS  A  FUNCTION  OF  SEASON  AND  DEPTH 


_ Depth  (cm) _ 

Surface8 

15 

30 

60 

90 

Depth  Independent  Rate 
Constant  (on  d  ')  in  July  LAAP 
Water'’ 

Depth  Independent  Rate 
Constant  (cm  d_1)  in  December 
LAAP  waterb 


Rate  Constant  (*103  d  ')  in 


July  LAAP  Water 


Fall 

Winter 

Spring 

Summer 

390 

200 

710 

930 

2.0 

1.0 

3.8 

5.0 

1.0 

0.52 

1.9 

2.5 

0.51 

0.26 

0.95 

1.3 

0.34 

0.17 

0.63 

0.83 

(30) 

(16) 

(56) 

76 

(39) 

19 

(75) 

(101) 

aRate  constant  at  surface  not  used 

bRate  constants  in  (  )  are  for  comparison  of  seasonal 

effects  on  photolysis  and  are  not  environmentally  relevant. 
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The  above  calculations  showed  that  the  light  intensity  for  a  par¬ 
ticular  LAAP  water  changed  by  a  factor  of  five  going  from  winter  to 
summer  months.  The  depth-independent  rate  constant,  however,  changed 
by  only  a  factor  of  three  due  to  the  greater  absorptivity  of  the  LAAP 
water  In  the  summer  (cc^  =  9.0  at  300  nn-)  than  in  the  winter  (a^  m  3.0 
at  300  nra) . 

To  determine  how  well  the  calculated  depth-independent  rate  con¬ 
stant  compared  to  a  measured  value,  we  performed  an  RDX  sunlight 
photolysis  experiment  in  a  19  *  3-cm  wide  crystallizing  dish  (darkened 
on  sides  and  bottom)  over  a  seven-day  period  (starting  9/18/82)  using 
lagoon  water  collected  on  24  July  1982  and  filtered  through  a  0.25-um 
filter.  A  plot  of  In  RDX  versus  time  gave  a  straight  line,  from  which 
a  first-order  rate  constant  (slope)  was  calculated  tc  be  4.77  x  10~5 
rain-*.  In  the  same  experiment,  TNT  photolyzed  with  a  pseudo-first- 
order  rate  constant  of  2.4  ±  0.3  x  10~3  min  * ,  or  approximately  50 
times  faster  than  RDX.  The  depth  independent  rate  constant  was  cal¬ 
culated  to  be  0.20  cm  day-*. 

This  value  is  about  three  times  larger  than  the  calculated  value 
for  July  LAAP  water.  Using  the  GC  SOLAR  program,  we  calculated  the 
RDX  photolysis  rate  constant  in  pure  water  to  be  2.0  cm  4ay~*  during 
the  summer,  which  is  in  good  agreement  with  an  experimental  value  of 
2.2  cm  day-1  measured  in  the  summer.  Therefore,  the  differences  in 
the  measured  and  calculated  photolysis  rate  constants  for  LAAP  water 
could  be  due  to  sensitization  reactlous  that  are  not  understood  at 
this  time.  For  the  modeling  study,  we  used  the  measured  rate  constant 
in  LAAP  water  corrected  by  the  light  intensity  factors  for  the  various 
seasons . 

C.  Blotransf ormat ion 

1 .  Bacterial  Count 

The  identification  of  aminodlnitrotoluenes  in  the  lagoon  water 
suggested  that  biotransfonna t ion  was  occurring  in  the  lagoon  water  and 
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sediment.  To  obtain  an  Indication  of  aerobic  and  aneroblc  microbial 
populations  In  the  lagoon  and  sediment,  bacterial  plate  counts  were 
performed  on  waters  and  sediment  collected  during  the  December  field 
study.  In  the  lagoon  water  (Table  22),  the  aerobic  bacterial  countB 
averaged  1.7  x  106  CFU  ml  *  and  the  anaerobic  bacteria  averaged  1.7  x 
103  CFU  ml  In  the  sediment  (Table  23),  the  averages  were  1.1  x  10 
CFU  ml-1  for  the  aerobic  bacteria  and  1.5  x  lO5  CFU  ml”1  for  the  anae 
robic  bacteria. 


Table  22 

TOTAL  AEROBIC  AND  ANEROBIC  BACTERIA 
IN  LAGOON  9  WATER  COLLECTED  IN  DECEMBER 


Sample 

Site 

Aerobic  Bacteria 
(CFU  ml"1) 

Anaerobic  Bacteria 
(CFU  ml"1) 

1 

1.6  x  106 

2.9  x  io3 

2 

• 

00 

X 

t— 

o 

<J\ 

1.1  x  103 

3 

1.7  X  106 

0.48  x  103 

4 

1.7  x  I06 

2.2  x  103 

Table  23 

IN 

TOTAL  AEROBIC  AND  ANAEROBIC 
LAGOON  9  SEDIMENT  COLLECTED 

BACTERIA 

IN  DECEMBER 

Sample 
SI  te 

Aerobic  Bacteria 
(CFU  ml"1) 

Anaerobic  Bacteria 
(CFU  ml"1) 

1 

0.51  x  LO7 

1.4  x  105 

/. 

1.4  x  107 

1.9  x  10s 

3 

0.76  x  IO7 

0.71  x  105 

4 

1.9  x  IO7 

2.0  x  IO5 

The  tea-colored  water  and  cedlraent  contained  unexpectedly  high 
numbers  of  aerobic  bacteria.  The  numbers  of  anaerobic  bacteria  were 
substantially  smaller,  reflecting  the  high  dissolved-oxygen  content  (7 


to  8  ppm)  of  the  water  in  the  winter.  Total  aerobic  bacteria  in  a 
water  sample  collected  on  25  May  1982  averaged  1.9  x  106  CFU  ral~*. 

A  lagoon  water  sample  collected  on  24  July  1982  was  also 
evaluated  for  viable  cells  by  bacterial  plate  counting.  Because  the 
lagoon  water  temperature  had  risen  above  35°C,  plates  were  incubated 
at  37 °C  and  25 °C.  The  number  of  bacterial  colonies  formed  at  each 
temperature  is  shown  in  Table  24.  These  results  indicate  that  more 
bacteria  will  grow  at  25 °C  than  at  37 °C.  However,  the  bacteria  that 
do  not  grow  at  37°C  are  still  viable. 

Table  24 


BACTERIAL  PLATE  COUNT  OF  LAAP  LAGOON  WATER 
COLLECTED  IN  JULY 


Incubation 

Temperature 


Average  Colonies  Formed 
(CFU  ml-1) 


2.60  x  10° 


37®C  1.19  x  106 

These  counts  do  not  differ  significantly  from  those  found  in 
December  1981  and  May  1982,  thus  suggesting  that  the  microbial  pop¬ 
ulation  within  the  lagoon  remains  fairly  constant  although  populations 
of  specific  organisms  may  change. 


2.  Biotransformation  Screening  Tests 


Because  of  the  presence  of  substantial  microbial  populations  in 
the  water,  we  initiated  screening  tests  for  the  biotransformation  of 
TNT  and  RDX  in  lagoon  water  under  aerobic  and  anaerobic  conditions. 

Lagoon  waters  from  the  four  sites  (north,  south,  east  and  west 
quadrants)  were  mixed  in  equal  parts  and  three  preparations  of  2-liter 
samples  were  each  placed  in  4-liter  bottles.  The  three  preparations 


were  1)  water  alone,  2)  water  plus  SO  ppm  sterile  yeast  extract,  and 
3)  water  plus  IX  dry  weight  of  bottom  sediment.  For  the  anaerobic 
incubation,  1  liter  of  each  of  the  above  solutions  was  placed  in  a  1- 
liter  Erlenmeyer  flask,  which  was  flushed  with  N2  gas  and  sealed  with 
a  rubber  stopper.  Flushing  with  N2  gas  was  repeated  after  each  sam¬ 
pling.  The  samples  were  analyzed  periodically  for  total  organic  car¬ 
bon  (TOC),  TNT,  and  RDX.  The  initial  concentrations  were  75  ppm  TOC, 
5.0  ppm  TNT,  and  13.0  ppm  RDX.  No  additional  TNT  or  RDX  was  added  to 
the  water. 

a.  Aerobic  TNT  Biotransformation 

After  1  day  of  incubation,  the  bacterial  population  in  the  lagoon 
water  under  aerobic  conditions  increased  slightly  to  3  *  10°  CFU  ml-1, 
then  gradually  decreased  to  1.2  x  10®  CFU  ml-*  after  16  days  of  incu¬ 
bation,  and  remained  constant  at  l.l  x  l6  CFU  ml-1  from  34  days  to 
70  days.  The  TOC  remained  constant  (75  ppm)  over  this  period. 

In  the  lagoon  water  alone,  the  TNT  concentration  slowly  decreased 
from  5  ppm  to  2  ppm  over  90  days  of  Incubation  (Figure  10).  The 
transformation  followed  first-order  kinetics.  The  pseudo-first-order 
rate  constant  was  calculated  to  be  4.1  *  10  4  hr  *,  and  the  second- 
order  rate  constant  was  calculated  to  be  3.7  x  10  10  ml  cell  1  hr  *, 
based  on  microbial  counts  of  l.l  x  10®  CFU  ml-*. 

The  addition  of  yeast  extract  to  the  lagoon  water  accelerated  the 
transformation  rate  by  increasing  the  microbial  population.  The 
microbial  count  increased  from  0.3  to  1.6  x  107  CFU  ml  *  after  one  day 
of  incubation.  However,  it  decreased  rapidly  to  1.6  x  10®  CFU  ml-* 
after  16  days,  and  the  transformation  rate  also  decreased.  The  rate 
accelerated  again  when  yea^t  extract  was  added  at  Day  41  (Figure  10). 

In  the  lagoon  water  with  added  sediment,  the  total  concentration 
of  TNT  increased  to  74  ppm  because  the  bottom  sediment  contained  large 
amounts  of  solid  TNT.  Total  TNT  in  thi*j  water  decreased  from  74  ppm 
to  59  ppm  during  the  90  days  of  incubation  in  a  similar  fashion  to 
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FIGURE  10  AEROBIC  BIOTRANSFORMA1  ION  OF  TNT  AND  RDX  IN  LAGOON  WATER 

•  TNT  in  water  alone;  *  RDX  in  water  alone;  O  TNT  in  water  plus  yeast  extract; 
a  ROX  in  water  plus  yeast  extract;  f  50  ppm  yeast  extract  added. 


that  observed  In  lagoon  water  alone.  Bacterial  counts  averaged  3  x 
106  CFU  ml-1  in  the  water- sediment  mixture. 

b.  Anaerobic  TNT  Biotransformation 

Under  anaerobic  conditions  In  the  lagoon  water  alone,  TNT  was 
transformed  at  a  slower  rate  than  under  aerobic  conditions  (Figure 
11).  This  may  be  a  result  of  a  lower  population  of  anaerobic  micro¬ 
organisms.  The  pseudo-f irst-order  rate  constant  was  calculated  to  be 
1.6  x  10'4  hr-1. 

The  addition  of  yeast  extract  to  the  water  also  accelerated  the 
disappearance  of  TNT.  Yeast  extract  was  added  at  Day  0,  Day  41,  and 
Day  76;  there  was  a  significant  drop  in  TNT  concentration  after  each 
addition. 


FIGURE  11  ANAEROBIC  BIOTRANSFORMATION  OF  TNT  AND  RDX  IN  LAGOON  WATER 

•  TNT  in  water  alone;  *  RDX  in  water  alone;  O  TNT  in  water  plus  yeast  extract; 
a  RDX  in  water  plus  yeast  extract;  |  yeast  extract  added. 


In  the  sediment-added  lagoon  water,  TNT  decreased  from  74  ppm  to 
57  ppm  during  76  days  of  Incubation.  Therefore,  It  appears  that  ana¬ 
erobic  biotransformations  in  anaerobic  waters  are  similar  to  the  aer¬ 
obic  biotransformation  of  TNT. 


c.  Aerobic  TNT  Biotransformation  at  High  and  Low 
Concentrations  of  TNT 

The  biotransformation  of  TNT  was  investigated  at  high  and  low 
concentrations  of  TNT  and  in  the  presence  of  a  high  microbial  pop¬ 
ulation.  The  lagoon  water  (stored  at  <  10°C)  was  exposed  to  room 
temperature  for  one  week,  centrifuged,  and  resuspended  in  one-tenth  of 
the  volume  of  lagoon  water  in  two  bottles.  Pure  TNT  was  added  to  one 
bottle  to  bring  the  aqueous  concentration  to  ~50  ppm  without  changing 
the  other  chemical  components  of  the  lagoon  water.  The  other  bottle 
contained  the  original  TNT  concentration  of  lagoon  water  (A. 2  ppm). 

The  average  bacterial  counts  in  these  waters  were  originally  7.8  x  106 
CFU  ml-1,  but  they  decreased  to  4.0  x  106  CFU  ml-1  after  25  days  of 
incubation. 

The  TNT  in  the  water  with  a  low  concentration  of  TNT  transformed 
with  a  pseudo-f irst-order  rate  constant  of  1.9  x  IO”3  hr*1  (4.6  x  10  2 
day"1;  correlation  coefficient,  r2  -  0.98).  The  second-order  rate 
constant  was  calculated  to  be  3.2  x  IO-10  ml  cell  1  hr  1  (7.7  x  10  9 
ml  cell-1  day-1),  based  on  a  median  bacterial  count  of  6.0  x  106  CFU 
ml”1.  This  rate  constant  is  close  to  that  observed  in  the  original 
water  (3.7  x  io  10  ml  cell  1  hr  1). 

The  TNT  in  the  water  with  a  high  concentration  of  TNT  decreased 
from  48  ppm  to  33  ppm  in  25  days.  If  the  data  are  treated  as  pseudo- 
first-order,  the  rate  constant,  calculated  by  a  least-squares  method, 
is  6.3  x  IO'"  hr-1  (1.5  x  10”2  day-1;  r2  «*  0.95).  This  rate  constant 
is  only  one-third  of  that  observed  in  the  water  with  a  low  TNT  con¬ 
centration.  If  the  data  are  treated  as  following  zero-order  kinetics, 
which  is  expected  for  high  chemical  substrate  biotransformation,  the 
zero-order  rate  constant  is  2.5  x  io”2  ppm  hr”1  (6.0  x  10”1  ppm  day  ^ 
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r  ■  0.96).  Because  the  correlation  coefficients  are  nearly  equal, 
more  data  points  are  required  to  prove  the  zero-order  kinetic 
behavior.  However,  since  the  calculated  pseudo-first-order  rate 
constant  is  much  lower  than  that  observed  in  the  low-TNT-concentration 
experiment,  a  zero-order  kinetic  process  may  better  describe  the 
transformation  rate  at  high  TNT  concentrations  (such  as  water  in  the 
sediment).  The  calculated  rate  constant  for  6  x  106  cell  ml-1  is  4.1 
x  10”9  ppm  cell”1  hr”1  (9.8  x  10“8  ppm  cell”1  day”1). 
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d.  TNT  Metabolitles 

During  the  laboratory  biotransformation  study  of  TNT  in  lagoon 
water,  we  observed  that  the  concentrations  of  2-araino-4,6-dlnitro- 
toluene  (2-A-4,6-DNT)  and  4-amino-2 , 6-dinitrotoluene  (4-A-2 ,6-DNT) 
remained  nearly  constant  throughout  the  study.  To  investigate  whether 
2-A-4.6-DNT  and  4-A-2.6-DNT  are  products  of  lagoon  water  biotrans- 
formatlon,  we  centrifuged  lagoon  water  and  resuspended  the  pre¬ 
cipitate-containing  microorganisms  in  l/20th  the  volume  of  lagoon 
water.  Ring-labeled  ^C-TNT  (1.0  pCi  in  50  ml)  and  unlabeled  TNT 
(final  concentration  40  ppm)  were  added  to  the  water  and  the  water  was 
Incubated  in  the  dark. 

After  3  days  and  19  days  of  incubation,  the  water  was  extracted 
with  ethyl  acetate,  then  acidified  to  pH  2,  and  re-extracted  with 
ethyl  acetate.  In  the  3-day  sample,  more  than  90%  of  the  radio¬ 
activity  was  found  in  the  neutral  fraction.  The  extracts  were  dried, 
concentrated,  and  chromatographed  by  TLC  using  silicia  gel  plates  and 
chloroform  as  the  developing  solvent.  The  plates  were  analyzed  by 
radioautography.  Most  of  the  activity  was  associated  with  TNT; 
however,  radioactive  2-A-4,6-DNT  and  4-A-2,6-DNT  were  produced  and 
identified  by  their  yellow  color  and  by  co-chromatography  against 
authentic  standards. 


In  the  19-day  sample,  the  distribution  of  the  radioactivity  among 
the  neutral,  the  acid  extract,  and  the  remaining  aqueous  phase  was 
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80X,  16%,  and  4%,  respectively.  The  radioautograph-TLC  plate  deve¬ 
loped  In  chloroform,  In  benzene:ethyl  acetatetacetlc  acid  (6:3:1),  and 
In  chloroform:methanol  (6:1)  solvent  systems  showed  the  major  activity 
at  spots  with  values  corresponding  to  TNT,  4-araino-2,6-dinitro- 
toluene,  and  2-amlno-4,6-dinltrotoluene  (Figure  12).  The  activity  at 
s pots  corresponding  to  2,4-diaraino-6-nltrotoluene  and  2 ,6-diaraino-4- 
nltrotoluene  were  Insignificant.  Several  minor  spots  with  weak 
activity  had  Rf  values  similar  to  those  for  ti initrobenzyl  alcohol  and 
trinitrobenzaldehyde.  However,  further  identification  is  needed.  The 
acid  extract  showed  several  radioactive  components  with  low  values. 


Chloroform 


Benzeneiethyl  acetate: acetic  acid 
(6:3:1) 


Chloroform :  methanol 
(6:1) 


LA-7934-37 

FIGURE  12  RADIOAUTOGRAPH  OF  THIN-LAYER  CHROMATOGRAM  OF  TNT 

METABOLITES  IN  LAGOON  WATER  AFTER  19  DAYS  OF  INCUBATION 


Samples;  N  =  neutral  extract;  A  =  acidic  extract;  r  *=  reference  chemicals. 


Compounds:  1  ■  TNT;  2  *  2,2',6,6'-tetranitro-4-8Zoxytoluene;  3  ■  4-amino- 
2.6-dinitrotoluene;  4  ■  2-amino-4,6-dinitrotoluene;  5  *  2,6- 
dinitroaniline;  6  =  2,4-dinitroaniline;  7  =  trinitrobenzyl 
alcohol;  8  =  trinitrobenzaldehyde;  9  *  1,3,6-trinitrobenzene; 

10  ■  picric  acid. 
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The  extracts  were  streaked  on  preparative  TLC  plates.  Several 
unknown  compound  spots  were  scraped  from  the  plate  and  subjected  to 
probe  mass  spectrometrlc  analysis;  however,  they  could  not  be  iden¬ 
tified  because  of  large  amounts  of  Interfering  phthalates  and  binders 
from  the  TLC  plates. 

On  the  basis  of  this  study,  we  conclude  that  the  formation  of 
amlnodlnitrotoluenes  In  lagoon  water  is  the  result  of  microbial 
processes. 

e .  Biological  Oxygen  Demand  of  the  Lagoon  Water 

The  biological  oxygen  demand  of  the  LAAP  lagoon  water  was  mea¬ 
sured  as  an  indicator  of  the  oxygen-utilization  rate  of  the  mic¬ 
robes.  The  refrigerated  lagoon  water  was  left  at  room  temperature 
overnight,  then  shaken  to  saturate  the  water  with  air,  and  poured  Into 
four  biological  oxygen  demand  (BOD)  test  bottles,  which  were  then 
Incubated  at  20°C  in  the  dark.  The  dissolved  oxygen  (DO)  concen¬ 
tration  was  monitored  with  a  DO  meter  after  0,  3,  5,  7,  and  10  days  of 
incubation.  Thu  DO  decreased  at  each  sampling  point,  yielding  a 
first-order  rate  constant  of  9.6  x  lO-*  day  We  therefore  expect 
that  in  the  LAAP  lagoon  (DO  **  8  ppm)  the  microbes  will  use  oxygen  at  a 
rate  of  0.77  ppm  day-1. 


f .  RDX  Aerobic  Blotransformatlon 

The  original  concentration  of  RDX  lr.  the  lagoon  water  was  13 
ppm.  No  loss  of  RDX  was  observed  during  the  aerobic  experiment — 
whether  RDX  was  in  the  lagoon  water  alone,  with  added  yeast  extract, 
or  with  1%  of  bottom  sediment  (Figure  10).  Therefore,  it  appears  that 
RDX  is  biologically  persistent  in  the  lagoon  water  under  aerobic  con- 
di t  ions . 


g .  RDX  Anaerobic  Blotran3formatlon 


Under  anaerobic  conditions,  no  change  in  RDX  concentration  was 
observed  in  the  lagoon  water  alone  or  in  water  plus  sediment.  RDX  (13 
ppm)  in  lagoon  water  with  50  ppm  of  added  yeast  extract  did  not  show 
any  loss  of  RDX  after  46  days.  An  additional  50  ppcn  of  yeast  extract 
was  added  on  Day  46  and  100  ppm  on  Day  76.  After  83  days  of  incu¬ 
bation,  the  RDX  concentration  remained  at  13  ppm  and  then  dropped  to 
2.9  ppm  at  Day  90  and  1.4  ppm  at  Day  92.  t  ->rhaps  the  repeated 
addition  of  yeast  extract  acclimated  RDX-utilizing  organisms. 

On  Day  92,  the  microorganisms  in  the  lagoon  water  were  inoculated 
into  basal  salts  medium  (BSM)  plus  500  ppm  of  yeast  extract  and  9.1 
ppm  of  RDX.  The  RDX  concentration  decreased  to  0.6  ppm  after  5  days 
of  anaerobic  incubation.  The  organisms  were  transferred  to  a  second 
flask  containing  BSM,  50  ppm  yeast  extract,  and  9.4  ppm  of  RDX;  the 
RDX  was  transformed  to  3.4  ppm  after  4  days  of  incubation. 

The  RDX-acclimated  organisms  were  then  inoculated  into  two 
flasks:  one  contained  BSM  plus  50  ppeu  yeast  extract  and  8.0  ppm  RDX, 
the  other  contained  lagoon  water  (12.9  ppm  RDX)  and  50  ppm  of  yeast 
extract.  The  results  of  the  anaerobic  incubation  are  shown  in  Table 
25. 


Table  25 

RDX  CONCENTRATION  ON  DIFFERENT  DAYS  AFTER 
ANAEROBIC  INCUBATION  IN  LAAP  LAGOON  WATER 


Day 

RDX 

(ppm) 

0 

12.4 

5 

9.8 

10 

5.8 

17 

2.5 

20 

<0.1 

V  TV 


ID 


When  the  above  organisms  were  inoculated  into  BSM  with  8.0  ppm  RDX  and 
50  ppm  of  yeast  extract,  the  RDX  was  transformed  to  nondetectable 
levels  (<0.01  ppm)  within  5  days. 

In  a  second  experiment,  the  above  organisms  were  inoculated  into 
samples  of  lagoon  water  with  50  ppm  or  300  ppm  yeast  extract.  A  third 
flask  was  inoculated  with  RDX-transforming  organisms  [obtained  from 
the  HMX  waste-line  sediment  at  Holoton  Army  Ammunition  Plant  (HAAP)], 
lagoon  water,  and  50  ppm  of  yeast  extract.  The  RDX  concentration  was 
monitored  over  a  20-day  period.  The  results  (Table  26)  indicate  that 
the  anaerobic  biotransformation  of  RDX  can  be  accelerated  by  extra 
organic  nutrients  (300  ppm  yeast  extract)  or  by  a  different  source  of 
RDX-transformlng  organisms. 

Table  26 

ANAEROBIC  TRANSFORMATION  OF  RDX  WITH  SELECTED 
ORGANISMS  AND  YEAST  EXTRACT  CONCENTRATIONS 


Source  of 


Yeast  Extract  RDX  Concentration  (ppm)  at  Day 


m  . 

Organisms 

0 

7 

13 

17 

20 

> 

LAAP  Lagoon  Water 

50 

13.0 

11.3 

7.0 

5.3 

3.6 

LAAP  Lagoon  Water 

300 

12.7 

11.1 

8.4 

7.8 

<0,1 

« 

HAAP  Sediment 

50 

12.6 

11.2 

4.8 

1.0 

<0.1 

For  comparison,  RDX-tranforming  organisms  from  the  -AAP  water  and 
the  HAAP  sediment  were  each  inoculated  into  flasks  containing  BSM  with 
9.0  ppm  RDX  and  50  ppm  yeast  extract;  the  RDX  disappeared  within  three 
days.  This  result  and  the  data  in  Table  26  show  that  the  anaerobic 
biotransformation  of  RDX  is  much  slower  in  lagoon  water  than  in  BSM 
even  in  the  presence  of  large  yeast  extract  concentrations  (300  ppm) 
or  with  an  alternate  source  of  RDX-transfonning  organisms.  These  data 
suggest  that  lagoon  water  components  have  an  inhibitory  effect  on  RDX 
b lot rana format  ion. 
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To  test  the  possibility  that  the  slow  RDX  anaerobic  biotrans- 
formatlon  in  lagoon  water  was  due  to  the  lack  of  mineral  nutrients 
rather  than  to  inhibition  by  some  chemical  present  in  the  lagoon 
water,  we  investigated  the  biotransformation  of  RDX  in  lagoon  water 
with  yeast  extract  and  BSM  components.  The  results  showed  that  the 
addition  of  mineral  salts  to  lagoon  water  collected  either  in  December 
or  May  did  not  significantly  accelerate  the  RDX  biotransformation, 
suggesting  that  the  lack  of  salts  is  not  responsible  for  the  inhibi¬ 
tory  effect  of  RDX  anaerobic  biotransformation  in  lagoon  water. 

3.  Detailed  Biotransformation  Rate  Study 

The  biotransformation  of  TNT  depends  on  enzymes  resulting  from 

microorganisms  utilizing  organic  nutrients  other  than  TNT  for 

growth.  Under  these  conditions,  the  biotransformation  rate  will 

follow  Mlchaelis-Menton  kinetics  at  different  concentrations  of  TNT. 

dc 

The  transformation  rate,  — ,  can  be  represented  as  shown  in  Equation 
35. 


jf  ’  u  '  VmC/(K»  +  C>  -  (35) 

where  C  is  the  TNT  concentration,  V„  is  the  maximum  transformation 
rate  at  a  given  enzyme  concentration  or  microbial  population,  and  Km 
is  the  Michaelis  constant,  which  is  the  TNT  concentration  that  yields 
half  the  maximum  transformation  rate. 

When  the  TNT  concentration  is  very  low,  Km  is  much  greater  than  C 
and  the  transformation  rate  expression  reduces  to  Equation  36: 


dC 

dt 


(36) 


which  is  the  pseudo-first-order  transformation  rate  expression. 
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In  environments  where  the  TNT  concentration  Is  high,  such  as  near 
the  sediment  in  the  LAAP  lagoon,  C  is  much  greater  than  K  and 
Equation  35  reduces  to  Equation  37: 


dC 

■  V 

dt  m 


(37) 


When  V  is 
m 


which  is  a  zero-order  rate  expression  with  respect  to  TNT 
divided  by  the  bacterial  concentration,  the  maximum  rate  per  organism 
can  be  obtained. 


Equation  35  can  be  rewritten  as  shown  in  Equation  38,  which  is 
the  Lineweaver-Burke  equation.  A  plot  of  1/u  versus  1/C  yields  a 
straight  line,  from  which  the  values  Km/Vm  (slope)  and  1/Vm  (y-inter- 
cept)  can  be  determined. 


I 
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K  i 

JS-.+ 

VC  V 
m  m 


(38) 


To  determine  Km  and  Vm,  we  Incubated  LAAP  lagoon  water  with  100 
ppm  of  yeast  extract  for  two  days  to  generate  a  high  cell  popu¬ 
lation.  The  water  was  dispensed  into  five  flasks  and  TNT  was  added  to 
give  concentrations  ranging  from  2.39  to  37.7  ppm.  TNT  was  monitored 
for  24  hr  and  the  transformation  rates  were  determined.  These  data 
are  shown  in  Table  27. 


Using  the  Lineveaver-Burke  method,  we  plotted  the  above  data 
(Figure  13)  and  calculated  Kra  and  Vra.  Linear  regression  of  the  plot 
gave  values  of  Km  “  19.3  ppn  and  Vm  *  3.5  ppm  day-1.  The  initial 
microbial  population  was  1.0  x  108  CFU  ml  1  and,  therefore,  the 
maximum  rate  per  organism  will  be  3.5  x  1O-0  ppm  cell-1  day”1  (1.5  * 

10  9  ppm  cell  1  hr  1).  The  value  is  slightly  lower  than  the  9.8  x 
10  8  ppm  cell  1  day  1  obtained  with  centrifuged  lagoon  water  mico- 
organisms  and  50  ppm  TNT  determined  in  a  previous  experiment.  This 
may  be  the  result  of  a  modification  of  the  mixed  culture  population 
when  grown  on  yeast  extract. 

The  above  results  indicate  that  the  biotransformation  of  TNT  will 
follow  pseudo-first-order  kinetics  at  TNT  concentrations  less  than  19 
ppm  and  zero-order  kinetics  at  TNT  concentrations  greater  than  19  ppm. 

A  second  experiment  was  conducted  with  lagoon  water 
microorganisms  Incubated  without  the  addition  of  yeast  extract.  From 
the  lagoon  water  collected  in  May,  we  prepared  high  cell  populations 
by  incubating  the  water  at  20°  to  25°C  in  the  dark  over  a  three-week 
period,  centrifuging  the  microorganisms,  and  resuspending  the  cells  in 
l/35th  of  the  volume  of  lagoon  water.  This  water,  containing  2.1  x 
107  CFU  ml  1 ,  was  dispensed  into  flasks  and  TNT  was  addsj  to  give 
concentrations  ranging  from  8  to  47  ppm. 


The  biotransformation  of  TNT  was  followed  over  five  days.  The  initial 
TNT  concentrations  and  biotransformation  rates  are  shown  in  Table  28. 


Table  28 

INITIAL  TNT  CONCENTRATIONS  AND  BIOTRANSFORMATION 
RATES  OVER  5  DAYS 


Initial  TNT  Concentration 
_ (ppa) _ 

8.04 

10.7 

16.3 

21.3 

34.7 

47.3 


Biotransformation  Rate 
(ppm  day-1) _ 

0.33 

0.40 

0.51 

0.55 

0.61 

1.02 


Using  the  Lineweaver-Burke  plotting  method,  we  calculated  Km  (the 
Michaells  Constant)  to  be  19.1  ppm  and  Vm  (the  maximum  transformation 
rate  at  a  given  enzyme  concentration  or  microbial  population)  to  be 
1.1  ppm  day  *.  The  maximum  rate  per  organism  is  5.3  *  ID-0  ppm  cell"1 
day'1  (2.2  x  10  9  ppra  cell”1  hr”1),  which  is  close  tc-  the  value 
obtained  for  the  yeast  extract-generated  lagoon  wacer  organisms. 


VI.  DEVELOPMENT  OF  THE  SIMULATION  MODEL 


A.  General  Equations 

A  simulation  model  was  developed  based  on  a  mathematical  des¬ 
cription  of  material  balance.  The  basic  equation  for  material  balance 
is  shown  in  Equation  39: 


dM 

dt 


(39) 


where  M  *  total  amount  of  chemical  in  the  system 
1^  »  total  input  of  chemical  to  the  lagoon 
Oj  ■  total  removal  of  chemical  from  the  lagoon. 

The  change  in  total  mass,  dM/dt,  can  also  be  related  to  the  chemical 
concentration,  0,  and  volume,  V,  of  the  lagoon  as  shown  in  Equation 
40: 


dM 

dt 


(40) 


Thus,  from  Equations  39  and  40,  Equation  41  can  be  formulated: 


(41) 


Becaus '  no  munitions  wastes  have  been  dumped  since  early  in  1981, 
the  only  inpit  source  of  munitions  arises  from  the  diffusion  of  accum¬ 
ulated  chemicals  in  the  sediment  layer  (D<jif)  and  the  dissolution  of 
solid  chemicals  lying  on  the  bottom  sediment  surface  (Dj^g).  There¬ 
fore,  the  chemical  input  can  be  described  as  shown  in  Equation  42: 


(42) 
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Because  the  bottom  sediment  and  surface  are  heterogeneous  with 
respect  to  the  amounts  of  TNT  and  RDX  present,  the  bottom  sediment  was 
divided  Into  25  compartments  (30*  x  30')  from  which  the  diffusion  flux 
and  dissolution  rate  were  calculated  Individually  (1  -  1  to  25). 


Chemicals  are  removed  from  the  lagoon  through  transformation 
processes  (photolysis  and  biotransformation)  and  seepage.  The  seepage 
transport  process  has  been  accounted  for  in  the  calculation  of  net 
diffusion  flux  (D^if);  therefore  the  total  removal  of  chemical  can  be 
expressed  as  Equation  43: 


E  0  -  Ek  CV  ,  (43) 

j  J  * 


where  represents  the  transformation  rate  constants. 
Equation  41  can  now  be  written  as  Equation  44: 
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A  numerical  integration  scheme,  using  Euller's  prediction  and 
correction  method,  was  used  to  integrate  Equation  44  over  selected 
time  intervals  to  provide  simulation  results. 


B.  Program  Structure 

A  flow  chart  of  the  computer  program  is  shown  in  Figure  14. 

System  parameters,  transformation  rate  constants,  weather  data,  lagoon 
descriptions,  and  initial  conditions  are  read  in  when  the  program  is 
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FIGURE  14  SIMULATION  FLOWCHART  OF  THE  LAGOON  MODEL 


initiated.  The  main  program  is  linked  to  a  number  of  subroutines  such 
as  VOLUME  and  SURFACE,  which  calculate  the  depth  above  each  sediment 
compartment  and  the  number  of  compartments  covered  by  water.  This 
number  and  the  corresponding  locations  of  sediment  covered  by  water 
are  then  used  separately  in  FLUX  and  FLUX  1  subroutines  to  calculate 
the  diffusion  flux  and  dissolution  rate  of  chemicals  from  the  bottom 
sediment,  respectively.  Material  balance  is  performed  in  the  sub¬ 
routine  COMPS  when  Euller's  prediction  and  correction  method  is 
employed  to  predict  the  chemical  concentration  in  the  lagoon.  The 
above  procedures  are  repeated  at  every  time  step  until  the  simulation 
period  is  complete  or  the  lagoon  goes  dry.  The  status  of  the  lagoon 
is  obtained  at  specified  print-out  time  intervals,  and  a  plotting 
subroutine  GRAPH  4  is  used  to  plot  the  simulation  results.  The 
program  Interactions  are  shown  in  Figure  15,  and  a  complete  listing  of 
the  program  appears  in  Appendix  C. 

C.  Simulation  Results 

1.  Simulation  of  the  Water  Surface 

The  effect  of  seepage  rate  on  the  lagoon  depth  was  discussed  in 
Section  V.A.l;  an  estimated  rate  of  0.12  cm  day-1  gave  good  agreement 
with  measured  values.  Changes  in  the  water  surface  were  also  simu¬ 
lated,  as  shown  in  Figure  16.  As  the  lagoon  dries,  the  southeast  cor¬ 
ner  is  the  first  to  dry  out  and  the  lagoon  recedes  to  the  northwest 
corner.  Observations  of  the  water  remaining  in  the  lagoon  in  July  and 
August  were  in  good  agreement  with  the  model  simulation. 
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FIGURE  15  PROGRAM  INTERACTIONS 
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FIGURE  16  APPROXIMATE  WATER  CONTOUR  LINES  DURING  THE  SIMULATION 
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2. 


TNT  Simulation 


a.  Case  1:  Complete  Mixing 

Five  scenarios  were  investigated  in  which  the  assumption  was  made 
that  the  lagoon  is  a  completely  mixed  reactor.  Figure  17A  shows  the 
TNT  concentration  as  a  function  of  time  using  inputs  from  diffusion 
flux  and  dissolution,  a  photolysis  rate  constant  of  0.78  cm  day  *,  and 
a  biotransformation  rate  constant  of  0.012  day**.  The  results  indi¬ 
cate  that  the  TNT  concentration  will  increase  until  Day  200,  with  the 
peak  concentration  being  17.5  ppm.  After  Day  200,  TNT  concentration 
decreases  due  to  decreases  in  the  lagoon  depth,  which  increases  the 
magnitude  of  the  photolysis  rate  constant. 

To  test  the  impact  of  the  photolysis  rate  constant,  a  simulation 
was  run  under  conditions  identical  to  those  above  except  that  the 
photolysis  rate  constant  was  reduced  by  half — to  0.39  cm  day  *. 

Figure  17B  shows  the  simulation  result.  In  this  case,  the  TNT  con¬ 
centration  builds  up  faster  than  in  the  previous  case  and  reaches  a 
maximum  of  23  ppm  on  Day  200.  After  Day  200,  the  TNT  concentration 
decreases  ac  a  slower  rate  than  before. 

The  photolysis  rate  constant  is  also  seasonally  dependent  due  to 
high  light  Intensity  in  the  summer  months  relative  to  that  in  the 
winter  months.  The  magnitude  of  the  intensity  change  is  usually 
between  a  factor  of  2  to  3.  In  Figure  17C,  we  adjusted  the  photo¬ 
chemical  rate  constant  by  a  factor  of  2,  starting  in  June,  and 
observed  an  accelerated  decrease  in  TNT  concentration  in  the  late 
summer  months.  The  overall  effect  appears  to  be  a  lowering  of  the 
maximum  TNT  concentration  buildup  followed  by  a  rapid  loss  of  TNT. 

We  performed  a  simulation  in  which  the  biotransformation  rate 
constant  was  doubled  for  every  10®C  rise  in  lagoon  temperature,  based 
on  a  number  of  assumptions  involving  the  enzymatic  reduction  of  a  TNT 
nitro  group.  The  result  is  plotted  in  Figure  18.  In  the  early  part 
of  the  simulation,  the  TNT  concentration  rises  above  that  observed  in 
Figure  17a  due  to  the  low  biotransformation  rate  constant.  However, 
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FIGURE  17  SIMULATED  TNT  CONCENTRATION  WITH  RESPECT  TO  TIME  OF  YEAR 

Curve  A:  kp  =  0.78  cm  day'1,  k1  -  0.012  day'1;  input  from  diffusion  and  dissolution. 

Curve  B:  kp  =  0.39  cm  day’1,  k£  =  0  012  day'1,  input  from  diffusion  and  dissolution. 

Curve  C:  kp  =  0.78  cm  day'1  plus  radiant  correction;  kp  =  0.12  day'1;  input  from 

diffusion  and  dissolution. 


the  effect  of  biotransformation  becomes  evident  between  Days  ISO  and 
200.  In  the  latter  phase  of  the  simulation,  photolysis  becomes  the 
dominant  transformation  process. 

Finally,  we  considered  a  case  in  which  dissolution  was  the  only 
input  source  of  TNT  (Figure  19).  The  simulation  results  show  TNT 
concentration  decreasing  in  the  early  stage  of  simulation  and  then 
increasing  in  May,  reaching  a  maximum  at  Day  200  and  then  decreasing 
as  the  lagoon  dries  out. 

The  concentration-tirae-dependent  profile  measured  in  the  lagoon 
for  TNT  appears  in  Figure  20.  None  of  the  above  scenarios  exactly 
match  the  measured  data.  The  low  concentration  observed  in  March  (Day 
74)  does  not  appear  in  any  simulation  except  in  the  case  of  the  sole 
input  arising  from  dissolution  (Figure  19).  The  maximum  TNT 
concentration  was  observed  in  May  (Day  145),  whereas  the  simulation 
results  predicted  that  the  maximum  would  occur  in  July.  To  reach  a 
maximum  in  May,  a  large  input  from  the  bottom  sediment  would  be 
expected.  The  present  model  probably  overestimates  the  diffusion  and 
dissolution  processes  in  the  latter  months.  Also,  the  lagoon  begins 
to  dry  out  at  the  end  of  June;  this  process  should  concentrate  chem¬ 
icals  due  to  volume  reduction  and  concentrate  the  TNT  in  July  unless 
the  transformation  processes  are  greatly  accelerated.  Both  the 
photolysis  and  biotransformation  of  TNT  in  lagoon  water  collected  in 
July  give  no  suggestion  of  dramatic  rate  increases  compared  to  other 
experimental  data. 

The  discrepancy  between  the  observed  data  and  the  simulations  may 
be  explained  by  the  sampling  technique  and  the  lagoon  location  from 
which  the  samples  were  taken.  Because  the  plant  personnel  did  not 
have  access  to  th^  center  of  the  lagoon  to  collect  samples,  all  sam¬ 
pling  done  by  them  was  performed  at  the  southeast  edge  near  the  dump 
ramp,  where  high  levels  of  TNT  were  found  in  the  sediment.  As  the 
lagouu  dried,  samplet,  were  taken  progressively  closer  to  the  center  of 
the  lagoon  approaching  the  northwest  corner.  Poor  mixing  of  the 
lagoon  water  may  therefore  obscure  the  simulation  results. 
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FIGURE  19  SIMULATED  TNT  CONCENTRATION  FROM  THE  CASE  OF  SOLE 
INPUT  FROM  DISSOLUTION 
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Case  II:  Zero  Mixing 


If  the  sampling  locations  become  critical,  the  assumption  of  com¬ 
plete  mixing  may  not  be  a  valid  representation  of  the  distribution  of 
chemicals  in  the  lagoon.  Although  the  data  on  dissolved  oxygen 
obtained  in  December  indicated  that  the  lagoon  was  well-mixed  in  the 
vertical  direction,  discrepancies  in  TNT  concentrations  at  the  edge  of 
the  lagoon  relative  to  the  center  were  noted.  This  suggests  that  the 
horizontal  mixing  of  lagoon  water  may  not  be  complete  and  may 
intensify  as  the  lagoon  becomes  shallow. 

We  performed  simulations  for  the  case  of  zero-mixing  for  com¬ 
partments  1  (dump  site)  and  13  (center).  In  the  case  of  zero-mixing, 
the  lagoon  waters  do  not  interact  with  neighboring  compartments — only 
with  the  bottom  sediment  directly  below  the  water.  The  simulation 
results  are  shown  in  Figure  21.  In  compartment  l,  the  TNT  concen¬ 
tration  continuously  increases  until  maximum  solubility  is  reached. 

In  compartment  13,  the  TNT  concentration  continually  decreases  due  to 
low  TNT  levels  in  the  sediment  and  transformation  rates  that  surpass 
those  of  diffusion  and  dissolution  fluxes.  In  the  real  environment, 
partial  mixing  probably  occurs  that  dissipates  the  TNT  from  compart¬ 
ments  l  to  other  compartments  in  the  lagoon.  These  results  do  imply, 
however,  that  sampling  from  the  southeast  corner  and  following  the 
water-coverage  profile  (Figure  16)  for  further  samples  may  give  skewed 
versions  of  the  TNT  concentration  profile  that  are  difficult  to 
model.  In  this  respect,  the  model  simulation  may  be  adequately 
predicting  the  loss  and  movement  of  TNT  in  the  dissolution-input-only 
case  (Figure  19),  when  compared  to  the  actual  data,  that  may  be  higher 
than  in  a  completely  mixed  reactor  (Figure  20).  An  overlay  of  the 
latter  profiles  appears  in  Figure  22. 
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FIGURE  21  SIMULATION  RESULTS  FOR  TNT  IN  COMPARTMENTS  1  AND  13 
FOR  THE  CASE  OF  ZERO-MIXING 
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3.  RDX  Simulation 


The  laboratory  studies  showed  that  there  would  be  little  bio¬ 
transformation  of  RDX  in  the  lagoon  and  that  photolysis,  although 
slow,  probably  would  be  the  dominant  transformation  process.  Knowing 
the  range  of  the  depth-independent  photochemical  rate  constant  for  the 
winter  through  summer  seasons,  we  derived  a  linear  expression  relating 
the  rate  constant,  kp,  to  the  day  of  the  year,  D: 

_2  _2  _1 
Kp  =*  (3.0  x  10  cm  day  )  (D)  +  1.6  x  10  cm  day 


(45) 

Thus,  at  Day  0  (January  1)  the  winter  rate  constant,  0.016  cm  day-1, 
is  applied  to  the  model  and  linearly  increases  to  Day  200  (July  19), 
where  the  summer  rate  constant  is  applied.  At  each  simulation  inter¬ 
val,  the  rate  constant  is  adjusted  according  to  the  depth  of  the 
lagoon  for  that  time  period. 

Three  scenarios  were  simulated  representing  the  observed  rate 
constants  and  rate  constants  representing  75%  and  50%  of  the  observed 
rate  constants.  The  latter  rate  constants  provide  a  sensitivity 
analysis  of  the  impact  of  the  photochemical  rate  constant  on  RDX  loss 
in  the  lagoon.  The  simulation  results  are  plotted  in  Figure  23.  In 
each  case,  the  simulated  results  indicate  that  the  removal  of  RDX  is 
balanced  by  the  diffusion  and  dissolution  input  over  the  first  150 
days.  As  the  water  gets  wanner,  the  concentration  increases  up  to  Day 
225;  then  the  effect  of  photolysis  is  noted  only  when  the  lagoon  is 
very  shallow.  Also,  halving  the  photolysis  rate  constant  (Curve  C) 
does  not  have  a  large  effect  on  the  PDX-loss  profile. 

Superimposed  on  these  curves  are  the  actual  measured  values  for 
RDX  in  the  lagoon.  Considering  the  problem  of  partial  mixing,  the 
simulation  results  and  the  actual  values  are  comparable  except  for  one 
data  point  where  a  very  high  level  (42.1  ppm)  of  RDX  was  observed. 
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FIGURE  23  SIMULATION  RESULTS  FOR  RDX  CONCENTRATION  IN  LAAP  LAGOON  WATER 

Curve  A:  k'  =  k  ;  Curve  B:  k'  =  0  75  k  ;  Curve  C  =  k'  =  0.50  k 
p  P  p  p  p  p 

•  =  Found  Value 


The  remaining  data  are,  however,  inconsistent  with  this  point.  A 
curve  for  the  profile  of  RDX  can  be  approximated  (dotted  line)  that 
shows  a  peak,  in  RDX  concentration  in  late  July  and  a  decrease  in  late 
August.  These  results  suggest  that  the  diffusion/dissolution  pro¬ 
cesses  become  faster  than  predicted  by  the  model  in  May  and  June  (Days 
145  to  175)  and  that  the  loss  rate  is  faster  than  predicted  in  late 
August  (Day  245).  This  is  probably  a  result  of  (a)  warmer  waters  near 
the  edge  of  the  lagoon  where  the  samples  were  taken,  which  could  give 
high  concentrations  in  May-June,  and  (b)  the  low  water  level  and  the 
recession  of  the  water's  edge  from  the  muni tions-enr lched  sediment 
area,  which  could  give  low  concentration?  in  late  August.  If  complete 
mixing  were  to  occur,  the  simulation  results  might  be  fairly  accurate. 


95 


VII.  DISCUSSION 


The  results  of  this  study  show  that  the  lagoons  at  LAAP  are 
highly  complex  environments  in  which  the  fate  of  TNT  is  controlled  by 
photolysis  and  biotransformation  and  the  fate  of  RDX  probably  19  con¬ 
trolled  by  photolysis.  If  there  were  no  chemical  inputs  to  the 
lagoon,  we  would  estimate  a  TNT  half-life  ranging  from  8  days  in  the 
winter  months  to  <4  days  in  the  summer  and  an  RDX  half-life  ranging 
from  173  days  in  the  summer  to  >1000  days  in  the  winter  in  a  lagoon  50 
cm  deep.  However,  there  are  chemical  inputs  to  the  lagoon  from  solid 
and  sorbed  munitions  in  the  bottom  sediment  and  on  the  bottom  sur¬ 
face.  Thus,  both  diffusion  and  dissolution  kinetics  play  Important 
role9  in  the  input  of  chemicals  to  the  lagoon.  These  processes  are 
complicated  by  1)  a  heterogeneous  distribution  of  munitions  in  and  on 
the  sediment  and  2)  a  temperature-dependence  that  Increases  chemical 
inputs  as  the  lagoon  warms.  Thus,  chemical  input  is  not  a  constant 
process,  but  one  that  can  change  dramatically,  depending  on  the  time 
of  year  and  the  portion  of  the  lagoon  that  is  covered  by  water. 

Furthermore,  each  transformation  process  has  additional  com¬ 
plications.  In  the  case  of  biotransformation,  the  enzymatic  reduction 
of  the  nltro  group  should  he  temperature-dependent,  and  therefore  the 
biotransformation  rate  constant  could  vary  by  a  factor  of  8  or  more 
(assuming  a  doubling  of  the  rate  constant  for  every  10°  rise  in  tem¬ 
perature)  between  the  winter  and  summer  months.  Also,  we  cannot  be 
sure  that  the  population  of  organisms  at  10°C  is  the  same  as  that  at 
40®C.  The  temperature  effect  could  have  an  impact  on  the  number  and 
type  of  microorganisms  that  are  available  to  transform  TNT.  To 
average  this  effect,  we  used  a  rate  constant  determined  at  25°C  as  an 
average  value  in  one  example  and  doubled  the  rate  constant  for  every 
10°C  rise  in  temperature  in  another  example.  The  results  of  these 


examples  indicated  that  changes  in  the  biotransformation  rate  constant 
did  not  have  a  major  impact  on  the  loss  profile  of  TNT  because  the 
photochemical  transformation  process  overwhelmed  any  such  changes. 


There  is  also  a  high  degree  of  uncertainty  in  the  photochemical 
rate  constant  as  a  function  of  the  time  of  year.  Photochemical  rate 
constants  measured  in  lagoon  water  collected  in  December  showed  an 
inverse  relationship  with  depth,  which  allows  the  prediction  of  a  rate 
constant  at  any  depth,  provided  that  the  absorptivity  of  the  water  is 
the  same.  However,  as  the  lagoon  waters  evaporate,  much  of  the  UV- 
absorbing  materials  become  concentrated,  reducing  the  ability  of  light 
to  penetrate  into  the  water.  The  effect  of  the  more  concentrated 
waters  will  be  to  lower  the  photochemical  rate  constant,  but  sim¬ 
ultaneously  the  shallower  depth  will  tend  to  increase  it.  Because  the 
photochemical  transformation  of  TNT  occurs  by  indirect  mechanisms  and 
is  governed  by  unknown  sensitizers  (both  humic  materials  and  TNT 
transformation  products)  in  the  waters,  it  is  difficult  to  correlate 
the  photochemical  rate  constant  dependences  on  the  absorptivity  and 
depth  of  the  lagoon  water  from  the  experiments  that  were  performed 
during  this  study. 

Sampling  was  a  problem  because  of  poor  access  to  a  center 
location  for  the  bulk  water  and  because  the  sites  from  which  the 
samples  were  collected  changed  as  the  lagoon  dried.  Finally,  the 
lagoon  dried  up  in  late  August. 

In  spite  of  these  difficulties,  we  were  able  to  develop  a 
simulation  model  that  describes  the  loss  profile  of  TNT  and  RDX  as  a 
function  of  time  of  year.  Of  primary  importance  in  this  model  is  the 
dissolution  rate  of  pure  chemicals  from  the  bottom  sediment.  This 
process  is  temperature-dependent  and  therefore  we  expect  the 
concentrations  of  TNT  and  RDX  to  rise  as  the  lagoon  warms  in  the 
spring  and  summer  months.  In  the  absence  of  transformation, 
dissolution  will  continue  until  the  solubility  limit  for  each  chemical 
is  reached. 
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A  concentration  buildup  in  the  lagoon  is  counteracted  by  the 
photochemical  and  biological  transformation  of  TNT  and  the  photo¬ 
chemical  transformation  of  RDX.  The  rate  of  photochemical  trans¬ 
formation  of  TNT  and  RDX  follows  pseudo-first-order  kinetics  and 
first-order  kinetics,  respectively.  Thus,  for  a  given  depth,  TNT  will 
photolyze  approximately  50  times  faster  than  RDX.  The  rate  constant 
is  dependent  on  the  solar  radiance;  therefore,  rates  will  be  two  to 
three  times  faster  in  the  summer  than  in  the  winter  months  for  TNT  and 
up  to  five  times  faster  for  RDX. 

The  photochemical  rate  constant  is  inversely  proportional  to 
depth;  hence,  rates  of  transformation  will  Increase  as  water  is  lost 
to  evaporation  and  seepage.  Usually,  the  lagoon  is  shallower  in 
summer  than  in  winter;  therefore,  rates  of  photochemical  trans¬ 
formation  are  predicted  to  be  greater  in  the  summer.  The  blo- 
transfonnatlon  rate  constant  for  TNT  is  an  order  of  magnitude  smaller 
than  the  photochemical  rate  constant;  so  a  comparison  of  the  dis¬ 
solution  flux  and  photochemical  rate  constant  will  indicate  whether 
the  TNT  or  RDX  concentration  will  increase  or  decrease  over  any  par¬ 
ticular  time  period.  Furthermore,  if  the  input  source  is  eliminated 
(as  was  observed  when  portions  of  the  lagoon  dried  out) ,  the  TNT  and 
RDX  will  transform  at  rates  defined  by  the  depth  and  mitigated  by  the 
absorptivity  of  the  water.  Thus,  an  overall  profile  of  munitions  in 
Che  lagoon  will  show  a  buildup  of  munitions  when  the  lagoon  is  full, 
periods  of  equilibrium  when  the  dissolution  and  transformation  fluxes 
are  equal,  and  periods  when  transformation  flux  exceeds  that  of  dis- 
eolutlon--usually  in  the  summer  months. 

The  above  theme  has  a  number  of  variations  that  can  be  postulated 
and  for  which  modeling  can  serve  as  a  tool  to  predict  chemical 
behavior.  For  example,  if  TNT  and  RDX  were  uniformly  distributed  in 
the  bottom  sediment,  as  was  observed  at  the  dump  ramp,  they  would 
reach  their  solubility  limits  in  the  lagoon  water.  Similarly,  if  the 
sediment  had  a  uniform  distribution,  as  was  found  at  the  center  of  the 
lagoon,  the  transformation  of  TNT  and  RDX  could  be  predicted  from  the 
transformation  rate  constants  developed  in  this  study. 
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In  conclusion,  Che  persistence  of  TNT  and  RDX  in  IAAP  lagoons  is 
controlled  by  dissolution  and  photochemical  processes,  with  dis¬ 
solution  being  complicated  by  a  heterogeneous  distribution  of  solid 
munitions  on  the  bottom  surface  and  photolysis  being  complicated  by 
changes  in  absorptivity  of  the  water  that  affects  the  photochemical 
rate  constant  as  a  function  of  time  of  year.  This  effect  could  not  be 
Investigated  in  detail  because  of  project  time  and  money  contraints. 

A  simulation  model  was  developed  to  describe  a  first  approximation  of 
the  loss  and  movement  of  TNT  and  RDX  in  the  lagoon  water,  but  vigorous 
validation  could  not  be  achieved  due  to  sampling  and  mixing  con¬ 
straints  within  the  lagoon. 


Appendix  A 

ANALYSIS  OF  LAAP  LAGOON  NO.  9  SEDIMENT 


Sediment  Preparation 

Twenty-four  sediment  samples  scraped  from  the  bottom  surface  of 
Lagoon  No.  9  were  sieved  through  a  large  mesh  screen  to  remove  twigs, 
leaves,  and  large  pebbles.  Each  sample  was  homogenized  with  a  mortar 
and  pestle  until  an  even  consistency  was  obtained. 

Moisture  Determination 

A  known  weight  (~20  g)  of  sediment  was  placed  In  a  tared  beaker 
and  allowed  to  dry  at  75°C  for  48  hr.  The  percent  moisture  in  the 
sample  was  calculated  as  the  dry  weight  minus  the  wet  weight,  divided 
by  the  wet  weight  times  100. 

Extraction  of  Sediment 

A  sample  of  wet  sediment  (10-20  g)  was  placed  in  a  50-ml 
screw-cap  centrifuge  tube  and  20  ml  of  neutral  ethyl  acetate  was 
added.  The  tube  was  vigorously  shaken  for  2-3  min  and  then  cen¬ 
trifuged  at  2000  rpra  for  10  rain.  The  ethyl  acetate  was  removed  by 
plpet.  This  process  was  repeated  twice.  The  sediment  was  then 
acidified  by  adding  5  ml  of  2N  HC1  and  re-extracted  with  ethyl  acetate 
(3  x  20  ml).  The  combined  extracts  were  dried  over  anhydrous  sodium 
sulfate,  filtered,  and  rotary-evaporated  to  dryness.  The  residue  was 
dissolved  In  methanol,  and  an  internal  standard  was  added  for  HPLC 
analysis . 


The  HPLC  analyses  of  the  sediment  extracts  were  performed  under 
the  following  conditions: 

Instrument:  Spectra-Physice  Model  350QB  Liquid 

Chromatograph 

Column:  250  mm  x  4.6  mm  Altex  RP600-C.o  @  50°C 

Solvent:  A  -  water,  B  -  50%  methanol:  502! 

acetonitrile 

Program:  30%  B  in  A  to  50%  B  in  A  in  15  min;  linear 

gradient 

Flow  rate:  2.0  ml  min-1 

Detection:  UV  @  254  nm 

Retention  times:  5.68  min,  HMX 

6.65  min,  RDX 

10.55  min,  3, 5-dlnitroaniline 
12.40  min,  TNT 

13.77  min,  amlnodlnltrotoluenes 
15.65  min,  3,5-dinltrotoluene  (Internal 
standard) 

A  typical  chromatographic  profile  of  the  sediment  extract  appears 
in  Figure  A-l . 


Results 

The  results  of  duplicate  analyses  of  the  24  sediment  extracts  are 
presented  in  Table  A-l  (algo  see  Figure  2  In  main  text).  Table  A-2 
shows  the  percent  water  in  the  sediment  samples. 


Table  A-l 


PERCENT  COMPOSITION  OF  MUNITION  COMPONENTS 
IN  LAAP  LAGOON  NO.  9  SEDIMENT 


Sample  TNT  RDX  KMX  ADNTSa  3,5-DNA 
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aADNTS  -  aminodinitrotoluenes. 


Table  A-2 


PERCENT  WATER  IN  LAGOON  NO.  9  SEDIMENT  SAMPLES 


Sample  No. 

Percent  Water 

A-5 

51.6 

A-4 

33.6 

A-3 

34.1 

A-2 

43.3 

A-l 

48.5 

B-5 

41.6 

B-4 

32.4 

B-3 

31.0 

B-2 

33.9 

B-l 

42.1 

C-5 

33.3 

C-4 

38.6 

C-3 

30.2 
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32.9 

C-l 

58.2 
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D-l 

43.3 

E-4 

32.2 

E-3 

29.5 

E-2 

30.1 
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45.7 
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Appendix  B 

MASS  TRANSFER  PROBLEMS— ASSUMPTIONS  AND  SOLUTIONS 


Description  of  Mass  Transfer  Problems 

A  body  of  water  containing  some  concentration  of  one  or  more  con¬ 
taminants  is  in  contact  with  a  porous  sediment  (Figure  B-l).  Water 
residing  in  the  pores  of  the  sediment  (interstitial  water)  also  con¬ 
tains  some  concentration  of  one  or  more  contaminants.  The  sediment 
solids  may  contain  the  same  contaminants,  presumed  to  be  at  equi¬ 
librium  with  the  interstitial  water  concentration  of  contaminants. 
Contaminants  will  diffuse  through  the  water  both  in  the  free-water 
phase  and  the  sediment/water  phase.  In  addition,  chemical  reactions 
(photolysis,  biodegradation,  hydrolysis)  may  occur,  and  contaminants 
may  be  transported  by  bulk,  convection  (seepage  of  water  into  the 
sediment) . 


Far  from  the  water/sediment  Interface,  one  expects  the  con¬ 
centrations  of  contaminants  to  reach  a  constant  (known)  value.  Also, 
on  the  free-water  side  of  the  Interface,  a  region  should  exist  where 
the  macroscale  mixing  processes  of  the  pond  are  negligible.  With 
these  ideas  in  mind,  one  can  pose  a  set  of  simultaneous  differential 
equations  for  the  concentration  profile  of  a  contaminant  in  the  water 
and  sediment/water  phases.  ,  . 


The  change  in  concentration  of  chemical  with  time 


a2c 


ac 


sw 


at 


is  equal 


f  flv 

to  the  sum  of  the  diffusion  and  seepage  rates  (D^ — and  v 


oz 


ac 

_ r 

az 


sw 


respectivelyj  minus  the  biotransformation  rate  (kQ)  and  the  amount  of 


/ 


chemical  that  is  sorbed  to  the  sediment: 
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Water 


2  >  0 


Water 


.  -  -  :  \  •  ! ,  .  V;  ^'V'.  \  Sediment/Water ->’ 

-V. ‘  ‘  "  '/• .  V-°  '  V ; 


Seepage 


Diffusion 


Sorption 
'  Desorption 


Transformation 


Dissolution 


Seepage 
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FIGURE  B-l  LAGOON  WATER  RESTS  ON  SEDIMENT 


Pond  and  sediment  are  considered  to  be  semi-infinite. 


Rearrangement  of  Equation  1  yields 


,  SC  d2C  SC 

,  .  „  _n-l  sw  sw  ,  sw 

1  +  KnC  — r—  "  D  - —  +  v  — - k 

sw  St  1.2  Sz  o 

Sz 


where  Cgw  -  chemical  concentration  in  the  se  i  w 

K  -  sediment  partition  coefficient  (C8  ■  Kc£w) 


*  diffusion  coefficient  in  the  sediment /water  phase 

*  seepage  velocity 

*  distance  from  the  sediment 

*  zero-order  rate  constant  for  biotransformation  in 
sediment . 


Similarly,  for  the  free-water  pha^e,  the  change  in  concentration 
of  a  chemical  with  time  can  be  described  by 


SC  s2c  SC 

w  _  w  ,  w  _  _ 

TT —  *  D-  — X—  +  V  - - Ek  C 

St  2-2  Sz  i  w 

Sz 


where  C 


EkiCw 


*  concentration  of  chemical  in  the  free-water  phase 

■  time 

•  diffusion  coefficient  in  the  free-water  phase 

■  seepage  velocity 

“  first-order  rate  expression  for  transformation  process 


The  boundary  conditions  for  these  equations  are  that 
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c 

8W 


(0) 


Cw  <°> 


(5) 


C 

8  W 


(-  ») 


C  ■  known 
sw,® 


(6) 


C  («)  *  C  *  known 
w  w,® 


(7) 


C  (z)  ■  known  at  time  zero  ■  C  ,  (8) 

sw  sw,o  ’ 


C  (z)  ■  known  at  time  zero  ■  C  .  (9) 

w  w,o 

Because  of  Equations  5  and  6,  Equations  2  and  3  must  be  solved 
simultaneously  and  are  not  amenable  to  a  closed-form  analytical  solu¬ 
tion.  The  equation  could  be  solved  on  a  digital  computer,  but  the 
additional  complexity  would  greatly  increase  the  expense  of  running 
the  overall  SRI  pond  simulation  model.  Therefore,  we  have  made 
several  simplifying  assumptions,  described  below. 


Comparison  of  Diffusion,  Convection,  Reaction  Rates 


Because  there  is  a  strong  incentive  to  reach  a  closed-form, 
analytical  solution,  the  first  assumption  is  that  the  concentration  of 
chemical  is  constant  everywhere  in  the  free-water  phase.  This  assum¬ 
ption  allows  Equation  3  to  be  eliminated. 
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Further  simplifications  could  be  made  if  either  the  convective 

(  ac  \ 

terra  (vr—^-jor  the  reaction  term  (k  )  could  be  neglected.  By  non- 

diraensionalizlng  the  equations,  one  can  evaluate  the  relative  impor¬ 
tance  of  diffusion,  reaction,  and  convection.  Unfortunately,  even  the 
single  differential  equation  cannot  be  solved  analytically  if  there  is 
no  boundary  on  the  sediment  phase.  Therefore,  to  compare  the  impor¬ 
tance  of  diffusion,  reaction,  and  convection,  we  pose  a  problem  where 
the  sediment  phase  has  a  boundary  layer  of  finite  thickness,  i,  near 
the  water/gediment  interface  (Figure  B-2).  Furthermore,  we  assume 
steady-state  conditions.  This  assumption  leaves  a  simple  differential 
equation  describing  the  concentration  in  the  sediment-water  phase.  In 
dimensionless  form,  this  equation  is 

d2C*  dC* 

sw  .  „  sw  „  .... 

— *r +  5.7?  -  Yo  0  -  (10> 

dz  dz 

where  ZQ  *  ratio  of  convective  transport  rate  to  diffusive  transport 
rate 

y Q  *  ratio  of  chemical  reaction  rate  to  diffusive  transport  rate 

with  boundary  conditions 

*  C„  * 

Csw  (0)  "  C -  Cw  "  kn°Wn  (U) 

8W,® 

and 


C*  (-1)  -  1 

SW 


The  dimensionless  parameters  are 


* 

C 

sw 


C 

sw 


C 

sw,® 


(12) 


(13) 
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Water 
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FIGURE  B-2  BOUNDARY  LAYER  APPROXIMATION  IN  SEDIMENT/WATER  PHASE 


Beyond  a  depth  of  2,  the  sediment  is  undisturbed  by  the  mass  transfer. 
Hence,  flux  can  be  computed  as  if  sediment/water  phase  extended  only 


The  parameter  La  a  measure  of  the  ratio  of  convective  flux  to 
diffusive  flux,  and  yQ  is  a  measure  of  the  ratio  of  the  rate  of 
reaction  to  the  rate  of  diffusion.  The  flux,  N,  in  dimensional  form 


N  *  '  ♦Dl<fc>z-0  ‘ 


Defining  a  dimensionless  flux, 


"l  ♦dTc  » 

1  sw,« 

* 

we  find,  by  solving  Equation  10  for  Cgw  and  substituting  into  Equation 


17,  that 


5  C  +  ~  e  -  £  -~ 

°  *  h  0  5o 


1  u 
1  -  e 

To  evaluate  the  relative  importance  of  the  various  phenomena, 
Equation  8  was  solved  with  Y0  "  0  (no  chemical  reaction,  N*),  with 
£0  *  0  (no  convection,  N^),  and  with  y0  *  50  “  0  (no  chemical  reaction 
and  no  convection).  The  dimensionless  fluxes  for  each  of  these  cases, 

it  It  it 

N„,  N_ ,  N. ,  are  given  as  follows: 


To  compare  N^  to  Nj,  N^,  and  N^,  we  need  to  pose  some  reasonable 
values  of  the  parameters  C^,  5C,  and  yQ.  SRI  has  taken  data  from 
field  samples  of  the  waste  disposal  pond  and  found  the  following 
numbers : 


V 

m 

2.2  x  10  6  cm  s-1 

(23a) 

ko 

m 

1.2  x  10  ^  yg  cm-'*  s-*- 

(for  TNT) 

(23b) 

D1 

wt 

3.05  x  10"6  cm2  s'1 

(for  TNT) 

(23c) 

C 

sw,® 

m 

50  yg  cm  J 

(for  TNT) 

(23d) 

* 

C 

w 

m 

0.1 

(for  TNT) 

(23e) 

Because  the  boundary  thickness,  l,  exists  only  In  our  model  (not  In 
the  real  pond),  reasonable  estimates  for  Jt  must  be  made.  Therefore, 


In  Table  B-l,  the  results  are  presented  for  various  boundary  thick- 

* 

nesses.  The  table  shows  the  dimensionless  flux  (N  )  and  flux  values 


*  * 

neglecting  chemical  reaction  (N«),  convection  (N_),  and  reaction  and 


convection  (N^),  along  with  the  ratios  ^  /N^  ,  N^/N^  ,  and  N^/N^ 


These  ratios  show  how  much  difference  It  makes  If  the  convection,  the 

* 

reaction,  or  both  are  neglected.  It  can  be  seen  that  N-  Is  not  much 
*  *  1  * 


different  from  N,  and  that  N,  Is  not  much  different  from  N.  . 

13  A 


Hence,  the  chemical  reaction  term  can  be  neglected.  In  contrast,  N. 


,  *  *  *  J 
can  be  much  larger  than  N.  and  N,  can  be  much  larger  than  N-.  These 


^  ^  van  wv  u*wvii  *a  *  O  "  ^ 

results  Imply  that  the  convective  term  cannot  be  neglected.  Hence, 
the  reaction  can  be  neglected  but  convection  must  be  considered.  Thi 
conclusion  was  expected  from  the  magnitude  of  the  dimensionless 
numbers  ( £0  -  1  and  yQ  <  1  (Table  B-l)]. 
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Unsteady  Solution  with  Convection,  Without  Reaction 


By  neglecting  the  chemical  reaction  term,  we  can  develop  a 
closed-form  solution  for  flux  without  appeal  to  some  arbitrary  boun¬ 
dary  layer  thickness,  l.  In  this  regard,  the  solution  developed  below 
Is  the  most  realistic  situation  for  which  a  closed-form  solution  can 
be  developed.  The  solution  will  apply  to  TNT  and  RDX  because  the  RDX 
was  found  not  to  react  In  the  sediment/water  phase  (l.e.,  kQ  -  0  for 
RDX). 

Physically,  we  envision  a  situation  like  that  In  Figure  B-l.  We 
retal:;  the  assumption  that  the  concentration  In  the  water  phase  Is 
constant  and  known.  For  this  situation,  the  following  equations  apply 
In  the  sedlroent/water  phase: 


.  dC  d2C 

.  .  ..  „n-l  sw  _  sw  . 

1  +  Kn  C  ■  ...  -  D.  - =— +  v 

sw  ot  1.2 

bz 

ac 

sw 

az  » 

(24) 

with  the  following  boundary  conditions: 

C  (-  ®,t)  -  C  -  known  , 

8W  8W,® 

(25a) 

C  (0,t)  ■  C  ■  known 

8W  ’  w 

» 

(25b) 

C  (z,0)  -  C  •  known 

sw  sw,® 

• 

(25c) 

With  n  ■  1  (linear  Isotherm  assumption)  and  In 

these  equations  become 

dimensionless 

form, 

*  2  *  * 

ac  a  c  ac 

SW  8W  n  sw 

*  *2  Po  * 

dt  dz  dz 

> 

(26) 

where  (30  "  dimensionless  ratio  of  convective  transport  rate 

to 

diffusive  transport  rate. 

and 


A  A 

C  (-  ®,t  )  **  0 

SW 

(27a) 

A  A  A 

C  (0,t  )  -  c  -  1 

SW  *  w 

(27b) 

A  A 

C8W  U  ’0)  "  °  * 

(27c) 

where  the  dimensionless  variables  have  been  defined  as  follows: 


* 

C 

sw 


C  -  C 

8W  SW,” 

C 

8W(® 


* 

fc 


l  <14*> 


(28a) 


(28b) 


* 


z 


*  z/z 

o 


(28c) 


P0  *  <28d> 

The  quantity  zQ  is  an  arbitrary  length  scale  that,  if  reasonably 
chosen,  can  give  some  idea  of  the  relative  importance  of  diffusion  and 
convection.  (The  final  answer  will  not,  of  course,  depend  on  zQ.) 

Equation  26,  subject  to  the  boundary  conditions  of  Equations  27a- 
27c,  can  be  solved  by  using  LaPlace  transformations.  One  finds 


*  *  * 

C  (z  ,t  )  - 
sw 


erfc(ri1)  +  exp  |-Pqz  jerfc(n2)| 


(29a) 
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where 


and 


(29b) 


n2 


(29c) 


The  flux,  N,  la  given  by 


N 


z_0 


-  $v  C  (o) 
sw 


In  dimensionless  form. 


* 

N 


z  N 
o 


♦Dj 


C 

BW,® 


* 

-pc 
Ko  w 


(30) 


(31) 


Differentiating  Equation  28  and  substituting  Into  Equation  31,  one 
f  lnds 


* 

N 


*2  1/2 

<l  -  C  )  -  Bt*  P  *  P  e*  ' 

— 172  — >  ’  r<1  -  Cw)orfc(f")  - 


(32) 


Equation  32  can  be  cast  In  a  fore  that  explicitly  removes  zQ  frorc  the 
equation  for  the  flux.  If  we  define  a  new  dimensionless  flux. 


N 


N 


$vC 


SW,“» 


(33) 


and  a  new  dimensionless  time, 


then  Equation  32  becomes 


*  *  1/2 

(1  -  C  )  _  1-C  ~  * 

N  -  - exp(-  t/4)  -  (-y-*)  erfc( — j~)  -  Cw 


Relationship  To  SRI  Lagoon  Model 

Equation  32  is  the  basis  for  computing  how  much  material 
transfers  into  or  out  of  the  sediment  during  a  given  time  interval/ 
If  M  is  the  total  mass  of  material  per  unit  area  transported  across 
the  water/sediment  Interface  in  a  time  Interval,  x, 


M  -  /  Ndt 


$D.  <14K)C 
i  sw,® 


Note  that  the  z  direction  was  chosen  positive  in  Figure  B-l. 
Therefore  N  >  0  implies  that  material  moves  from  the  sediment  into 
the  water. 
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and 


2 

V  T 


D  L ( 1+K) 


(38) 


Substituting  Equation  35  into  Equation  36,  one  finds 


M  "  (l"Cw)  \  erf  \ — 2 — /  '  2“^rfc  V — 2 — /  +  ■7^pexP("'c*/4)/  ~  C*T*  . 


(39) 

The  functions  erf  and  erfc  are  well-known,  tabulated  integrals.* 

These  functions  are  available  on  any  computer  that  has  the  standard 
IMSL  math  library  software.^ 

In  the  SRI  lagoon  model,  one  knows  C*  and  CflW  „  at  the  beginning 
of  a  time  step  (and,  by  implication,  at  time  zero).  During  a  time 
step,  one  allows  mass  to  transfer  between  the  water  and  the  sedi¬ 
ment/water  phases.  The  total  amount  of  mass  transferred  is  given  by 
Equation  39.  At  the  end  of  the  time  step,  the  new  values  of  C*  and 
CflW,*  are  computed  based  on  the  addition  or  depletion  of  contaminant 
mass  in  the  water  or  sediment/water  phases.  Computations  are  repeated 

as  desired  with  these  new  values  of  C*  and  C  ,  . 

w  aw1® 


Carslaw,  H.S.  and  J.C.  Jaeger  (1959)  Conduction  of  Heat  in  Solids. 
Oxford-at-the-Clarendon,  London. 

'Abramowitz,  K.  and  I. A.  Stegun  (1964)  Handbook  of  Mathematical 
Functions,  National  Bureau  of  Standards,  Applied  Mathematics  Series 
55,  June,  1964. 
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Sample  Calculation 

The  flux,  N,  and  the  total  amount  transformed,  M,  will  be  com¬ 
puted  for  the  particular  parameters  relevant  to  the  TNT  and  RDX  dis¬ 
posal  lagoon.  The  values  of  the  parameters  are  given  by  Equations 
23a,  23b,  23c,  and  23d.  In  addition,  we  must  specify  a  value  for  K 
and  $.  For  TNT,  K  —  10;  no  data  are  available  for  RDX.  Consequently, 
the  calculations  below  apply  only  to  TNT.  The  porosity,  $,  which  has 
been  measured  by  SRI,  is  0.2. 

With  those  values  of  the  parameters.  Tables  B-2  and  B-3  can  be 
generated  from  Equations  33  and  39.  It  Is  worthwhile  to  record  the 
values  of  the  groups  of  parameters  appearing  in  Equations  33,  34,  and 
37.  Each  group  is  a  "characteristic"  value,  and  we  define  the 
following  parameters: 


-2  -1 

N  *  4>vC  ■  1.872  ug  cm  day 
o  sw,®  ro  3 


(40) 


M 

o 


(1+K)  C 


8W,® 


V 


154.8  pg  cm 


-2 


(41) 


11; 


Table  B-2 


I 


I 

I  FLUX  CALCULATIONS  FOR  SAMPLE  LAGOON  PROBLEM 


Time 


Flux 


t 

t  (hr) 

N 

N  (uk  cm”2  day  *) 

0.01 

19.8 

4.57 

8.55 

0.04 

79.3 

2.01 

3.77 

0.09 

178 

1.18 

2.21 

0.16 

317 

0.769 

1.44 

0.25 

496 

0.528 

0.989 

0.36 

714 

0.371 

0.695 

0.49 

971 

0.263 

0.491 

0.64 

1269 

0.184 

0.344 

1.00 

1982 

0.080 

0.149 

Table  B-3 

TOTAL  MASS  PER  UNIT  AREA  TRANSFERRED  DURING 

A  TIME  INTERVAL  FOR  SAMPLE  LAGOON  PROBLEM 

Time 

Interval 

Mass 

Transferred 

* 

T _ _ 

T  .(Hr) 

M 

M  (ng  cm-2) 

0.01 

19.8 

0.0961 

14.88 

0.04 

79.3 

0.182 

28.14 

0.09 

178 

0.257 

39.86 

0.16 

317 

0.324 

50.10 

0.25 

496 

0.381 

58.95 

0.36 

714 

0.429 

66.48 

0.49 

971 

0.470 

72.78 

0.64 

1269 

0.503 

77.89 

1.00 

1982 

0.548 

84.86 
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1982  hr 


(42) 


D^l+K) 


Equation  42  shows  that  a  steady-state  solution  for  the  flux  would 
be  Inadequate  unless  time  scales  of  Interest  were  much  longer  than 
1982  hr  (~  83  days).  Time  steps  taken  In  the  SRI  lagoon  model  are 
variable;  however,  time  steps  shorter  than  83  days  necessitate  the  use 
of  an  unsteady  solution  as  developed  herein.  Equation  40  states  that 
a  characteristic  value  of  the  flux  is  approximately  2  pg  cm  2  day  * . 

Part  of  the  significance  of  NQ,  MQ,  and  tQ  Is  that  any  situation 
with  the  same  values  of  NQ,  tQ,  and  Cw  will  be  the  same.  MQ  is  not  an 
Independent  characteristic  value  because 


M  -  N  t 
o  o  o 


(43) 


Addition  of  Dissolution  to  the  Mass  Transfer  Problem 

The  above  solutions  and  assumptions  were  derived  for  the  case 

where  the  dissolution  rate  constant  (kr)  equals  zero.  We  must  also 

account  for  the  case  where  kr  *  0,  which  occurs  when  the  lagoon  Is 

full  or  when  the  lagoon  waters  cover  the  area  near  the  dump  ramp.  As 

previously  developed,  the  dissolution  rate  of  solid  TNT  or 
dw 

RDX,  *7-,  can  be  written  as 
’  dt* 


dw 

dt 


k  (C  -  C.  ) 
r  sat  1 


(44) 


where  kr  Is  an  experimentally  determined  dissolution  rate  constant, 
Cgat  Is  the  saturation  concentration  of  the  dissolving  species  in  the 
liquid  phase,  and  Is  the  concentration  of  the  dissolved  species  at 
the  solld/liquld  Interface.  If  diffusion  from  the  Interface  to  the 
surrounding  liquid  is  fast  compared  to  dissolution,  then  Equation  44 
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* 


can  be  represented  as  shown  in  Equation  45: 


-  k  (C  -  C  ) 
dt  r  sat  sw 


where  Cgw  is  the  concentration  of  TNT  or  RDX  in  the  sediment/water 
phase  at  any  depth  in  the  sediment.  Therefore,  Equation  45  can  be 
incorporated  into  the  transport  equation  developed  in  the  previous 
section,  as  shown  in  Equation  46: 


bC  0  C 

<l  +  K)"sr '  “i~rf +  '  c»"’  * 

Oz 


where 


K  -  (^-)K  F  , 

$  eq  oc 


and  where  FQC  is  the  organic  fraction  of  the  sediment  solids,  S^gt  is 
the  surface  area  of  solid  per  unit  mass  of  solid  (typically  determined 
by  the  BET  absorption  method),  pfcs  is  the  mass  density  of  the  sediment 
solids,  4>  is  the  void  fraction  and  F  is  the  fraction  of  the  solid  sur¬ 
face  area  that  is  covered  by  the  dissolving  species. 

A  linear  isotherm  has  been  assumed  in  deriving  Equation  46,  with 


C  -  K  C 
s  eq  sw 


where  Cg  is  the  concentration  of  TNT  or  RDX  dissolved  in  the  organic 
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fraction  of  the  sediment  solids  and  KfiCj  is  experimentally  determined. 
The  boundary  conditions  for  Equation  46  are: 


C  (_®,t)  ■  known  a  C 

sw  sw 

t  >  0 

i® 

(50a) 

C  (0,t)  -  C  t  > 
sw  w 

0 

(50b) 

c  (z,0)  -  C  . 

sw  sat 

» 

(50c) 

where  C  is  the  concentration  of  TNT  or  RDX 
w 

in  the  bulk-water 

phase 

The  flux,  N,  is  given  by 

foe  \ 

N  ■  -♦d4-» r)z  -  0  +  (l-*>krF(C,at  '  C«)-  *vc„  • 

(51) 

If  Equation  46  can  be  solved  subject  to  Equations  50a-50c,  then  the 
flux  can  be  computed  from  Equation  51.  Equation  46  will  describe  the 
flux  in  the  laboratory  experiment  and  in  the  lagoon.  With  laboratory 
data,  we  can  find  the  value  of  kr  and  use  it  to  predict  the  flux  in 
the  lagoon. 

Solution  to  Flux  Equations  for  Laboratory  Experiment 

In  the  previous  sections,  Equation  46  was  solved  with  kf  ■  0 

(i.e.,  no  dissolution).  In  the  laboratory  experiment,  however,  v  ■  0 

and  kr  *  0  (i.e.,  there  is  no  seepage  but  there  is  dissolution  of 

solid  TNT  and  RDX).  Also  in  the  laboratory,  C  ■  C  . 

sw,®  sat 

Consequently,  the  mass  transfer  equations  become 
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(1  +  K) 


ac 


sw 


d2C 


sw 


at 


dz 


+  WCsac  ‘  CS»’ 


withC  (-*,t)  ■  C 

sw  sat 


C  (0,t)  -  C 
sw  w 


and  C  (z,0)  ■  C 

sw  sat 


The  flux  Is  given  by 


At  steady  state,  we  can  solve  Equation  52  subject  to  the 
conditions  53a  and  53b.  We  find 


C  -  C 
sat  sw 

C  -  C 
sat  w 


»  exp 


f/k 

^  D,  / 


1/2 


The  flux,  N,  Is  then  given  by 


r  .  <1  +  pD)  (1  -  cu)  , 


"lth  "o  '  Csat  <*DlkrrSbet',be‘1-»»1/2  • 
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(52) 

(53a) 

(53b) 

(54) 

(55) 

boundary 

(56) 

(57) 


(58) 


(59) 


.  y->> 

D 


*  w 

and  C  -  7 •  (60) 

w  C  _  v  7 

sac 

In  Equations  58,  59,  and  60,  $,  D,  and  C8ac  are  known.  The  surface 

fraction,  F,  approximately  equals  the  (known)  mass  fraction.  The 

product  p.  can  be  estimated  by  assuming  a  reasonable  particle 
bet  bs 

size  for  the  sediment  solids.  If  the  particle  diameter  is  d^,  then 


Sbetpb8 


6_ 

d 


(61) 


Sediment  that  can  be  described  as  fine  soil  or  clay  has  particle  sizes 
in  the  range  of  10  to  100  pm.  (The  square  root  dependence  of  the  flux 


on  Sj)etP^ja  de-emphasizes  the  Importance  of  a  proper  choice  for  dp.) 


Because  all  parameters  in  Equations  58,  59,  3nd  60  except  kr  are 
known  or  can  be  reasonably  estimated,  the  laboratory  flux  data  can  be 
used  to  find  a  best  value  of  kr.  To  make  this  procedure  a  bit  more 
transparent,  we  rearrange  Equations  57  through  62,  substitute  Equation 


61  for  S.  p,  ,  and  obtain 
bet  bs 


N  -  (C  -  C 
sat  w 


1/2 


+  krF(l-<fr)> 


(62) 


In  the  flux  experiments,  N  and  Cw  are  measured  and  a  plot  of  N  versus 
CBat-Cw  will  have  a  slope  equal  to  the  bracketed  expression  in  Equa¬ 
tion  17. 

A  least-squares  fit  of  the  data  to  the  form  y  ■  mx  yields  the 
following  values  for  the  bracketed  expression  in  Equation  62: 


~-y 
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& 


,N  / 
A  ' 


ta 


£6 

V-*.' 

I 

>>>} 

*  ■>' 
AYS 

v-v- 


TNT 


L(i-O)  (^-)  kr: 


1/2 

F |  +  krF(l-^) 


-4  -1 

1.0776  x  10  cm  sec 


RDX 


♦(!-♦) 


(?)  V 


1/2 


+  krF(l~$)  *  7.752  x  105  cm  sec1 


(63) 


(64) 


The  values  of  the  known  parameters  In  Equations  63  and  64  are 

listed  In  Table  B-4  for  TNT  and  RDX.  The  diffusion  coefficient  is  the 

* 

average  result  from  Equations  8-34  and  8-36  of  Reid  and  Sherwood. 

The  void  volume  fraction,  $,  was  calculated  from  the  dry  weight  of 
sediment  solids,  measured  prior  to  the  flux  experiments.  The  sat¬ 
uration  concentrations  are  known  from  solubility  data  reported  in 
Section  V.A.,  Physical  Transport.  The  surface  fraction,  F,  is  assumed 
equal  to  the  mass  fraction  of  TNT  or  RDX  in  the  sediment  solids,  a 
quantity  measured  prior  to  the  flux  experiments.  The  value  of  dp  is 
an  estimate. 

From  Equations  63  and  64  and  the  values  of  the  parameters  in 
Table  B-4,  we  find  that  the  kr  is  1.172  x  10“5  cm  sec  for  TNT  and 
1.164  x  10~5  cm  sec  *  for  RDX. 


*Reid,  R.D.  and  T.K.  Sherwood  (1958)  The  Properties  of  Gases  and 
Liquids.  McGraw-Hill,  New  York. 
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The  flux  derived  from  Equation  62  with  these  values  of  kf  has 
been  placed  on  Figure  B-3  with  the  data  points  to  shov>  the  goodness  of 
fit.  One  can  see  that  Equation  62  Is  within  a  factor  of  2  in  pre¬ 
dicting  the  flux  over  the  range  studied  experimentally.  This  degree 
of  accuracy  can  be  disturbing  unless  one  allows  for  the  inherent 
difficulty  of  performing  the  flux  measurements  and  of  describing 
analytically  the  mass  transfer  in  a  multicomponent,  heterogeneous 
system.  One  must  also  refrain  from  the  temptation  to  fit  the  data  in 
Figure  B-3  to  a  straight  line  of  the  form  y  ■  mx  +  b.  The  result  of 
this  seemingly  intuitive  approach  is  a  model  that  predicts  negative 
flux  values,  a  physically  impossible  situation  in  the  laboratory 
experiment. 


Table  B-4 


PARAMETERS  FOR  ANALYZING  FLUX  DATA 


Parameter 


(25®C)  2.13  x  io“5  cm2  sec"1  2.24  x  10' 


C8at  (25°C) 


131  tig  cm" 


50.6  jig  cm" 


0.156 


0.078 


50  )im 


50  tun 


S  ^  aJIll 


TNT  FLUX  109-N  (gem 


RDX  12345678 
CONCENTRATION  IN  BULK  WATER  (pg/ml) 

LA-7934-27 

FIGURE  B-3  TNT  (•)  AND  RDX  (O)  FLUX  DATA  COMPARED  WITH  CURVE 
FIT  (EQUATION  62) 


RDX  FLUX  10y  *  N  (g  cm 


Appendix  C 


SRI  COMPUTER  PROGRAM  FOR  MODELING  THE  FATE  OF  TNT  AND  RDX 
IN  LAGOONS  CONTAINING  MUNITIONS  WASTES 


This  appendix  presents  the  computer  program  developed  to  simulate 
TNT  and  RDX  in  a  waste-lagoon  environment.  It  encompasses  a  number  of 
subroutines  to  determine  changes  in  volume,  surface  area,  and  depth  of 
the  lagoon  and  the  transport  of  munitions  (flux)  from  the  bottom  sedi¬ 
ment  to  the  free-water  phase.  The  subroutine  COMPS  performs  the  cal¬ 
culation  of  munition  concentration  in  the  bulk-water  phase  and  in  the 
suspended  sediment.  The  data  output  can  be  conveniently  displayed 
graphically  by  the  subroutine  GRAPH  4. 
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10501  0019 

10502  0020 

-10503 _ 0021 

10504  0022 

10505  0023 


10507  0025 

105u8  0026 


30600 

30610 

JO/Ofi... 

10/10 

10720 


30900  0034 


go  to  500 
16*06291. 
go  to  500 
temp«303, 
go  to  500 


l  M1—  mUI'II 


boo  ta=sir'»10.**(8.)»exp(-ts/temp) 
oo  100  1=1 ,n) 

let.  I  (  NS.l  )- 


oml  J )  =  rxec«t ( 3 )*area/l 000. *  c  ta*looo ,»c) 
c  om (33  is  in  ug  per  day 

_L0.a__.UMA55=DHAS6ttJ«UJ. _ 

•kITt(6,2O0)DMAS5 
200  POKMAT l 7rt  DHASSc ,  1 1  2 . 4  ) 

KN  _  _ 


thO 


<< 


1  75uo 

2  fttUU 

-  .2  77  00 . 

2  Jluu 
jJaio 

o  jam, 

OJD  1 

OoO  2 

0004 

OOOS 

-X _ _ 

c 

SUHMjUllfck.  tLUA 

1  n±i>  .SUfchuUllhL  IS  IU  CAlXULAIt  JHi:  IvlAl  llAtoBlilU  _U^ 
ChtrtJCAL  ImmouGH  SEDlMf.Nl  lhltklACE 
co*"'on/p6i  «/n»  ,  n«  ,at  ,tl«ie,prt,t,»t 

rriMank.  /chi  si  /o.i  r  Ai  Bui  .  m«  .  or  .  n  1  k  .  PH  1  . SI) .  sf  IJliR  .  Mil. 

- - 

2buuq 

0007 

COMMUN/CUNO/C,  TSS.SSP  ,CS»(50)  ,CkT  150) 

JBUiU 

0000 

COMMON /PHI  S3 /D*< SO  >,UMA5S,ThaS(S0 J 

COP<tnfc/pH»s'2/pi<,£,v,5P,vOi.(PD3PlHi6U)(.4rea.Qeptri(50l 

26055 

0010 

CO»»on/»6ler/n(50 ) 

26057 

001 1 

co*K>on/*olu/t6 

001  ? 

2«0bu 

0013 

1000 

lotnuiri  tlux.lS) 

20IOU 

0014 

THASS*0. 

001 

fin  loo  1*1 r n J 

26300 

0016 

3«h((Ns«n-n 

26310 

0017 

c 

2.T  \2u 

0018 

CM! 

t«mperftt.ur*  id^ultiitnt  on  solubility 

J0J  JO 

0019 

C 

26340 

0020 

te*t«/e»*>(  j )« looo. 

28  JSll 

0021 

iiiaml  .It.  ta»i000.)  tpsi. _ 

26400 

CO^OlCC/ lc**( 

2BoU0 

0023 

ID*SP»»2,»07/ ID  If »(1,«AUPHA) ) 

ffi;K^TT7ITl 

'nTTTH 

aael .-CGNDL 

2  B  BOO 

0025 

t>b»exp(*ia/4.0) 

28900 

0026 

a«l.*eriua»M0.5)/2.) 

0027 

tnail  l  i«i«a»bD/l3. 1416»td)»»( 0.5) -6t/2. 0»a-conal j 

29010 

0026 

c 

29020 

0029 

C*** 

eeta  in  nicrograa  per  iquare  cenlieeier  per  aay 

2901Q 

0030 

c 

29100 

0031 

DCtexpnl •tp»c*6())*tb»10,»4(*J) 

291  06 

0032 

29100 

0013 

TMASS»TMASS41MA6(J)4*BtA 

29400 

0034 

100 

CONTINUE 

29600 

0035 

RETURN 

WUN  o 


2o«lo  uobd 
26820  0059 

-2691UI - OOell - 

26910  Oool 
269J2  U062 

-26SLU - Ou6j - 

26914  0064 

20916  0065 

-26916 - 0066 _ 

26920  0067 

26930  0068 

■2.6950 _ 01169 _ 

27000  0070 

27100  0071  C 

-mill _ 01)12 - 

27200  0073 

27400  0074 


t»tma*«/v6lsco*( t  opt -roev«ro*P )/wol )/ 1 1 .•i***pc*10.*« l-J.l ) 

lf(C  .Qt.  0.)  oo  to  500 

_ C6C0  4  014.1  •S  UBdCtC - — - 

$it««6*/vol«c»(r6pr«ro*w-ro*p)/vol)/il .«t»s»pc*io.*» i«3. ) ) 

500  11 (C  .Cl.  0.)  CO  TO  700 

- 1>  CUl  -.KU.-l.)  £U_ XQ— 7.00 - - - 

DT»DT/7. 

Pt'T»RHT/7. 

_ CU-XCJ— 6011 - - — --  ...  - 

700  CONTINUE 

Ulc  ,qt.  ta*1000.)  catMlOOO. 

_ me.  JU-DJ  c»o. _ 

6SP«C/100u.*PC 
SSP  IN  pg/c 

_ tmMeJbui _ _ — - - - 

RETURN 

END 


PROGRAM  SECTIONS 


1 tABt 

Rvtee 

Attribute! 

(I  SCOPE 

688 

PIC  CUN 

REL 

LCL 

SHK 

EXE 

HP 

NUhRT 

LONG 

1  SPUATA 

158 

PIC  CON 

HF.k 

LCL 

ShR 

HUEXE 

HP 

NQmRT 

I’TlM'BHNkBRIII 

2  SLuCAL 

52 

Pic  Con 

REP 

LCL 

NUSHP 

NOEXE 

RD 

hRT 

long 

3  PARA 

28 

PIC  OVR 

REL 

CBL 

SHR 

hUEXE 

HD 

HRT 

LONG 

*  4  RATE 

20 

PIC  OVh 

KF.ii 

Cbb 

SHR 

HUEXE 

PD 

hHI 

LONG _ 

5  PHISl 

36 

PIC  OVR 

REL 

CBL 

SHh 

NOEXE 

HD 

HRT 

LUNG 

6  PHIS2 

420 

PIC  OVH 

HEP 

CBL 

SHR 

HUEXE 

RD 

HRT 

LUnG 

PF.li 

QBk 

SHH 

HD 

HRT 

8  CONCN 

412 

PIC  OVN 

REL 

CBL 

SHR 

HUEXE 

RD 

HRT 

LUNG 

9  MASS 

4 

PIC  OVR 

KEL 

CBL 

SHR 

nuexe 

RD 

HRT 

LONG 

_ 10  RAI61 _ r 

8 

PIC  OVR 

IWl 

BT1 

rpiTr*i 

RD 

Ktl 

|  q  i]  1 1 1 

11  CONST 

8 

PIC  OVR 

REL 

CBL 

SHR 

HUEXE 

RD 

»RT 

LUNG 

12  *AltR 

200 

PIC  OVh 

HEL 

CBL 

SHR 

HOEXE 

RD 

HRT 

LUNG 

_ LI _ nEAlrl _ 

12 

PIC  OVR 

REL 

CBL 

SHR 

HUEXE 

RD 

HRT 

LUNG 

14  SOPU 

4 

PIC  OVR 

REL 

GBL 

SHR 

HUEXE 

RP 

HRT 

LONG 

15  COhlRL 

12 

PIC  OVH 

REL 

GbL 

SHR 

HUEXE 

RP 

HRT 

bON(j 

Total  space  Allocated  2466 


EM  HI  PU1NTS 

_ A  d  ares* _ type _ n  aE«_ 

o-oooooooo  COMPS 
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23300  OJ  u2 
2b400  uoOJ 
2bbo0  0004 
2bouo  ou05 

-  -2b  740 _ OUUe__ 

2bhoo  0Uo7 
2b9oo  ooob 
_  _2ouuo._.J)i).Q9.. 
2o010  0010 

2b  1  Oy  ooll 

2ol 10  0011 

2o  1 20  0014 

- 2o3-Ui - OOJJi _ 

26132  OOle 
2ell4  0017 

- 263-16 _ 001  a.,. 

2ol40  0019 

26160  0020 

- 2nJ_Ui _ 002 1. 

26100  0022 
26102  0023 

— 261314 _ 0.024 

2olbe  002b 
2616b  oo26 

— 263jL9 _ 0027 

26190  0028 

26192  0029 

- 26200  0030 

26400  0031 

26410  0032 

— 2nA2n _ 0033 

26430  0034 

20440  003b 

—26.4.45 _ VO  36 _ 

2o4b0  0037 

2o45b  0038 

_2e4bJ _ 0039 _ 

2o4o0  0040 

26462  0041 

-26.46b _ 004  2 

20407  Ou43 
2040b  0044 

-2  047(1  0045 

2b47b  0046 

2b4b0  0047 

2648b _ 0048 _ 

2o4b7  0049 

2b440  0060 

.2  6  4  9_2 _ 0.0  bj _ 

2o49b  0052 

2bb00  00b  3 

2bboe _ 00b4 _ 

2o7  00  oubb 
2b l 10  0  u  b  6 

26800 _ 00b  7 _ 


3bO 

C 

C  ♦_»  ♦ _ 

C 

oOO 


SUbhuL'lTNi  CO"PS  .  ...  _ _ _ 

This  SUBROUTINE  IS  10  CALCULATE  THE  CONCENT H*1 luV  IN  FREE  WATER 
A  NO  0I>  SUSPEuOEt’  SOLID 

- COHKL>U/PAhk/l,S.,N**UUllML<PRlsUia _ _ _ 

CUIKOn/KAl t/*KU ,»RF  »«KH|6K*,6Kb 

CUi'inON/PHTSl/Nj,  ALPHA  ,lHASS,RC.DIt  , Pnl ,  Mj  , SCOUR ,SETL 

- COl*l,'U!|i/Pii>S2/PK«L<  »SP  j  VOLj  POLP.lHj.bO  1 .4te*.aept>itb()) _ 

COMMON /PHY S3/DR l bC  ) #0RA5S » 1 R ASlbU ) 

COMMON /CUNCh/C,lSS,SSP,CS»t bo  1 ,CMT(bo) 

- co«'B«ii/inai*/3B. _ _ _ 

CClMNUN/ RATE1/A6KP  ,  A»Kb 
COMRON/CONST/A ,  B 

- COAf'.01'-/*Ali;K/>  LbiU _ _ _ _ 

coflinoii/bbatn/iopr  ,ro*v,ro*p 
eoniAion/solu/ta 

- cafflffloiucontxlyKA-,  xJuxn _ _ _ 

av6*voi«io*»(3. )/iarea»n}l 
wr  1 1 e ( o , 400 )  avh,a«Kp 

8*10,4)  **'  af"  ***  -,*i0*‘t'-2Bh  dO}usleo  pootolvs.is  rates^. 

il(a»od(t .be.)  ,n«,  o.)  go  to  900 

- U  C  L.  ,.gt^_2  2 ,4*J_a.»  -3 _ _ _ 

p«B*a 

a«Kbsfckb*2,M  it) 

- Ltix3_ea._U  a*xb«4xt) _ _ _ 

)  «ru*(b,800)  awXb 

i  torn.at  (bh  abk.o*,*12,4) 

- &UBX»JbX04_a*Kp.4M>31i)iXiiAXXB _ _ _ 

cO*c  ______ 

- c.al3_t  oi_tl  u  x3_ _ _ _ _ _ 

call  tluxl 

call  tor  UuT  “  — “ — " — 

i _ call  _.tl  ox _ _ _ 

do  3oo  1*1 ,nj  ~  ~  - - - 

Jsh{ (NS4 l )•!) 

_ tI*3MASLUl*AREA4DHlJ) 

lt(kk  .eg.  2)  ttsamljT  ~  - — 

It  UMUD(I,7.)  ,NE.  0.)  00  TO  300 
_ *lUltibj2  OJ.)  _  J  ,  TT  _ _ 

FURMAllZOh  U It  FUSION  RATE  FKUM,lbf4H  IS  ,E12.4)  - 

CONTINUE 

_ IHASS®3M.A8StDXAS5 

ltock  ,eq,  2)  titass*amats 
wRlTtl6,3bO)TMASS 

FORK AT (  4oh  TOTAL  DI FF  USlOfTTATE  THROUGH  BOTTOM  SEDIMENT* ,E1 2  , 4 ) 

_ S5_l N_C / LITER, C  IN  UG/L1TER _ 

continue 

_ c=cO  +  (dt /2  )  *  (  »sutpi?*cO 


ini**  003H  1 1  (Kii  ,ne,  V)  gc  to  000 

Ulii  o ob*  ItCt  .tit.  16R.)  »NRt*2, 4*k,pt 

—  ---44  23*-..  -OuoO -  .Sou  ...  continue . . . . . 

W4/<*0  0061  lt(»Kpt  ,9t.  N*P46.)  whptiVkpeV. 

2**45  0062  d«kp>a**p**kpt 


i - U3±a~ 

_ UObl 

- LHAHUHLl^.1^1.  _HE_  0.1  C(1  IL  -J 00 

24000 

0064 

*rlt*lb,ibOj  nUns*l)-L),poepthU),aeua  . 

24010 

0065 

200 

continue 

- - 24alx- 

_0£l6b _ 

—  iiXnJ.— eu_  01  .ao  to  350 

24020 

0067 

AnKPsAnKP/NJ 

24400 

0068 

150 

tor»«t(15,2el2.4) 

- 1U9SL2— 

—111169 _ 

.350 

- continue 

24410 

0070 

*riteib,looojnj 

24920 

0071 

1000 

tormat(8n  surface. lb) 

_ 25uuO 

_ 00X2 _ 

-  — Hi_XllRU-. 

2S100 

00  73 

END  ~~  ~  "  ~~ 

HKuGhAh  stTTin«f 

Name 

Bytes 

Attributes 

0  scout 

1  spuata 

- 2-SUtCAL _ 

645 

31 

56 

PIC  CUN  REL  LCL 
PIC  cun  REL  LCL 

SHR  LXE 

SHR  NUEXE 
NtlSHR  NUFX1 

HO 

RD 

RD 

NUNPT 

NUnRT 

LONG 

LUNG 

i  PARA 

4  PhtSl 

- 5_ms3 _ 

28 

36 

_ 420 

PIC  OVP  REL  GbL~ 
PIC  UVK  REL  CBL 
PIC  CIV R  REL  CRL 

SHR  NOEXE 
SHR  NOEXE 
shr  hurts 

RD 

RD 

Rll 

NRT 

NRT 

LONG 

LONG 

6  NAlbK 

7  HAltl 

,  B  rnuTKI, 

200 

B 

12 

PIC  UVP  RtL  GbL 
PIC  OVR  RtU  GbL 
PIC  CiVR  REL  GbL 

SHR  NUbXE 
SHR  NOEXE 
shr  hnrxj 

RD 

RD 

Oil 

- Hri . 

WHl 

NRT 

WUT 

LONG 

LOj*G 

9  R  Alb 

20 

PIC  OVR  REL  GbL 

SHR  NOEXE 

PD 

WHT 

LONG 

— - lfltal  Space  Allocated _ 

- I4.6fe.__ 

-taxtxjfculius _ 

Adore**  Type  Name 
U-OOUUOUOO  SURFACE 


VARIABLES 


-  Aflg.rjL»JL_ 

-IYB«_ 

Name 

Addrex* 

TYD» 

4-0000O004 
— - - 2-uoodouj_i_ 

R»4 

R»4 

ALPHA 

DELIA 

5-OOUOOODB 

2*000000 1 R 

KM 

p*  4 

AKLA 

5*00000004 

R«4 

EV 

2-o'ooo'ooob 

1*4 

1 

b-OOOUUOOO 

1*4 

KK 

8*00000004 

1*4 

KL 

2-UUuil0i)l£ 

1*4 

M 

4-00000000 

1*4 

NJ 

4-OOOuUOOC 

R*4 

PC 

4-00000014 

P*4  " 

PHI 

4-OUOUOOlC 

R*4 

SCOUR 

4-00000018 

P»4 

SD 

.2-00000000 

R*4 

SUM 

2-00000010 

R*4 

SUM  1 

Address 

Type 

Name 

Address 

jype 

Name 

7-00000004 
_ _ 4-(KiU0UO10 

R»4 

R  *4 

A*Kb 

DIF 

7 -OOOOOOOO 
J-OOOOUOUB 

H*4 

H*4 

AN  Up 
D1 

2-00000004 

J  *4 

J 

2-OCVddOOC 

1*4 

K  " 

B -OOOOOOOB 

1*4 

KN 

2-OoOoOo24 

144 

L 

3-00000000 

1*4 

NS 

3-00000004 

1*4 

NN 

b-OOOuOOUO 

R»4 

PP 

j-oooudofo 

_  R*  4 

PKT 

4-UOOOOU2U 

R  *  4 

SETL 

b-UOOOOOOB 

R*4 

SP 

i-OoUOoO 1 4 

P»4 

1 

3-OOOOOOOC 

P  •  4 

TIME 

136 


218WU  00  01 

21»uo  0002 

42000 - 4)UUJ. 

22100  0004 

22200  OOOS 

_ 2X3  oo - 0uQ6_ 

424UO  000? 

22*00  OOO# 

_ 2250* - 0003- 

22*07  0010 

22*00  0011 

_ 22*08. _ 00X2- 

22*10  0013 

22*1*  0014 

— 32&2u _ OOXA. 

22*2*  0016 

22* JO  0017 


22* J2  0010  400 

22* J*  0020 

32Sj(i _ 0021 _ JOtL 

22600  0022 
227uo  002 J 

_ 22000 _ 0024 _ 

22000  0025 

23000  002b 

21P1Q . QQ27 _ 

2 JlOo  0028 

23 JOo  0020 

_2JJ10 _ 0J13.Q _ 

.23400  0031 

23410  0032 

_ 2J*0u _ OOJJ _ 

2J600  0034  110 

23000  0035  100 


_ 23  SiU _ 00.3.* _ 

2 Jb02  0037  C 

23004  0038  C»*« 

_ 2  JO  Ob _ 00  JS _ e _ 

2J810  0040 

23820  0o4 1 

_ 2LJp  JO _ 0042 _ 

23840  0043 

23000  0044  12C 

_ 2  4  0  On _ 0  04* _ C _ 


241U0 

0046 

24110 

00*7 

_ 24120 _ 

-0.048. 

24130 

0049 

24140 

0050 

2420U 

0051 

24400 

C052 

24/00 

00*3 

_  24710— 

.  00*4. 

24/20 

O0SS 

247  JO 

0056 

_  2i/J5-. 

0057. 

. . suoioutine-suitaca - - - - 

C  this  SUbKOUUNt  IS  TU  CHiCF  THE  SURFACE  AhlC«  1*  COVERED 

C  tlTR  FREE  HATER 

_ C0HHUIi/J'AKA/Jt2i*J(»i,0I^IlHX.  PJIXjI^JlZ - - - — 

COMMON/F'HISI  /  NJ  « ALPHA.  '1  MASS » PC  ,  DIF  ,  Phi  ,SL  .SCOUR  .  Sil  I. 
CUM«UN/PNiS2/Pf  , EV, SP, VOL, PUtPTii (*0)  .area, deptnl*0) 

_ eomon/kbUi/lUSOJ - 

COMMON /RATE 1/A*AP, ARM) 
conit.on/contri/KK .  k  1  ,*n 

_ COMMUN/RAU/RKU,l*Kf«,*MU  MIUU«M- - 

C 

C 4 64  »HEN  pONli  IS  DRIED 

— c - - 

IFIVOL  .Cl.  0.)  CO  TO  300 
Nj«0 

_ waiECtu40.0J_ltJJ  VOL _ _ _ 

400  30«>-Am*,tl2.4) 

CO  10  200 

_ 300 _ COf'UJUlE _ 

su*«o. 

DO  100  jal .  (NS»1 ) 

_ LsNUtlSjU-JJ _ 

N«NINS-J) 

SUM«SUM4 IDtPlhll J-DtPlHlK)  )4AHEA*J*10.*4 (-3.) 

_ lUJL_*iQ--U  _$U*1  »S1>* _ 

IF (VOL  .GT,  SUM)  GO  TO  110 
0ELTAaVUL"5UMl 

_ lrXOtLtA.-.tt.  0.1  DEi,TA»Xfii _ 

DtLTti»OELlA/(APEA4l)4l)*10.4*l-3.)) 

NJ«J 

_ GD_TOJJO _ 

110  SUMlaSUK 

100  continue 

_ h J5J _ 

C 

c»M  «nen  pond  is  toil 

_ e _ 

ltlvol  .le.  sum)  go  to  120 
deltaavol'sonil 

_ oe  1  t.a* de  1 1 a / 1 a r ea*.[j ♦  l )  *10.  »♦  t-3 .  j ) _ 

120  CONTINUE 

C _ 

C *  * 4  DETERMINE  WATER  UEFTH  AUOVE  EACH  SEDIMENT  COMPARTMENT 

C 

_ ni«n.t  Cns4.U*nJJ _ 

AWKRSO, 

•  KH«WNP 

_ DO  2  0_o  _L*  1 ,  nj _ 

J«M(NS4i)-L) 

PDEPTH  H)  =  4(UEPTRll)“0EPTH(K)),DELTU 

_  _£ _ _ _ 

C  ADJUSTING  PHUUlLVSlS  RATE 

C 

_ *KP.t«*Ep42f  B/PdepthU  )  ..-~n _ 
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SUBROUTINE.  VOLUME 

•THIS  SUbKOUHAE  .15  10  tSllMAli  ,IHt..E..Kfct_iUllE_.VULjJ*lL*_ 
conMOh/pH»sk/Ph,K.v,sP,voL,i»t'tPrH(i»y)  »«r«*,oeptnit>0) 
coR»on/«A t*r /n ( SO ) 


comwon/pai  4/n»  *M,<Jt,tlinefpri,t.»'t 
co»«'OP/«A*th/ropr  »ioev ,ro*p 


CALCULATE.  THE.  TOTAL  SEEPAGE  UUklNG  THE  TIME  INTERVAL 


IF(DT.E0.  1.)  GO  TO  200 
0T«DT/7. 


t  »UUU 

0001 

1  Uvu 

0007 

174UU 

OoOJ  . 

SUbHOUUNL  1MTL  _ _  .... - 

- ....... - 

1  7300 

0004 

C 

THlb  6UbHCOnlNE  16  TO  INITIALIZE  lnt  SIMULATION  6161EM 

l  Mug 

0005 

COMMON/CONCH/C.TSS.SSH.CSMbUl.CE'T  150) 

^  i joob 

XuKMOh/Pnl 52  APR ,  LV  ,  SP , VOL*  PDtp lnlbALL*  Ul  1 4*.atl)iliL5JLl _ 

nun 

0007 

co"'B'on/p«t«/n8fn»iOI,'hiii8,pr'i  ,T,wi 

1  Miu 

0008 

coiP«iori/ni8«»/tB. 

1  7SUD 

0009 

c 

17600 

0010 

C*M 

READ  IN  1NU1AL  CONCENTRATIONS 

177UU 

0011 

c 

WMyo 

OOIU 

NIADI 5,1001  f  ,  TKS  ,  JiSp 

r/ooo 

001  3 

100 

FORMAT < 1 OEW ,2 ) 

1  7010 

0014 

C 

17^40 

m)j  S 

cm 

read  lr.  initial  cpncpnr rat i ons  in  sediment 

1  7*30 

00 1  b 

c 

1  79  JS 

0017 

do  200  1st, ns 

17941) 

ooia 

MkAD  ( 5 , l5o )  CS»m,rHT(Il 

17047 

0019 

*8111.(6, 1601  CSN(1),CM1(1) 

17945 

0020 

200 

continue 

1  79?iu 

0021 

15u 

tnrmat l helo.4 ) 

18000 

0022 

c 

1.8100 

0023 

HtAD  IN  INITIAL  VULUME 

18700 

■qiVTH! 

c 

18710 

0025 

c 

VOL  IN  LITER,  AHEA  IN  (CK)*»2 

18700 

0026 

Ft.AD  (  5 , 1  Ou  )  VnL.AHbA 

1 8710 

0027 

tmacvvol 

18400 

0028 

*811X16,110) 

18500 

0029 

no 

KORMATUhl  ,20h  INITIAL  CONDITIONS  ) 

_ 18 600 _  0030 _ 

WRITE  16 . 1201  VUL 

Ib/uo 

186O0 

18900 

0031 

0032 

_ 0033 _ 

120  HJRKATlbR  VOL* ,612.4) 

»HJTt(6,130) 

130  KuPMaI  1 81l  Ctmp  .10h  Ch  IUG/L)  .  1  Oh  SS  IUG/D.10H  Cs  (UG/L)) 

19000 

19100 

_ 192  UO _ 

0034 

00J5 

0036 

WHITE  1 6 , 1 40 )  C.TSS.SSR 

140  KUPMAT110X.3E10.2) 

RtlUHN 

19300 

0037 

END 

»  '•*  U  I 


135ou  0056  C 

iiooo  oobs  C  CONPAkiU  lu  Iht  UNt  btkORk  Iht  PktvlUUS  ui<L 

_ liluo - uVeU - C -  - - 

1J600  Oool  JKDEPTHIMH))  OiPThlN IH-l ) ) )  Ull  1U  I1U 

13900  O0b2  L«Wlk-l) 

J _ 14  ODD - OUoJ - WH-1 3  aUXHJ - - - — - 

14100  0004  N ( k  )*L 

14200  0065  M«h-1 

_ 14400 _ 00 to _ 00  IU  44Q - - - 

14400  0067  120  NU*1)«1 

14O00  0068  130  CONTINUE 

_ 14X00 _ 0063 _ - - 

1 4bu0  0070  410  CONTINUE 

14300  0071  C 

_ 16UU0 _ 0012 _ till _ DEFX6t_MJiaA_JUiD_QELIJ _ _ _ _ 

15100  0073  C 

15200  0074  ALPHA«1.8*10.»*(-3.J*PC»(1.-PH1)/PHI 

_ L5A00 _ 0015 _ wmiUJiiu _ 

15500  0076  300  EOkMAT  ( 1 H  1 , 29H1NPUT5  OF  KINETIC  1HF0RMAI  ION  //  9H  COMPAHT- 

15600  0077  S»4X,5HPAPT.,9X,4H  XI, 8X,  JHPHUTUE-,bX,7hh](DROL-,5x,8HVULATl,- 

_ 16100 _ UT 18 _ S.,6A.,.6ftN.lCKtl.  /  3X» Blj-bHCQtf  K..» , 7Ai.OHQA.I.l.UR.t2.L.ilA-»...41lX&l&j- 

15BU0  0079  $1X) ,  6X#7HIZAHON,8X,bHblUDEG  1 

15900  0080  6RlT£(6«12b)NS,PC,WKG, mKP ,6KH,*iKV»6KB 

_ 16000 _ OMJ _ U5 _ ItiRMATl/  15,3x.9t3X,L10.3U _ 

lolOO  0082  C 

lo200  0083  C*M  kRITE  DEPTH 

_ laion  0084  C _ _ _ 

lb 310  0085  Mrlte(6,500) 

163  ■'O  0086  500  (or«at(lhl(7x,4hcoiiiF'>2x>Shdcptn/2xt5norder>4x,10il3x,2nc6) 

_ L6A0Q _ OJULl _ DtLl-7.1)-.K«l«K5 _ 

16500  0088  MK1TL(6,160)K(N(KJ , DEPTH  In ( X ) ) 

•lbbOO  0089  160  tornat(215,eio.2) 

_ lulOU _ flaifl _ l_XO.COM  II  HUE _ 

16800  0091  RETURN 

16900  0092  END 
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v*  nuo  ouo2  suohOgiihI  1NPU7 

Go7tu-  .-UUuJ _ .  C  .  THIS  SutiPuUlinE  IS  1U  ULtlM  SirtiLLAJllhl  >AP  A^tlERS 

O8600  0004  CONPUI./PANA/PS,  N«  ,UT  ,TlNt  ,PK7  ,  1  ,»7 

08900  «00b  CO‘*ML>|,/|<Alt/iiKOf»KP,kftn,aAV,i<lki: 


_ mtumi 

U9100 
091  10 

MUilh _ 

0007 

0008 

Oiuo 

COMMON /PH  J  52  IP*”  ,  t V  , sP,  90L, PUEPHI >1  ( 50 1 , *f  K ,  OCptlH 50  ) 

coR*sii/«*te(/ntSO) 

rilHMOh/ilal  f  l/A«NPrA»Nh 

G91  JO 
0V1  40 

0010 

oon 

c 

co»Mon/4olul/rK(»m,H 

co>raion/eontri/AA>Kl.kn 

0  94uu 

uoi  1 

C«M 

ki.AU  TllLt  1 N  P  OMMAT  J  Oh 

095U0 

0U1  4 

c 

oymu _ 

llil|5 

*KAfH5r«4(Jl  TULL  _ 

oVolu 

0016 

«rltr(b,90j  title 

097  00 

U01  7 

90 

roPHAt C9A10) 

001  fi 

C 

oVVvO 

001  9 

c*«« 

READ  PUN  CONTROL  DATA 

louoo 

OOJO 

c 

1  L  1  UO 

LEAD!  5 , 100  1  NS,1 lMt,DT.  PRT  .m  t  .JCA.A.  1  .  An 

10200 

100 

lOrnatdS. 4(5.1.315) 

10210 

Ou2  J 

■  r  1 1* (6 ,  loo  )  n* , t !■« ,dt ,pr t , »t 

10  400 

^VTTiTT^H 

c 

10400 

0025 

READ  IN  SETTLING  AND  SCOURING  NATL 

lObOo 

002* 

RKADIS.lSOl  SETL. SCOUR 

10700 

002* 

150 

FORMAT ( 2£ 1 0, J ) 

loeoo 

0029 

C 

109tl<l 

_0VAC _ 

c  •  »  • 

READ  IN  SYSTEM  PARAMETERS 

1 10UO 

OuJl 

c 

lllOo 

0032 

HEAO(5,2oO  JPC,  phi,  dll,  SO,  *nD,  r.KP,»NH,  »AV,«Nb 

11110 

vOJi 

r*«alS.200)xk,»«.  ,tA 

11200 

0044 

200 

FORMAl (9f  0,4) 

11210 

00J5 

A»N6*«PB 

luoo 

00  Jo 

c 

11400 

00J7 

c*»« 

REAL  IN  CEPTM  OATA 

ilboo 

00  JO 

c 

11000 

oijiy 

DU  490  J*1 j  NS 

11700 

0040 

490 

HEAU(S,400)  DEPTH ( J ) 

11000 

0041 

400 

PUPHAT (8E12*47 

_ uvoil 

0042 _ 

C 

UOoo 

004  J 

CM» 

ARRANGE  COMPARTMENTS  IN  ACCORDANCE  mIIH  DEPTH 

12100 

00«4 

c 

1  22oO 

0045 

P*1 

12400 

9046 

N(M«1 

12400 

0047 

N(M1 J«2 

12bOO 

0048 

DO  410  I«2,*5 

1  2ovu 

0049 

c 

12700 

0050 

c 

compared  to  the  previous  depth 

_ 12*00 

c 

12000 

0052 

450 

IF1DEPTHU)  .CT.  UEPTHINUD)  SO  TO  120 

1  JoOO 

005  J 

L«N(R) 

0054 

NU  7*1 

1  J200 

0055 

h  (  P  ♦  1 ) *L 

moo 

0056 

M*P 

_ L  140.0 _ 

0057 

HO. 

1KM  .EO.  n  GU  TU  110 _ 
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UBlHiS  0  11S  U'(A^UU(I,7.)  .tit.,  u.)  uu  lb  1«U 

ubuIu  ulle  C 

ilaluu _ UJ-U _ C _  SlOHt  1H1.  UALA _ 

ueuio  Ullt)  C 

0*040  0 1  1 V  Mltxt 


o  *2))  oiio  rROG«o. 

o*2)S  oi)i  ou  10  1*1  $  too 


Q62\Q 

ttu-iA  1 

01)) 

01)4 

JUS 

IfUKUG  .LI.  20.)  mUCt20. 

11*11-1 

CALL  ClkllL(O) 

002M 

01)6 

CALL  GXAPH4(X(1 ) ,1(1 >.11,1 ,)00.,0, ,10. 

,0..1Tu),JA(D) 

00201 

01)7 

C 

_ 0*071 

01)0 

c 

Ti-RHINATIUH  Ilf  GRAPH 

00201 

01)9 

c 

002*1 

0140 

CALL  0PAPH4U(!),U1  ,0.,O.,O.  ,0. 

,11(1 ) , 1A1 1  )) 

0*20) 

0141 

CALL  GPLUT  10. r  0  . r  V99  ) 

uoioo 

0142 

(top 

00400 

014) 

tHD 

1  SO 


Uiiuu _ oubu. 

y  4  JuO  OOtol 
0**4yo  0062 

— ilAiJifl _ UUbX- 

046UO  0004 

y4?vO  0005 

_ u«ttlU) _ Oii-Oo- 

04*00  0067 

049)0  0068 

.iJQvrt 


04940  0070 

Oiuuo  vo7i 

05200  007  J 

05400  0074 

05700  U07o 
0‘jfeuO  0g7  7 
(IVaU.l  11078 


osvio  og79 
uouoo  oogo 
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ueug  oo82 

ucWli  008  J 


06500  008 J 
geoog  0086 
Ua7uQ  0087 
06800  0088 
Og»oO  0089 
i]  7 li u ii  0090 


07100  009J 

07200  0092 


07204  0094 

07206  0095 

-O  J2XO _ 009  a  . 

07220  0097 

07240  0091 


CALL  4  Ur  1 1« FIJI  _ _ _ 

CALL  INPUT 

_UI1Ul2»D1/.Z. _ _ _ 

CALL  FOR  INITIAL  CONCH.  TLAllONS 

call  mm 

•■l . 

nail 


ii«J 

T»T«DT 

lflT  .Cl.  HMD  GU  TO  110 _ 

T2«T* .00000001 
TPPlNT«,f ALSE. 

ill  AHoD  l  T2-.  PftT  J  .LT.  .000001  )  IPRiNTi.IRUE. _ 

H  LAO  IN  HAILS  Of  EVAPUKA110N,  SLLPAOL,  Ant  PHLC IP  1 1  AT  ION 


IF  IT  .NL.  >.  .AND.  AHUD(T2.nT)  .GT.  ,000001  7  GO  TO  160 
REAUlS.lSO)  EV,SP,PH 


H1  — ■  ^VJ 


determination  ot  tcpuring  *na  settling 


HAT10«PH*ARLA«NS/VOL* 10, »•(.),) 

IF1HAT10  .IT.  0.1)  CO  TO  1«0 

_SCR«  LL.lMllUrfi.*lJ  •SCflUi! _ 

SCLXUKAJ  JO-0.1)*SLTL 
CONllNUt 


CALL  FOP  VOLUME 
CALL  VOLUME 

cnccK  to  64*  tnt  pond  it  drlto  out  or  not? 

JU VOL  _ _ _ 

VR11 t ( 6 , 350 )T 

f  0RHA1 ( 26h  THt  POND  IS  PRIED  0U1  AT  ,F5.),5H  DAIS) 
11 


CONI INUE 


07500  0101 

0)600  0104 

yium _ 9ius _ 

07600  OlOe 
O/VOO  0107 

_u.mii _ VXU8 _ 

07920  0109 

07940  0110 

_v?Vlfi _ OLU  .. 

0/950  0112 

0/955  0114 

_U  ttVUU _ 0.U.1 _ 


CALL  COMPS 

PRINT  001  1  HE  SIMULATION  RESULTS 

_ LLlyAL_.9U_41»XJ0_lP_l.l.e _ 

*rlt«lb,200) 

200  ior*4t(22n  tn«  pond  It  drltd  oul) 

_ 90_lo_Ji2 _ 

2lo  conllnu* 

Ulprt  .«q.  1.)  tpr intc.tru*. 

_ IU 7 Pfil_NJJ_6 Hi TE ( 6 ,  170)  T.V OL, AREA , C_ 
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